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ABSTRACT 

The  Bom  approximation  and  the  Schwinger  variational 
principle  are  applied  to  two  one-dimensional  three-body 
scattering  problems.     These  problems  are  also  solved 
exactlj^  and  their  solutions  exhibit  the  physically  in- 
teresting processes  of  elastic,  inelastic,  and  exchange 
scattering.     Comparison  between  the  exact  cross  sections 
for  these  various  processes  with  those  obtained  from  the 
approximation  methods  serves  to  evaluate  the  range  and 
validity  of  these  metbjads.     The  results  are  then  exa- 
mined to  determine  which  are  of  general  applicability 
and  which  are  peciiliar  to  the  specific  problems  studied. 
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Chapter  I 

IITTEODUCTIOH 

In  order  to  describe  the  TDhenomena  which  occur  in  the  ionosphere, 
in  gaseous  discharges,  in  collisions  between  atoms  or  between  electrons 
and  atoms,  a  knowledge  of  the  cross  sections  of  various  physical  processes 
is  reouired.  Experimental  determinations  of  these  quantities  is  beset  with 
many  difficulties.   In  a  large  number  of  cases  it  is  not  possible  to  re- 
produce in  the  laboratory  the  conditions  for  which  the  values  of  the  cross 
sections  are  desired.   In  other  cases,  the  measurements  made  are  a  sum  of 
the  effects  of  different  processes.   Therefore  it  is  important  to  extend 
the  theoretical  methods  for  calculating  the  individual  cross  sections. 

No  exact  solution  has  as  yet  been  obtained  for  a  physical  problem 
which  involves  many  particles.   The  mathematical  complexity  of  such  prob- 
lems reouires  the  use  of  approximating  procedures.   These  in  turn  wust   be 
studied  to  determine  the  reliability  of  the  results  obtained. 

For  two  body  collision  problems,  Rutherford  scattering   for  example, 
an  exact  solution  is  known  and  the  results  obtained  from  the  application 

of  approximating  procedures  caja  be  compared  with  the  exact  results.  Thus 

2 
in  this  case,  as  is  well  known,  the  Born  approximation   yields  the  exact 

scattering  cross  sections.   In  the  absence  of  an  exact  solution  to  a  prob- 
lem, the  evaluation  of  an  approximation  method  has  been  carried  out  by 
comparing  the  results  so  obtained  with  the  available  experimental  data  or 
with  the  results  determined  by  a  numerical  integration  of  the  differential 
equations  involved. 


1.  N.  P.  Mott  and  H.  S.  W.  Massey,  The  Theory  of  Atomic  Collisions. 
(Clarendon  Press,  Oxford,  19^9).   Second  Edition,  Chapter  3,  Sections 
1.  2. 

2.  Ibid.,  Chapter  7,  Section  1. 
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3 
Bates,  Fundaminsky ,  Leech  and  Massey   consider  the  Born  and 

Oppenheimer  approiiniations  and  compare  the  results  with  experimental 

values.   Jost  and  Pais   study  the  Bom  approximation  for  a  Yukawa 

potential  and  compare  the  results  with  numerical  calculations.   Dalitz-' 

on  the  other  hand  ainalyzes  this  same  problem  in  the  limit  as  the  Yu]<Hwa 

potential  approaches  the  Coulomb  potential  so  that  a  comparison  can  be 

made  to  the  well  known  Bom  and  exact  results  for  the  latter  potential. 

Variational  procedures  which  have  been  especially  fruitful   have 

7  8 
been  similarly  evaluated.   In  particular,  the  Schwin^er   *   variational 

principle  has  received  a  great  deal  of  attention  in  the  literature.   Its 

Q 

value  has  been  further  enhanced  by  the  work  of  Kato  who  established  upper 
and  lower  boxinds  for  this  principle  for  two  body  problems.   The  extension 
of  this  variational  principle  to  three  body  problems  was  made  by  Borowitz 
and  Friedman   .   The  Sciiwinger  variational  principle  is  formulated  in 
terms  of  a  Green's  function  which  in  practice  leads  to  computational  diffi- 
culties. A  variational  principle  which  is  more  amenable  to  calculation  is 

1 1  12 

the  Hulthfen   procedure  as  modified  by  Kohn   .   For  two  body  problems  Kohn 

(see  footnotes  12,13)  was  able  to  establish  a  relationship  betveen  his  pro- 

14 
cedure  and  the  Schwinger  variational  principle.  Kato   has  been  able  to 


3.  D.  E.  Bates,  A.  Fundaminsky,  J.  W.  Leech  and  E.  S.  rf.  Massey,  Phil. 
Trans.  Roy.  Soc.  A2^3.  93  ( I950) . 

4.  R.  Jost  and  A.  Pais,  Phys .  Rev.  82.  840  (1951). 

5.  R.  H.  Dalitz,  Proc .  Roy.  Soc.  A  206.  509  (1950). 

6.  J.  M.  Blatt  and  J.  D.  Jackson.  Phys.  Rev.  26.  Ic  ( 194-9)  . 

7.  J.  Schvvinger ,  Hectographed  Notes  on  Nviclear  Physics,  Harvard  194? 
(unpublisjied)  . 

8.  N.  Marcuvitz,  Sec.  HID.  Recent  Developments  in  the  Theory  of  Wave 
Propagation.  New  York  University,  Institute  of  Mathematical  Sciences 
(1949). 

9.  T.  Kato,  Prog,  of  Th.  Phys.  6,  394  (I95l)- 

10.  S.  Borowitz  and  B.  Friedman,  Phys.  Rev.  82,  441  (1953)- 

11.  L.  Hulthen,  Kungl.  Pysio,  Sallskapets  Lund  Forhand  14.  21  (1944). 

12.  W.  Kohn,  Phys.  Rev.  Jk,   1763  (1^48). 

13.  W.  Kohn,  Phys.  Rev.  84,  495  (1951)^ 

14.  T.  Kato,  Phys.  Kev.  80,  475,  (1950). 
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exhibit  the  relations  among  the  Schwinger,  Kohn  and  Hulthen  variational 
principles  for  two  hody  problems. 

Altshuler^  has  applied  the  Schwinger  variational  principle  to  the 
study  of  the  static  field  model  of  the  hydrogen  atom  and  compares  the 

results  with  the  numerical  calculations  of  Chandrasekhar  and  Breen   and 

17 
the  estimates  of  Massey  and  Moiseiwitsch  .   Some  of  the  difficulties  and 

practical  limitations  which  arise  in  a  study  of  the  application  of  varia- 
tional principles  to  three  body  problems  can  be  appreciated  by  considering 
the  work  of  Massey  and  Moiseiwitsch.   These  authors  compare  the  results 

-1  Q 

obtained  from  the  Hulth^n,  Kohn  and  Huang   variational  principles  with  the 

19  20 

numerical  calculations  of  Macdougall   and  of  Morse  and  Allis   .  Because 

numerical  results  are  not  available,  however,  the  effectiveness  ol  the 

variational  methods  in  dealing  with  exchange  and  polarization  cannot  be 

21  22 

explored.  The  results  obtained  by  Verde  ~  and  by  Troesch  and  Verde   -using 

Hulthln's  variational  principle  are  similarly  restricted. 

The  approach  to  the  problem  of  evaluating  various  approximation  methods 

adopted  here  is  to  consider  three  body  problems  which  can  be  solved  exactly. 

Then  the  same  problems  can  be  treated  by  various  methods  of  approxima.tion 

and  the  results  obtained  in  this  way  can  be  compared  with  the  exact  results. 

Wildermuth  -^»^^  has  constructed  and  solved  a  class  of  three  body  problems 

which  am   be  used  for  this  ptirjiose.   These  problems  are  one-dimensional 


15.  S.  Altshuler,  Phys.  Rev.  82,  1278  (1953). 

16.  S.    Chandraseldaar  and  F.  Breen,  Astrophys.  J.  10^,  4l  (19^6). 

17.  H.  S.  W.  Massey  and  B.  I.  Moiseiwitsch,  Proc.  Boy.  Soc.  A  205. 
483  (1951). 

18.  S.  S.  Huang,  Phys.  Hev.  2^,   ^77  (19^9)- 

19.  J.  Macdougall,  Proc.  Eoy.  Soc.  A  I36 ,  5^9  (1932). 

20.  P.  M.  Morse  and  W.  P.  Allis,  Phys.  Eev.  44,  269  (1933). 

21.  M.  Verde.  Helv.  Phys.  Acta,  22,  339  (19^9). 

22.  A.  Troesch  and  M.  Verde,  Kelv.  Phys.  Acta  24,  39  (1951)' 

23.  K.  Wildermuth,  Research  Report  CX-3 ,  I4ath.  Res.  Gr.,  Hew  York  Univer- 
sity (1952). 

24.  K.  Wildermuth,  Z.  f.  Phys.  12^,  92  (1949). 
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three  body  problems  which  exhibit  all  the  essential  physical  phenomena 
such  as  elastic,  exchange  and  Inelastic  scattering  that  occur  in  actijal 
physical  problems.   These  problems  are  mathematically  tractable  so  that 
the  physical  situations  can  be  explored  in  great  detail.  They  also  serve 
as  a  means  of  gaining  some  insight  into  real  physical  problems  which  are 
much  more  comDlex.   One  of  the  main  goals  of  this  study  is  to  examine 
the  results  obtained  to  see  which  are  of  general  validity  and  which  are 
consequences  of  the  specific  problems  stTodied. 

Two  problems  will  be  considered.   The  first  of  these  is  a  three  body 
one-dimensional  problem  which  involves  distinguishable  particles.   Its 
exact  solution  the  calculation  of  the  various  scattering  cross  sections 
and  the  application  of  the  Bom  approximation,  the  Hulth^n-Eohn  and 
Schwinger  variational  principles  make  up  the  content  of  Chapter  H.  The 
second  problem  is  again  a  three  body  one-dimensional  problem  but  now 
identical  particles  are  involved.   In  addition  to  elastic  and  inelastic 
scattering  which  also  occur  in  the  first  problem  exchange  scattering  is 
possible.  The  overall  procedure  that  is  followed  for  this  problem  is  the 
same  as  for  the  preceding  one.  However,  the  methods  of  approximation  can 
be  applied  in  two  distinct  ways.   The  first  is  to  take  the  pertxirbation  in 
a  symmetric  fonn.   The  development  of  the  approximations  in  this  way  is 
dealt  with  in  Chapter  III.   The  other  approach  is  to  take  the  perturbation 
in  an  asymmetric  form.  The  application  of  this  procedure  which  is  the  one 
most  often  fo-und  in  the  literature  is  the  content  of  Chapter  IV. 

The  extensive  use  of  the  Born  approximation,  especially  in  many  particle 
problems,  has  led  us  to  examine  this  procedure  in  some  detail.  Specifically, 
the  Born  expansion,  which  results  when  the  procedure  for  obtaining  the  Bom 
approximation  is  itert.ted,  is  studied.  An  important  result  is  that  the  Bom 
expansion  is  shown  to  converge  whereas  the  Bom  approximation  may  have  no 
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useful  domain  of  applicability  whatever.  The  second  Bom  approximation 
which  is  a  substantial  correction  to  the  Bom  approxiination  does  however 
have  a  useful  domain  of  applicability.  Jost  and  Paia^  obtain  results 
similar  in  content  to  these  but  for  lack  of  numerical  results  are  unable 
to  evaluate  the  significance  of  the  contribution  from  the  second  Bom 
approximation.  The  magnitude  of  the  correction  due  to  the  second  Bom 
approximation  is  indicative  of  the  rapidity  of  the  convergence  of  the 
Bom  exjjansion.  A  similar  result  can  be  fo^und  in  the  work  of  Dalitz 
who  evaluated  the  Bom  expansion  up  to  the  third  Born  approximation.  The 
calculation  was  carried  out  to  show  the  convergence  of  the  Bom  expansion 
but  this  coiild  also  be  seen  from  the  calculation  of  the  second  Bom  approxi- 
mation. 

The  above  resiilt  implies  that  in  practice  additional  terms  in  the  Bom 
expansion  must  be  calculated  if  useful  res\ilts  are  to  be  obtained.  Actually 
such  calculations  are  not  feasible  and  other  methods,  in  particxilar,  varia- 
tional procedures  are  used. 

The  Schwinger  and  Hulthen-Kohn  variational  principles  are  applied  to 
the  two  problems  considered  here.   Trial  fomctions  are  restricted  to  be  those 
of  the  Bom  type  because  in  general  the  Schwinger  variational  principle  is 
mathematically  tractable  for  such  trial  functions.  However,  in  this  case 

the  Hulth^n-Zohn  variational  principles  yields  the  same  results  as  the 

27 

Bom  approximation  .  Nonetheless  the  study  of  this  variational  principle 

leads  to  a  generalization  to  many  body  problems  of  the  relation  between  this 

2ft 
and  Schwinger 's  variational  principle  which  was  first  given  by  Kohn   for 

two  body  problems. 


25.  loc.  cit. 

26.  loc.  cit. 

27.  Kohn,  loc.  cit. 

28.  ibid. 
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In  Chapters  II  and  III,  the  Schwin^er  variational  principle  is  applied 
to  the  problem  of  distin^juisha'ble  particles  and  to  the  prohlem  involving 
identical  xjarticles.  For  Born  trial  functions,  this  principle  was  expected 
to  yield  results  that  included  the  second  Born  approximation^°»30.  However, 
the  problem  treated  in  Chapter  II  shows  that  the  variational  calculations 
do  not,  even  for  very  large  energies,  approach  the  exact  or  Born  results. 
On  the  other  band  for  smaller  values  of  the  energy  the  variational  results 
are  excellent  where  the  Eom  results  are  completely  inadequate.  The  result, 
for  the  symmetric  problem  dealt  with  in  Chapter  III,  is  very  interesting. 
Here  both  the  Schwinger  and  Born  results  are  good  approximatic-is  to  the 
exact  values  for  the  antisymmetric  solution  of  the  problem  for  large  values 
of  the  energy.   But  now  for  smaller  values  of  the  energy  the  Bom  results 
are  very  good  whereas  the  variational  results  become  far  less  satisfactory. 
To  summarize  the  situation,  these  two  problems  indicate  that  if  Bom  type 
trial  functions  are  used  in  the  Schwinger  variational  principle,  the  results 
which  are  obtained  may  be  eo^uivalent  to  the  second  Born  approximation  in  that 
range  of  the  energy  where  it  is  a  good  approximation  to  the  exact  results 
but  tiiat  outside  this  energy  range  the  variational  principle  fails  to  extend 
the  domain  of  usefulness  of  the  Born  results.   Furthermore  the  variational 
and  Bom  results  may  not  have  a  common  domain  of  agreement  yet  one  may  give 
excellent  results  Just  where  the  other  is  unreliable. 

Chapter  IV  is  a  study  of  the  symmetric  problem  formulated  with  an 
asymmetric  perturbation.   The  main  consideration  here  is  to  see  whether  or 
not  exchange  scattering  can  be  proJ)erly  treated  by  this  approach  to  the 
problem.  An  iterative  procedure — the  Bom  expansion — is  examined  in  detail 
and  it  is  found  that  this  method  does  not  yield  exchange  effects.  Only  if 


29.  Schwinger.  loc.  clt. 

30.  Altshuler,  loc.  cit. 
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the  resiilting  series  is  evaluated  to  ottain  the  analytic  function  it 
represents  in  its  restricted  domain  of  validity  can  exchange  phenomena  appear. 
Hence  the  Bom  expansion  as  an  approximation  method  fails  to  yield  any 
information  about  exchange  scattering.  This  result  is  in  contrast  with 
that  obtained  by  the  same  method  for  the  symmetric  pert\arbation  formula- 
tion of  the  same  problem.   There  each  term  in  the  Born  expansion  displayed 
the  effects  due  to  exchange. 

A  procedure  due  to  Mott  and  Massey-^  which  was  devised  to  take  exchange 
scattering  properly  into  account  was  also  studied.   It  was  found  that  if 
the  initial  state  of  the  system  is  taken  as  a  trial  function,  which  is  the 
procedure  followed  in  the  Bom  approximation,  then  spurious  terms  arise 
in  the  calculation  of  the  exchange  scattered  wave.   Since  this  procedure  is 
an  approximate  one  such  terms  yield  contributions  which  assume  unwarranted 
importance  especially  if  the  results  are  to  be  used  in  an  iterative  procedure. 
Hov/ever,  if  trial  functions  are  used  which  represent  the  initial  state  of 
the  system  when  the  interaction  between  the  incident  and  bound  particles  is 
neglected  then  the  spiirious  terms  do  not  appear  and  the  resulting  approxi- 
mation is  also  improved.  These  results  are  directly  applicable  to  actual 
physical  problems. 

Chapter  V  is  a  discussion  of  the  results  obtained  from  the  preceding 
sections.  Its  object  is  to  point  up  those  results  which  have  general 
validity  and  those  whose  special  character  indicate  the  direction  in  which 
further  study  would  be  profitable. 


31.   Cf.  Reference  1,  Chapter  8,  Section  2,  4. 


Chapter  II 

DISTIIIGUISHABLE  PAETICLES 

2.1  Statement  of  the  Pro'blem 

We  shall  consider  a  system  composed  of  an  "atom"  and  aji  incident 
particle.   The  "atom"  consists  of  a  nucleus  with  infinite  mass  and  a 
"boimd  particle  with  finite  mass.   In  addition  it  will  he  assumed  that 
the  system  is  one- dimensional  so  that  each  of  its  constituents  can  he 
descrihed  in  terms  of  a  single  coordinate. 

The  time  independent  Schroedinger  equation  for  this  system  is 

/-  -i!  A!-iL^  '  ZBbix^)  *nx^x^)\     Kx-^T^)  '   SKx^x^) .         (2.1.1) 
I  2m-j^  3x|   2m2  S^|  J 

The  time  dependent  factor  has  heen  taken  to  be  exp(-iEt/li) .  The  poten- 
tial energy  of  the  hound  particle  is  -2B  5(12)  wherr  B  is  a  real  positive 
constant  and  5(xp)  is  Lirac's  delta  function.  Vd,!^)  represents  the  energy 
of  interaction  between  the  incident  particle  1  and  the  bound  particle  2. 
Its  explicit  form  will  be  given  later.   m,  and  m2  are  the  masses  of 
particles  1  and  2  respectively  and  -if  is  Planck's  constant  divided  by  2n. 
The  total  energy  of  the  system  is  E.  Particle  1,  the  incident  particle, 
does  not  interact  with  the  nucleus — the  center  of  force — and  is  thereby 
distinguished  from  particle  2.  In  wliat  follows,  units  are  chosen  sxich 
that  -^  =  1  and  m^^  =  m2  ■  1/2. 

The  unpert\irbed  eq-uation  which  is  obtained  from  (2.1.1)  by  setting 

V^XiXj)  eoual  to  zero  is  serierable.  As  will  be  shown  in  section  (2.2), 

2 

particle  2  has  but  one  botmd  state  with  energy  -B  and  all  other  states 

are  in  the  continuxom  with  positive  energies.   This  is  a  consequence  of  the 
Dirac  delta  type  potential.  Furthermore  the  system  is  capable  of  exhibiting 


-  9  - 

both  elastic  and  inelastic  scattering.   The  latter  is  particularly  simple 
and  represents  the  ionization  of  the  "atom."  The  short  range  of  the  delta 
function  has  as  its  co-unterpart  in  actual  physical  prohlems  the  short  range 
nuclear  forces. 

Equation  (2.1.1)  will  he  considered  in  momentiun  space.   The  eidvantage 
of  the  momentum  over  the  coordinate  representation  was  first  pointed  out 
lay  Dirac-^  and  the  subseouent  work  of  others33»3^«35  has  demonstrated  its 
usefulness  for  the  consideration  of  scattering  problems.   To  obtain  the 
momentiam  representation  of  eq-uation  (2.1.1)  let 

I     >  00 

^^Vz^  ■  ^^"  ^~    /  I  ^^\^Z^    eM^\y^^  i.\x^)^^^   .      (2.1.2) 

Substitute  this  into  (2.1.1)   and  m\iltiply  the  resulting  equation  by 

exp  (-ik'x     -  ik' X  ) .       Then  integrate  over  -co  <  x,  <  +  oo     and  -  oo  <  x«<*oo. 

Eqiiation  (2.1.1)    becomes 

(2.1.3) 

(k2*  k^  -  E)f(k^k^)    -      ■?/^f(  V2^<^2'^];^/y(^^2l^l^>P<^^^  ^(44^=  ° 


where 


(k^k^lvlk^kp  »y    yv(x^X2)expj-i(k^-k|)x^  -i(k^-k^)x2^dx^dx^     .        (2.1.4) 


The  identity  (2.1.4)  will  henceforth  be  referred  to  as  the  interaction 
between  the  particles. 

To  prevent  the  notation  from  becoming  too  cumbersome  v;e  shall  hence- 


32.  P. A.M.  Dirac ,  Quant^om  Mechanics.  Third  Edition  (Clarendon  Press,  Oxford, 
1947) . 

33.  W.  Kohn,  Phys.  Hev.  "Jk'   ^763  (19^8). 

34.  K.  Wildermuth,   cf.  Reference  23. 

35.  W.   Heisenberg,   Z.f.  Phys.   120,  513  (1943). 
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forth  use  a  single  integral  sign  to  replace  multiple  ones.  The  multi- 
plicity of  integrations  to  te  performed  will  be  indicated  by  the  n\imher 
of  variables  of  integration.  The  range  of  integration  will  always  extend 
from  -00  to  foo  for  all  the  variables  of  integration.  Any  deviation  from 
this  convention  will  be  explicitly  indicated. 

To  complete  the  nathematical  description  of  the  physical  system,  the 
interaction  between  the  jDarticles  must  be  specified.  This  interaction  will 
be  chosen  to  be  herraitian  so  that  the  eigenvalues  of  the  Schroedinger 
equation  are  real  and  the  form  of  the  eauation  for  the  conservation  of 
probability  density  will  be  preserved.  The  interaction  assucod  is 

(k  k„|V|k'k')  =  aL 1 ^  +   ., 1 ^  V  .   (2.1.5) 

'HL  2'  '^  2'  -    |(k^  +iT)(k^  -1^     (k^*  iT)(k|  -  iT)  j     ' 

where  A  and  T  are  real  constants.   It  follows  from  (2.1.5)  that 

(k^k^lTlkjkp  =  (k^k^lTlk^k^)*  .  (2.1.6) 

where  the  asterisk  denotes  the  complex  conjugate  of  the  quantity  so 
modified. 

The  form  of  the  interaction  (2.1.5)  is  dictated  by  the  following 
considerations:   Equation  (2.1.3)  is  readily  solved  if  the  last  term  on 
the  left  hand  side  can  be  evalxjated.   This  can  be  done  by  taking  the  inter- 
action  to  be  a  product  of  the  form  V(k.k  )  V(k'k')   .  Another  Important 
property  of  the  interaction  is  that  when  expressed  in  coordinate  space  it 
should  vanish  rapidly  for  large  separations  between  the  "atom"  and  the 
incident  particle.  This  is  satisfied  by  (2.1.5)  •   In  coordinate  space 
this  interaction  is  an  integral  operator  wit/i  a  kernel  given  by 
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a[5  (ig)  ^(-Xj^)  exp{']!x{)  5  (7:{)    ^(-zp  exp(  Txp 

where  /  0     x  <  0 

^         VI     X  >  0 


This  interaction  is  evidently  one  which  vanishes  rapidly  for  large  distances 
X,  and/or  x_  *  Physically,  the  two  particles  interact  only  when  one  of  them 
is  at  the  origin  and  the  other  is  to  the  left  of  the  origin.   The  interaction 
"between  them  is  greatest  when  the  particles  are  at  the  origin. 

The  interaction  (2.1.5)  may  "be  further  simplified  to  a  produc*  of  two 
factors.  However,  in  that  case,  the  Schwinger  variational  principle  gives 
the  exact  solution  for  arbitrary  trial  functions.  To  make  this  principle 
depend  on  the  choice  of  trial  functions  a  sTom  of  such  products  must  "be  chosen. 

The  physical  problem  which  is  mathematically  represented  hy  equation 
(2.1.3)  and  the  identity  (2.1.5)  is  the  following:  Particle  1  is  incident 
from  the  left  with  momentum  K^  and  interacts  with  the  "atom. "  Since  it  does 
not  interact  with  the  center  of  force,  exchange  scattering  is  impossible, 
but  it  does  undergo  elastic  and  inelastic  scattering.  The  exact  wave  func- 
tion and  the  cross  sections  for  elastic  and  inelastic  scattering  are  to  be 
determined. 

These  cross  sections  are  also  to  be  calculated  from  the  Bom  approxi- 
mation, the  Hulthen-Kohn  variational  principle  and  the  Schwinger  variational 
principle.  The  comparison  of  these  results  with  the  exact  result  will  enable 
us  to  evaluate  these  approximating  procedures. 

2.2  Mathematical  Preliminaries 

Equation  (2.1. 3)  with  the  interaction  (2.1.5)  can  be  solved  in  various 
ways.  For  purposes  of  reference,  the  methods  of  solution  will  be  designated: 
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the  direct  method,  the  expansion  method  and  the  Green's  function  techniaue. 

All  these  methods  of  course  lead  to  the  same  result  hut  the  details  of  each 

method  will  he  fo\ind  useful  when  the  various  methods  of  approximation  are 
considered. 

2.2.1  The  Direct  Method 

Initially  particle  1  and  the  atom  are  so  far  apart  that  the  interaction 
■between  them  is  negligible.  The  initial  state  of  the  system  is  then  a  solu- 
tion of  tlie  unperturbed  Schroedinger  equation 

(k^+  ^  -  2)  ^^W^   -   (S/")  f^o^\^   ^^2  "  °  •  (2.2.1) 

which  is  obtained  from  eajuation  (2.1.3)  by  setting  the  interaction  equal  to 
zero.   This  equation  when  transcribed  to  coordinate  space  is  separable  and 
its  solution  is  straightforward.   It  is  also  separable  in  momentum  smce  and 
its  solution  may  be  readily  obtained  by  exploiting  the  physical  situation 
which  it  represents. 

Since  particle  1  is  incident  from  the  left  with  momentum  K  and  does 

0 

not  interact  with  the  center  of  force  or  with  the  bound  particle  2  it  m\ist 

proceed  undisturbed.   The  state  of  particle  1  is  described  by  a  plane  wave 

-1/2 
(2  n)   '   exp(  iK  I, )  in  coordinate  sr>ace  which  is  represented  by  6(k,-Z  ) 

in  momentum  space.   The  normalization  of  the  plane  wave  was  chosen  as  in- 
dicated purely  as  a  matter  of  convenience.   The  inonent''jm  representation  is 
easily  verified  by  substituting  it  into  the  one-dimensional  equivalent  of 
(2.1.2).   Therefore  the  solution  of  (2.2.1)  has  the  form 

'o  ^^1^2)  =   6(k^-X^)f(k2).  (2.2.2) 

where  f(k  )  mrust  represent  the  bound  state  of  particle  2  since  the  latter 
is  undisturbed  by  the  passage  of  particle  1.   The  energy  of  particle  1  is 
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2  2 

K   ,  whence  the  energy  of  particle  2  must  he  I„  ■  3  -  Z      ,  where  E  is  the 

0  /i  0 

total  energy  of  the  system.   Substitute  (2.2.2)  into  (2.2.1).   Then  we 
find  that  f(kj,)  must  satisfy 

•   (k|  -  S^)  nia^)  -   (B/n)D  =  0,  (2.2.3) 


where  the  constant 

D  =    y  f(k2)    dJCg   .  (2.2.4) 

Prom  (2.2.3)  ajid.  the  reauirement  that  f(k  )  represents  a  bound  state, 
it  follows  that 

f(k2)  »  (B/n)  D/(k^  -  E^),  (2.2.5) 

provided  that  S  <  0.  From  (2.2.4)  and  (2.2.5) 


./  i-^]"  . 


The  only  solution  is  that  for  which  E_  =  -  B   ,  where  the  positive  square 
root  must  he  taken  since  B  is  positive.  The  constant  D  which  appears  to 
he  arbitrary  is  fixed  by  the  normalization  of  f(k  ).   Hence  the  solution  of 
the  unperturbed  Schrcedinger  equation  (2.2.1)  has  only  one  bound  state  with 
negative  energy  for  particle  2.  This  bound  state  is  given  by 

f(k2)  a  ll/(k^  +b2).  (2.2.6) 


where 


N  =  (2/n)^/^  B^/^  (2.2.7) 


N  is  a  constant  determined  from  the  normalization  of  f(k_) .  Combining 
equations  (2.2.6)  and  (2.2.2)  we  find  that  the  solution  to  equation  (2.2.1) 
is 

W2^  '    ^^h'  h^   V*^^2  *  ^^^  (2.2.8) 

with 

N  =  (2/n)-^''^  b3/2  and  K^  =  B^+  B  .  (2.2.9) 
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Now  represent  the  solution  of  eaiiatlon  (2.1.3)   ^s 


f(klk2^  =     fo(^^2^    "  ^1^^^2^ 


(2.2.10) 


where  f  (k,k  )  is  the  solution  (2.2.8)  of  the  unperturbed  equation  (2.2.1) 
and  f,(k,kp)  represents  the  outgoing  flux  of  particles  which  results  from 
the  interaction  between  the  incident  and  "bound  particles.   Substitute 
(2.2.10)  into  eo^uation  (2.1.3).  Then  we  find  that  f,(k,k_)  must  satisfj 

(2.2.11) 
(k^  .  k2  -  E  )  f^(k^k2)  -  F(k^>  ^ 


.=  0, 


^ 


+  11 


V 


iT 


in  which 


and 


F(k^)  =  {B/n)  ff^{]s:^lc^)    dk^ 

C  =  y  dk^  dk^  i(^\)  I  '\  -   il)  . 
D  =  f    dkj^  dk^  i^W)  f^\-  ^'^)' 


(2.2.12) 
(2.2.13) 


From  (2.2.11)  we  obtain 


f.dc^k^)  = 


F(k  )  -  -k, 


S^ 


il 


k2  +  iT 


J  ^ 


*^ 


(2.2.14) 


In  this  eq-uation  it  is  to  be  understood  that  a  small  positive  imaginary 
part  (see  footnotes  23,  26,  28)  has  been  added  to  the  total  energy  2  to 
insure  that  f  (k  k  )  represents  outgoing  waves  in  coordinate  snace.  A 
singularity  of  the  delta  function  type  does  not  appear  in  (2.2.14)  because 
the  condition  on  f,(k.,k  )  is  that  it  represents  outgoing  waves  only.  F(k.) 
can  be  evaluated  by  substituting  (2.2.14)  into  (2.2.12).   Then 


^I'^i^a'  =  -  V, 


\J   3-1^  -  Bi  k^  * 


+  D 


Bi 


iT 


+ 


I  Vs-k^  -Bl)(/S-k2  *   iT) 
(2.2.15) 


k^  *  iT 


k2  *  k^  _  E 
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The  constants  C  and  D  are  determined  from  the  evaluation  of  (2.2.13) 
Their  values  are 


=   In  J^  (i-     AjiN   _   _ 

B     ^  B+T     \  Un     •■■  /  kv. 


-^^    (BJ2*  1^3)  \   A  .  (2.2.16) 


D  r     in  J  _1 A   _     _A  J  A     ^        JA "^^    I      A 

B      I   K-^T\  ^T-     ly         /Jn-BflT         ~T~   /      ^  • 


where 


A''  =  (1-  ^Jii)  *  (5IJ     ^(^V^j)'  '^-'-^'^ 

and  Jj^,  (i  ■  1,  2,  3.  ^)  are  fuiictions  of  the  total  E  and  the  constants  B 
and  T.  The  explicit  expressions  for  these  fimctions  will  be  found  in 
Appendix  I. 

From  (2.2.15)  and  (2.2.10)  we  find  that  the  exact  solution  to  equation 
(2.1.3)  for  the  interaction  (2.1. 5)  is 

^  '^  (k2+  b2)      4nZ    k?+k2  -B 

2  12 

-/^  C 


D^  ^^ 


VT^  -  Bi   ki-  iT      \  ^/^  -^"^  ^<^  *  ^^^  '  ^2'  '^ 


2.2.2  The  Srpansion  Method 

The  solution  of  eo^uation  (2.1.3)  can  be  found  by  expanding  it  in 
terms  of  the  eigenfxmctions  of  the  unperturbed  Schroedinger  eouation 

(k2+k2  -  E)  g(k^k2)  -  (B/n)  f  g()^^)^^)dk^   -  0.  (2.2.19) 

The  eigenfunctions  of  this  equation  are  to  be  determined  subject  to  the 
condition  that  they  represent  incident  plus  outgoing  waves.  To  obtain 
these  eigenfunctions,  equation  (2.2.19)  can  be  solved  in  a  manner  analogous 
to  the  solution  of  equation  (2.2.1).   Thus  one  of  the  eigenfunctions  of 
(2.2.19)  is 
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(K^  ?  \)t^^^)   =  e<VK^)  N/(k^+  B^)  (2.2.20) 

2     ? 

with  K  =  B  +  E.  K  is  the  momentum  of  the  incident  particle  1.  The 

other  eigenf-unction  is  obtained  from  equation  (2.2.3)  which  when  written 
in  our  present  notation  becomes 

(k^  -  k2  j  g(j^j  _(B/n)D  „  0,  (2.2.21) 

<-    2      ^ 

where  K^  =  E  -  K^   and 
2        1 

D  =  y gdc^  dkg  .  (2.2.22) 

The  solution  of  (2.2.21)  which  satisfies  the  conditions  of  the  problem 

is  T,  -n 

g(kj  =   6(k2-K2)+  -^ 1 —-  .  (2.2.23) 

2    2 

A  small  positive  imaginary  part  iP<>('^>0)  has  been  added  to  the  total  energy 
S  to  insure  that  the  last  term  on  the  right  hand  side  of  (2.2.23)  represents 
outgoing  waves  in  coordinate  space.  Throughout  the  remainder  of  this 
work  when  outgoing  waves  are  considered  it  will  be  vuiderstood  that  the  total 
energy  E  will  have  associated  with  it  a  small  positive  imaginary  part. 

The  constant  D  is  determined  by  substituting  (2.2.23)  iato  (2.2.22) 
and  carrying  out  the  indicated  integration.   In  this  way  we  find  that  the 
second  eigenfunction  is 

(K.K,  |k.kJ  5  6(krt)/'6(k  -O*  2 l^^i 1 I    (2.2.2^) 

^  ^  ^  ^        1  1  1^   2  2    n    |j^^,  _  g^  kl-K^-  i^J 

with  Z^  +  K^  =  E. 
1    2 

The  eigenfunctions  (2.2.20)  and  (2.2.24)  satisfy  the  orthogonality 
relations 
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/(K|K^  ik^k^V^dk^  (Kj^K^lk^k^)  =  5  (K^  -  K|)  6{K^   -  X^)  .         (2.2.25) 
/(K|2|k^kp*dk^clk2  (Z^  2|k^k2)  =  6(K^  -  K') 

/  (K^Z2  K^^*^l*^2  ^^1  ^  V2^  "  ° 
Ike  closure  property  is  also  satisfied,  i.e., 

/(Z^Kglk^kp^dZ^dZ^CK^K^Ik^k^)  *  y(K^  2|k«_k«)*dK^(K^  2  |  k^k^)     (2.2.26) 

=  5(k^  -  kp  6(k2  -  kp   . 

Equations  (2.2.25)  and  (2.2.26)  assert  that  the  eigenfunctions  (2.2.20) 
and  (2.2.24)  constitute  a  coinplete  set. 
To  solve  equation  (2.1.3)  we  let 

f(k^k2)  -  fo(k;L^2^  +  f3^(k;L^2^.  (2.2.27) 

with  the  same  significance  attached  to  f  (k-j_lc2)  and  f2^(k2^k2)  as  was  given 
to  them  below  eqioation  (2.2.10).  Then  f  (k-k  )  must  satisfy 

(k^+k^  -  E)  i^i\^2^   -   (B/n)yf^(k^k2)dk2  +  V(k^k2)  =0.       (2.2.28) 

in  which 

V(k^k2)  f  (l/4n^)/(l^k2l^l4^p<"^<^2  ^^^^^^  (2.2.29) 

Kow  expand     f^(k  k  )    in  terms  of  the  eigenfunctions  of  the  xmperturbed 
eoTjation.     Then 

V^^^   =/f(l4)dZ^(K^  2|k^k2)*    /f(KiZ2)dZ^dX2(K^K2|k^k2).        (2.2.30) 

where  t{K^)   and  f(Z^Z2)  are  the  expansion  coefficients.  These  are  different 
functions  hut  to  simplify  the  notation  their  arg\aments  are  used  to  distingxiish 
between  them.   Substitute  (2.2.30)  into  (2.2.2.8).   Then  the  latter  eqijation 
becomes 
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+  /(k^  +  k^  -  E)f(K^K2)dX^dK2(KlI^hL^2^   "  ^^/"^  /  f (K^K2)dX^dIi2(K^K2l  k^k2)dk2 

=     -  VCk^kg).  (2.2.31) 


Since 


(k2+  ^     -  ^l'  Kp{K^K^\k^k^)   -  (B/n)  /(Z3^K2|kj^k2)clk2  =  0, 


and 


(k^  +     k^  -  k2+  B2)(K^2|k^k2)   -   (B/n)/  (K^2  |k^k2)dk2  =  0, 

eoToation  (2.2.31)    reduces  to 

/(K^  -  b2  -  B)    f(K^)dK^(K^  Zlk^k^)  (2.2.32) 

+  /  (K^  +    K^  -  E)    f(K  K  )cLK  dX  (K  K    Ik  k  )   =  -  V(k  k  )  . 
•^        1         2  12121212  12 

Miiltiply  this  eqmtion  with  (K«  2  |k  k  )  dk^dk  and  (Z'H'  |k^k  )*dkLdk   in 

turn  and  integrate  over  the  indicated  variables.  Then  the  following  two 

eouations, 

(K^  _  K2)f(K^)  =  -  f{\   2  Ik3^k2)*dk^dk2  Y{k^]^) -,   Z^  =  B^  *  E 

and 

(K^  +  K2  "  ^^  ^^\^2^   =  ■/^V2'^lV*'"^^2  ^^^iV  (2.2.33) 

result  when  the  orthogonality  relations  (2.2.25)  are  used.   Since  f(K,) 
and  f(K^K  )  are  to  represent  outgoing  waves,  it  follows  that 

1  o 
and  (2.2.34) 

1    2 
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When  equations  (2.2,3^^)  are  substituted  into  (2.?.30)  we  find  that 


fidc^kg) 


-/ 


(K^Z2|k^k'T  6K^6K^{K^K^\)^y^)      (K^2|k^k^T  dK^(K^  Z\   k^k^) 


k2  +  k2  -  B 
1    2 


e2-  k2 
1    0 


dkjdk|  V(k|k^) 


(2.2.35) 


The  evaluation  of  the  right  side  of  this  equation  yields  the  same 
expression  for  f-(k  k  )  as  obtained  "by  the  direct  method — equation  (2.2.15). 

2.2.3  The  G-reen' 3,  Function  Technique 

The  Green's  fxmction  for  equation  (2.1.3)  is  that  solution  of 

(k^  +  k^  -  E)G<kj^k2;kjkp-(B/„)y*&(k^k2;  ]c^k^)dii^=     6(^1-4)  b(l^^-k^)  ,     (2.2.36) 


which  represents  outgoing  waves  in  coordinate  space.  This  equation  can 
he  solved  by  either  of  the  previous  methods  discussed.  The  expansion 
method  gives  the  result  immediately,  for  G(k^k  ;  k'k')  satisfies  the  same 
conditions  as  f,(k^k2) .  Thus  to  obtain  the  Green's  function  from  (2.2.35) 
we  need  only  replace  y(k,^k  )  by  -6(k  -  kJ)  6(k  -k')  .  It  follows  that 


G(k^k2:  k-kp  =y(Kig2|kik^!  dI^dKg(K3^K2|k^k3) 


2    2 


r    (Z^  2|k«_ky*  dK^(K^2  Ik^k^) 


(2.2.37) 


1    0 
On  evaluating  these  integrals  we  find  that 


^(k^kj;  kj^k^) 


,  kJkl^  =  6(ki"ki) 


k2+  k2-E 
1   2 


B   VE  -  kf  1 

"  VS  -  kf^  -  Bi   k2+  k^2  _j 


+  5(k2-kp|  . 


(2.2.38) 


Both  of  these  forms  for  the  Green's  function  will  be  needed  in  the  work 
which  follows. 
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It  is  important  to  note  that  equation  (2.2.38)  gives  the  relation 

CKk^k^;  k|kp  =  G(k^k^;  k^k^)  .  (2.2.39) 

This  symmetry  tetween  the  imprimed  and  primed  momenta  is  not  apparent  in 
equation  (2.2.37)'  However,  hoth  the  orthogonality  and  closure  relation- 
ships are  symmetric  in  the  interchange  of  the  eigenfxmctions  with  its 
complex  conjugate.   Therefore  the  expansion  may  be  made  in  terms  of  the 
complex  conjvigates  of  the  eigenfunction.   In  that  event  eqxaation  (2.2.37) 
would  display  the  property  given  hy  equation  (2.2.39) • 

The  solution  of  eq\mtion  (2.2.36)  hy  the  direct  method  leads  directly 
to  equation  (2.2.38),  so  that  there  is  no  need  to  repeat  the  details. 

Once  the  Green's  function  has  been  determined,  the  solution  to 
equation  (2.1.3)  nay  he  represented  by  the  following  intefral  eq-uation. 

(2.2.40) 
f(k^k2)  =  f^(k^k2)-(l/4n*)y'G(k^k2;k{k^)dk|dk«(kjk^  fV  |k3«'k2")dk3«'dk2"f(k^"k2") 

The  integral  on  the  right  hand  side  is  identical  with  that  obtained  for 
f,(k,kj,)  and  given  by  equation  (2.2.35).   In  (2.2.35)  the  term  in  the  square 
brackets  is  Just  the  Green's  function  and  the  definition  (2.2.29)  of  7(k-k  ) 
completes  the  identification.   Thus  the  integral  ea.uation  formulation  of 
the  problem  leads  automatically  to  the  exact  solution. 

2.3   The  5xact  Solution 

The  solution  to  eqxaation  (2.1.3)  given  by  (2.2.18)  represents  the  state 
of  the  system  in  momentum  ayiace.   To  obtain  the  solution  in  coordinate  space 
equation  (2.2.18)  must  be  inserted  into  (2.1.2)  and  the  resiiltr.ng  integral 
evaluated.   The  wave  function  in  coordinate  space  describes  the  state  of  the 
system  for  all  values  of  the  coordinates  i.  and  x^   .  This  information  is 
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much  more  than  we  require.  We  ere  interested  in  the  scc-ttered  wave  in 
coordinate  space,  I.e.,  for  large  values  of  x  end/or  x„.   Therefore  only 
that  part  of  the  wave  function  in  momentum  spece  which  contributes  to  the 
scattered  wave  in  coordinate  space  is  of  interest.   This  ;^i8rt  of  the  wave 
function  in  momentum  space  can  be  separated  from  that  portion  which  describes 
the  system  when  the  particles  are  close  to  each  other  by  applying  the 
Riemann-Lebesque  theorem  .   This  theorem  states  that  if  f (k)  is  bounded 

and  integrable  over  the  range  of  integration,  say,  a  ^  k  ^  b,  then 

b 
lim  J       f (k)  exp(ikx)  dk  =  0.  (2.3.1) 

i-v=oo   a 

An  immediate  consequence  of  this  theorem  is  that  non-zero  contributions 
from  integrals  of  the  form  given  in  (2.3.1)  for  large  values  of  x  come  from 
those  regions  of  the  variable  of  integration  where  f(k)  is  unbounded,  i.e., 
at  the  singularities  of  f(k).   Since  the  range  of  momenta  is  over  the  real 
axis,  the  singularities  that  are  important  are  those  that  lie  on  the  real 
axis.   Therefore  to  obtain  the  scattered  wave  in  coordinate  space  all  that 
need  be  known  is  the  value  of  the  wave  f\mction  at  its  singularities  on  the 
real  axis  of  momentum  space.   The  value  of  the  wave  function  at  its  singu- 
larities on  the  real  axis  of  momenttun  space  will  be  referred  to  as  the 
asymptotic  jjart  of  the  wave  function. 

The  asymptotic  part  of  the  exact  solution  f (k^k-)  is 

(2.3.2) 

^A    ,     N      ,A       rr   \         N  A     Ntt         C         ^       n  N  1 

f(kTk„)'%/6(k  -K   )      -5—2     -  — 2  —    -^ +  ~2~~2~T~~2. 

^  '^  ^     ^       k^+B'^         kv*^   B      k^+il  i(&fT)     k^+B     k^-K^ 

A  1  r    1^2!        _C_    ^jj     f  Bi ^     ll    " 

V     k^+k^E        Ikg^Ei     k^+lT  1   (|k2|-Bi)(|k2|+iT)  k^+ilj    ^ 

The  tilde  denotes  that  the  asymptotic  part  of  the  function  is  to  follow. 
Ihe  first  term  on  the  right  hand  side  of  the  equation  (2.3.2)  represents 


36.  H.  S.  Carsl?m,  Introduction  to  the  Theory  of  Fourier's  Series  and 
Integrals  (.Dover  Publications,  1930).  Third  Edition. 
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the  initial  state  of  the  system.  Die  second  and  third  terms  represent  the 
effect  of  the  interaction  on  the  initial  state  of  the  systeni.  To  exacnine 
the  significance  of  the  second  and  third  terms,  f (k-k^)  is  expressed  in  a 
mixed  representation  in  which  particle  1  is  descrihed  in  coordinate  space 
and  perticle  2  is  described  in  niomentvim  space.  This  representation  is 
obtained  by  evaluating  (2.1.2)  only  for  the  Integration  over  k, .  The  resiilt 


Is 


(2.3.3) 


•(x^.lTg) 


where 


e  =  - 


{'^  ?!?  e^CiK,^)  *  ^  exp(iK^x^)  *  b  exp-(i(^k^)l/2,J  ; 

(x^>0) 

t7o   "~5 — 9   expdK  x^  )  +  -§ — j  exp(-iE  x  )  +  d  exp/-i(E-k„)  '  x^  ( ; 

(x^  <  0) 

NlT  TTi       1      "" 


im^  B  K^      (2tt)^/^ 


^_  K +iT   i(B4-T) 


b  =  - 


TT  i 


i*!^  l^;;^  ^^::^ 


+  D-^ 


Bl 


(|k„|-Bi)   (/E-k^+iT)     (|k^l-Si)(|k^|+iT) 


k2+lT 


c  =  - 


kv 


A_  Ntt  TTi 
2 


B  1^   (2tt) 


l72 


K  -iT   i(B4-T) 

l_    0 


d  =  - 


A 


TTl 


^   (2t7) 


^2  (E-k2)l/2 


+  D^ 


Bl 


(Ikpl  -Bi)  (/&^-iT)    kikpl  -Ei)(|kpl+1T) 


kg+lT 
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The  one  dimensional  character  of  this  problem  is  indicated  by  equation 

(2.3-j)  where  J(x-,i:^)    is  given  for  forward  (x^  >0)  and  backward  (x,<  O) 

BCattering.   This  replaces  the  angular  dependence  which  arises  in  three 

dimensional  problems.   The  first  terms  in  (2. 3 '3)  represent  tine  incident 

wave  and  arise  from  the  factor  5(kT-E  )  in  moment-um  space.  The  second 

■'■  0 

terms  arise  from  the  factor  l/Ckr-K^)  and  express  outgoing  waves  with 
momentum  K  in  the  forward  and  backward  directions.  These  outgoing  waves 
are  the  result  of  elastic  collisions  since  the  energy  of  these  particles 
remains  the  same.   The  third  terms  which  arise  from  the  factor  l/(kr+  k~  -S) 
also  represent  outgoing  waves  but  with  momentum  |k^  |  =  l/E-k^  .  These 
waves  are  the  result  of  inelastic  collisions  since  the  energy  of  these 

particles  has  been  altered.   The  contributions  to  the  inelastically  scattered 

2 

waves  are  limited  to  values  such  that  E  -  k_^  0,  for  otherwise  these 

2 

scattered  waves  drop  off  exponentially.  It  is  evident  that  the  contributions 
to  the  scattered  waves  come  from  those  values  of  the  momenta  v;hich  satisfy 
the  principle  of  the  conservation  of  energy.  Of  course  the  physics  of  the 
situation  demands  this.   If  EO  there  is  no  inelastically  scattered  wave 
and  this  again  represents  the  physical  situation,  for  it  is  then  energetically 
impossible  for  both  particles  to  be  scattered. 

The  cross  sections  for  elastic  and  inelastic  scatter! rig  can  nov/  be  de- 
termined. Consider  the  inelastically  scattered  wave  in  eq.uation  (2.3.3)  • 
The  coefficient  of  the  exponential  term  is  the  probability  amplittide  and  its 
absolute  value  so-uared  gives  the  probability  of  finding  particle  1  at  x, 
with  momentum  |  k^^  |  =  Vs-k^  when  particle  2  has  momentum  in  the  range 

between  k  and  k  +  dk  ,  The  total  outgoing  current  density  is  then  given  by 
222  ^     -^ 
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The   total   incoming  c\irreiit  density  of  particles   1   is   determined  from   the 
prohabillty  amplltTide  of  the  incident  wave  and  is 

2 


/ 


K, 


m-i 


1        __r_ 

(2nlV2     ^^z 


^  =  Vn 


(2.3.6) 


The  cross  section  for  inelastic  scattering  is  ohtained  hy  dividing  (2.3'5) 
t>y  (2.3-6)  and  this  is 


^" 


\^\J^^ 


(2.3.7) 


In  a  similar  manner  the  cross  section  for  elastic  scattering  can  he  found 


and  it  is  given  by 

Q^  -  2TT|la|2  +  IcI^J  (2.3.8) 

Substitute  the  values  for  a,  b,  c,  and  d  given  by  (2.3.^)  into  (2.3.7)  and 
(2.3.8).  Then  the  cross  sections  become  eo_ual  to 


^ 


2 

8t« 


2^,r2  1     T        (B^ri)l72 


(2.3.9) 


K^+T 
0 


2T  +2T  E-B'^E 


T       i     B 

(B-KC)(Z^+T^)  l(B^+E)^/^    (2T^+B)T(T^+S)-'-/^^ 


CD*+  C*D| 


r 


BZ 


T(T^3)(T^+E)3/^   (T4B)^(T^+E)^/^   ( T+B)^(B^+E)^ 


for  S  >0,  and 


0  =  .^ 


r:^  t^*  ^  ^'^} 


Jl 


.   (2.3.10) 


K^+t2   (B-fT)(K^-tJT2)   L         .>'    (B-Hr)2 
From  (2.3.9)  it  is  easy  to  verify  that  Q^  ■  0  for  E  •  0  whatever  the  values 
C  and  D  might  have.  Physically  this  is  necessary  since  on  energetic  considera- 
tions alone  the  inelastic  scatteriiig  cross  section  must  vanish  for  all  values 
of  the  total  energy  which  is  less  than  or  equal  to  zero.  The  values  for  C 
and  D  shoxild  in  no  way  affect  this  basic  physical  req\ilrement. 
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It  is  clear  from  equations  (2.3-9)  and  (2.3.1-0)  that  the  cross  sections 
depend  on  the  constants  C  and  D.   The  approximation  procedures  considered 
below  are  used  to  deteraiine  the  values  of  these  constants. 

The  continuity  eq^uation  for  the  pro'ba'bility  density  leads  to  the  follow- 
ing relation  between  the  scattering  cross  sections 

q+n=--^In^l„  r  C   +   P   1  (2.3.11) 

®   ^      k'^     B      Z^  [T^+i-l!         i(B+T)J 

This  relation  is  an  instance  of  a  general  theorenr'  which  relates  the  total 
cross  section  for  scattering  to  the  imaginary  part  of  the  probability 
amplitude  of  the  wave  scattered  elastically  in  the  forward  direction.   This 
equation  has  been  independently  checked  by  inserting  the  values  for  C  and  D 
and  comparing  the  coefficients  of  like  powers  of  the  independent  parameter  A. 
The  right  hand  side  of  equation  (2. 3. 11)  may  be  used  to  calculate  the  total 
cross  section  but  we  are  interested  in  the  individual  cross  sections.  Never- 
theless eauation  (2. 3. 11)  serves  as  a  -useful  check  on  the  numerical  calcula- 
tions and  has  been  used  for  this  purpose. 

2.4  The  Bom  Approximation 

Consider  equation  (2.1-3)  written  in  the  following  way, 

(2.4-1) 

(k^*2.j;)f(i,_^j,^)_(B/^)y  f(k^k^)4k^  =  -(l/4n2)y(k^kj7|k'k')dk'dk^f(kjkp  . 

The  Born  approximation  to  the  exact  solution  is  obtained  by  solving  the 
equation  which  results  when  f(kLk  )  in  the  integral  on  the  right  hand  side 
of  the  above  equation  is  replaced  by  that  solution  of  the  unperturbed  equation 
which  renresents  the  initial  state  of  the  system-   The  vinperturbed  equation 


37.  B.  A-  Lippman  and  J.  Schwinger,  Phys.  Rev.  22..  ^9  (1950). 
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ia  the  left  hand  side  of  (2.^.1)  equated  to  zero  and  is  the  same  equation 
as  (2.2.1).   Its  solution  is 

^o^V2^  =  tCk^-K^)  N/(k2  +  B^).  (2./;. 2) 

The  Born  approxinuition,  f   (k^k^),  to  the  exact  solution  satisfies  the 
equation,  (2.^.3) 

(1.2  +  k^  -  B)f^^^k-k2)-(B/n)  /f^^\k,k2)dk2  =  -  A  |  ^+  -^— 


where 


and 


C^  =  ydk^dk2fQ(k^k2)  /  (kg-i^)  =  NTTi/B('iM-B)       (2.^.4) 
^B  =  /^Wo^'^1^2^  /  V^^^  =  NTT/B(K^-iT) 


Equation  (2.^.3)  is  solved  by  the  direct  method  and  yields, 
I    ».»^iK2''    ,  2  .  ^2      ,._2  ,  2.,  2 


k2  +  B 


/ITT    kg+kg-E 


(2./I.5) 


rs-lc^ 


+  Dt 


r 


Bi 


k^+iT   Vs_),2  _. 


Si 


Ve-I:^  -Bi )  (/l^k^+Ti  )     k2+ 


iT 


Its  asymptotic  part  is 
^(1) 


6(k,-K  )N      ,   ., 
T.  0       A   Nrr 


^^'^^'^  k^.B^   "  Utt^   B 


C       1^. 


k^+iT   i(BfT)   kg+B^  kJ-K^ 


(2./4.6) 


N 


?        2   2 
i4TT"  k^+kg-E 


Cp     Ikgl 


(k^+iT)  (ik2i-Bi) 


+  D 


Bi 


^  ^  (lk2l-Bi)(|k2i+iT)*k2+lT  I 
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If  thds  res^olt  is  compared  with  the  asymptotic  part  of  the  exact  solution, 
eouation  (2.3«2),  we  see  that  the  form  of  the  solutions  is  the  same.  The 
same  formulas  for  the  scattering  cross  sections  can  be  used  with  C  and  D 
replaced  by  C-g  and  D^  .   The  validity  and  range  of  applicability  of  the 
Born  approximation  depends  on  how  well  the  exact  constants  C  and  D  are 
approximated  by  C  and  D  .   This  can  be  determined  by  comparing  (2.4.4) 
with  (2.2.16).   It  is  clear  that  for  A^o  or  for  those  values  of  E  such 
that  the  prodxicts  of  A  and  J.(i  =  1,2,3,4)  are  very  much  less  than  1,  that 
the  results  obtained  from  the  Bom  approximation  will  be  useful  approxima- 
tions to  the  exact  results.   The  J.  are  functions  of  S  which  approach  zero 
for  sxiff iciently  large  E  so  that  the  products  AJ.  can  be  made  small  for  E 
sufficiently  large.  The  explicit  dependence  of  J^^  on  E  is  given  in  Appendix  I, 

These  results  are  well  known  and  represent  in  detail  the  motivation  for 
employing  the  Born  approximation.  Physically,  if  the  energy  of  the  incident 
particle  is  siofficiently  large  or  the  interaction  between  the  incident  and 
bound  particle  extremely  small  then  we  expect  the  initial  state  of  the  system 
to  ■undergo  a  negligible  amount  of  distortion.  Under  such  conditions  the 
initial  state  or  unperturbed  solution  should  be  a  good  approximation  to  the 
exact  solution.  There  is,  however,  one  point  that  needs  further  investiga- 
tion even  under  these  favorable  conditions.  This  is  the  problem  of  determin- 
ing what  is  meant  by  the  statement  that  the  energy  m\ist  be  sufficiently  large 
for  the  Bom  approximation  to  be  good.  Jost  and  Pais   in  an  analysis  of  the 
Born  approximation  for  a  Yukawa  potential  reached  the  conclusion  that  the 
value  of  the  energy  for  which  the  Born  approximation  gives  ^ood  results  is 
so  large  that  for  all  practical  purposes  the  approxiioation  is  useless  for  that 
potential.  A  similar  analysis  was  -undertaken  for  equation  (2.4.1)  and  the 

38.   loc.  cit. 
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results  lead  to  somewhat  the  same  conclusion.  Before  presenting  the  numerical 
results  of  this  work  an  analysis  of  the  iteration  procedure  which  starts  with 
the  unperturbed  solution  will  he  found  to  be  illuminating. 
The  iteration  procedure  is  equivalent  to  assuming  that 

f(k  k  )  =  r     a""   f^^'^C^k  )  (2.4.7) 

12      £—  12 

n-0 
This  power  series  in  A  will  he  henceforth  referred  to  as  the  Bom  expansion. 
Substituting  this  series  into  f2.4.l),  axid  equating  the  coefficients  of  like 
powers  of  A,  we  find  that  for  n  =  0, 

(k2+k|.3)  f^°^(k^kp  -  (BA)/f^°)(k^k2)dk2  =  0.  (2.i+.8) 


and  for  n  > 1 , 


(k2+k2-S)f^^^(k  k  )-(B/T,)/f(^)(k  kjdk,  «  -   1  /  ?ii=i-  +  ^JarL.    \ 
T.2         12       ^  1  2       2  i^n^      ^     j^_^^ij   ^^^^^   J 


(2.4. Q) 


in  which 

Vi  =  /^i'^2  f^^-'^C^^a^^V^^^ 

and  (2.4.10) 

Vl  =  /^^2  f^'^-'^(l^ik2)/(k^-iT)  . 

Eouation  (2.4,8)  represents  the  unpertxirbed  system  and  its  solution  is  the 
initial  state  of  the  system.   3q\aations  (2.4.9)  represent  the  additional 
contributions  to  the  initial  state  that  arise  from  the  interaction  between 
the  particles.   Hence  the  condition  which  must  be  placed  on  f^  '(k..k  ),   for 
n  >  1,  is  toat  they  represent  outgoing  particles. 
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Equations  (2.4.9)  can  be  solved  by  the  direct  metkod.  We  find,  for 
n  >  1,  that 


1   2  i^j^ 


Cn-l       V^^l 


(2.it.ll) 


k^+il     Y^Ik^    -Bi 

-  D^_,/ Bi ,  _ijn    — L_ 

((-j/E-k^  -Bi)  (Vs-k2  +  Ti)        k  +iTJ  k^*  k^- 

which  can  be  written  in  matrix  notation  as 

Q 

f('^^(k^k2)  =   (h.g)    ^    ^-^"j    ;        n   >1. 


(2.J^.12) 


where 


h  3  h(k  k  )  =  - 
1  Z 


Vffi-kf 


tl- 


4 1?    VB-k2  -Bi       k^+iT       k2+k2-B 


(2.if.l3) 


^  ^  ^(^^2^   =  -     6 


r 


Bi 


4"     1  (Vu-k?  -Bi)(V2-k?  +  iT)       k  +iT  1    k^*  k^-E 
Inserting  (2.4.12)    into   (2.4.10),  we  obtain 

C     »       c<C     .      +      P  D„   , 
n  n-1  n-1 

n  (>    n-1  n-1 


(2.4.14) 


where 


c<  =  y  clk^clk2h(k^k2)y'(k2-iT) ;     P  =  /  dk^dk2g(k^k2)/(k2-iT)  (2.4.15) 

^  -  /  dk^dk^h(k^k^)/(k^-iT) ;    6  =  /  dk^dk2g(k^k^)/(k^-iT)    . 


Equation  (2.4.14)  becomes  in  matrix  notation 


2r   6 


n-1 


i  n  >1, 


D 


n-1 


(2UI.16) 
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which  gives 

C_ 


y      6 


Hence  equation  (2.4.12)   "becomes 

.  n-l 

f^^^Ck^y  =  (h 


■"(;:)  ( 


;   n  >  1. 


Substitute  (2.4.18)  into  (2.4.?)  and  sum.  The  result  is 
f(k3^k2)  =  fQ(k3_k2)+  Ain£-[^(l-i£)Co  +  ApD^J  h(  k^kg) 
+  I^AKCp  +  (l-iPt)DoJ  gd^k^)  "] 


with  the  determinant 


m 


-1 


1  -  A<=<     -AP 
-  A  8"   1  -  A6 


and 


C^  =/dk^dk2fo(k^k2)/(k2-iT)  «  Nni/B(B+T)  =  Cj  . 
^0  =/<^i^2^o^^l^2^/^V^^^  =  NTT/B(Z^-iT)  =  D^   . 


(2.4.17) 


(2.4.18) 


(2.4.19) 


(2.4.20) 


From  equation  (2.4.8)  it  is  evident  that  f  ^^(kj^kg)  is  eo_ual  to  f^(k^k  ). 
The  evaluation  of  the  integrals  in  (2.4.15)  gives 

<^=  4=  iJj^/4TT   ;       ^C  =  J/^/^'rfr  (2.4.21) 

p  =  -  (BJ^  +  iJ3)  /  4tt  , 

where  J.  (i  =  1,2,3.4)  are  the  ssjne  quantities  which  occ\ir  in  (2.2.16). 
Combining  equations  (2.4.19),  (2.4.20)  and  (2.4.21)  we  obtain  the  exact 
solution  given  by  equation  (2.2.18).   From  (2.4.21)  it  follows  that  m  ■  /^  , 
the  latter  quantity  appearing  in  equation  (2.2.1?). 
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The  series  solution  obtained  by  the  iteration  procedure  converges  to 
the  exact  solution.  The  series  itself  represents  the  exact  solution  only 
for  those  values  of  the  parameters  which  lie  within  the  radius  of  convergence 
of  the  series.   In  practice,  where  only  a  few  terms  in  the  Born  exT>ansion  can 
he  eval-uated,  the  domain  of  applicability  or  validity  of  the  resulting  approxi- 
mation is  restricted  by  the  radius  of  convergence  of  this  series  representa- 
tion.  The  radius  of  convergence  of  the  Born  expansion, /O ,  taken  for  conveni- 
ence to  be  expressed  as  |A/4Tr|  ,  is  given  by  the  smaller  in  absolute  value  of 
the  two  roots  for  kjkv   of 

^"'  =  m-l  =  (1  -  ^  J,i)2.  (_A)'  h.   (BJ^.,  iJ3)  -  0.  (2.i..22) 

The  graph,  Figure  1,  is  a  plot  of  the  radius  of  convergence  P  as  a 
function  of  the  energy  in  units  of  B  i.e.,  Z  =  2/B''.  To  obtain  the  curve, 
T  Was  set  eatial  to  2B.   The  interpretation  of  the  curve  is  as  follows:   For 
a  given  Z  ,  the  series  converges  for  all  values  of  |A/^'rf|<  P.  where  P-,      is 
the  value  associated  with  Z,  ,  by  the  curve.  Of  if  P-i  is  given  and  Z,  is 
associated  with  tnis  value,  then  for  all  Z  >  Z,  ,  the  series  solution  converges. 
It  follows  that  for  a  given  a/4tt  the  best  that  can  be  expected  from  the  series 
representation  is  the  value  of  the  exact  solution  for  all  Z  >  Z^  where  Z,  ,  is 
the  energy  associated  with  the  absolute  value  of  the  given  a/4tt  .  For 
Z  <  Z, ,  the  series  representation  diverges.   For  large  values  of  the  total 
energy  E  the  radius  of  convergence  approximates 

/O  A'T(E)-'-/^  /tt   =  BT(Z)^/2/tt  (2.4.23) 

This  means  that  in  this  region  for  the  total  energy  a  small  change  in  A 
reojoires  a  very  large  change  in  the  value  for  E  for  which  the  series  would 
still  converge. 
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°  a 
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The  Born  approximation  is  obtained  from  the  series  expansion  (2.4.7) 
if  f^^'(k^k  )  is  set  equal  to  zero  for  all  n.  ^  Z   .     The  series  is  then 
represented  by  its  first  two  terms.   This  approximation  should  be  good  for 
"energies  stifficiently  large."  Using  the  expressions  obtained  in  Appendix  I 
for  J.  ( i  =  1,2,3,4)   for  large  values  of  S,  in  the  exact  formulas  for  C  and 
D,  we  find  that  C  is  approximated  very  well  by  C^  but  that  D  does  not 
approach  IL  unless 

a/^tt  «    T(B+T)/tt  (2.4.24) 

The  constant  B  Is  related  to  the  bound  state  energy  whereas  the  constant  T 
is  related  to  the  range  of  the  interaction  in  its  space  coordinate  representa- 
tion. Note  that  both  relations  (2.4.23)  and  (2.4.24)  depend  on  the  range  of 
the  interaction.  Equation  (2.4.23)  assures  the  existence  of  a  radius  of  con- 
vergence for  the  Born  expansion  no  matter  how  small  T  is  taken,  i.e.,  for 
long  range  interactions,  if  the  total  energy  E  is  taken  sufficiently  large, 
say,   E /v  kT"'-  with  k  constant.   However  (2.4.24)  asserts  that  no  matter 
how  large  E  is  taken  if  this  relation  is  not  satisfied  then  the  Born  approxi- 
mation has  no  range  of  validity.   Tor  the  Bom  approximation  to  yield  useful 
results  two  conditions  must  be  met:   first,  the  energy  must  be  large  and 
secondly  the  interaction  must  be  weak.  By  weak  we  mean  the  satisfaction  of 
equation  (2.4.24) . 

The  scattering  cross  sections  (2.3.9)  and  (2.3.IO)  for  large  values 
of  E  become 

\-     ^     IcIYsttTE  (2.4.2^) 

euid 


Ci  ■  {i^'ilW^   E) 
e      ' 


+ 


E      (B+T)' 


(2.4.26) 


-  3^  - 

Since  the  exact  C  arproaches  the  Bom  velue  C-g,  the  inelastic  scattering 
cross  sections  for  these  solutions  will  be  in  good  agreement  for  large 
values  of  the  energy.   However,  unless  (2.^.24)  is  satisfied,  D  does  not 
approach  Hg.  The  exact  and  Bom  calculation  for  the  elastic  scattering 
cross  section  can  not  possibly  approach  each  other  no  matter  how  large  E 
is  taken. 

Nezt  consider  the  second  Bom  approximr.tlon  which  we  define  as  the 
first  three  terms  in  the  Bom  expansion.   For  large  values  of  the  total 
energy  E  the  constants  C  and  D  are  approximated  by  the  expressions  C^o  ^nd 
D   which  are  given  by 

C^     -  C^  (2.4.27) 

^^  '  '  (2.4.28) 


D    s  D   /  1  +  ^ ~ 1 

2B      B  |_      5tt  T(B+T)  J 


These  values  are  precisely  the  values  for  C  ajid  D  for  large  values  of  E. 
Thus  the  scattering  cross  wections  calctolated  for  the  second  Born  apnroxi- 
mation  will  agree  with  the  exact  values  if  S  is  sufficiently  large.  Further- 
more equations  (2.4.27)  and  (2. 4. 28)  are  not  restricted  by  any  condition  on 
the  magnitude  of  A.   This  is  a  distinct  improvement  over  the  Born  approxima- 
tion.  As  a  matter  of  fact  the  ratio  of  the  elastic  scattering  cross  sections 
calculated  from  the  second  Lorn  approximation  to  that  calculated  from  the 
Born  approximation  is 

l_/i  +  n  4.  A     TT    .  \  (2.4.29) 

for  large  values  of  E.   This  ratio  is  a  constant  and  it  indicates  that  in 
general  the  second  Bom  approximation  is  a  substantial  correction  to  the  Born 
approximation. 

The  significance  of  the  second  Eom  approximation  resides  in  the  fact 
that  it  takes  into  account  the  second  and  third  terms  in  the  Born  expansion. 
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These  two  terms  are  the  leading  terms  in  a  power  series  expansion  for  the 
outgoing  scattered  wave.   Their  ratio  can  te   used  as  a  measure  of  the 
rapidity  of  the  convergence  of  the  series.  This  ratio  also  expresses  the 
pertinent  physical  aspect  of  the  problem.  The  first  term  in  the  Born 
expansion  represents  the  initial  state  of  the  system  and  is  independent 
of  the  interaction  between  the  particles.  It  is  the  second  term  in  the 
expansion  which  first  exhibits  the  effects  of  the  interaction  on  the  initial 
state.  Hence  to  get  some  idea  of  the  total  effect  of  this  interaction,  the 
third  and  second  terms  of  the  expansion  must  be  compared.  A  similar  point 
will  arise  when  Bom  type  trial  f-unctions  are  used  in  the  Schwinger  varia- 
tional principle. 

The  problem  treated  by  Jost  and  Pais  led  to  results  similar  in  content 
to  equations  (2.^.2^  and  (2. 4'. 29).  Furthermore  Dalitz  evaluated  terms  in 
the  Born  expansion  of  the  problem  he  considered  up  to  and  including  the 
third  order  terms.   The  terms  to  the  second  order  are  sufficient  to  indicate 
the  convergence  of  the  expansion  for  the  wave  function.   Hence  the  above 
results  are  not  peculiar  to  the  specific  problem  studied. 

Table  I  contains  the  exact  values  for  the  elastic,  0  ,  and  inelastic, 

(^  ,    scattering  cross  sections  for  various  values  of  the  total  energy  E  with 

the  parameters  A/^tt  set  equal  to  1  and  T  eq-ual  to  2B.  For  these  same  values 

of  the  parajmeter  the  radius  of  convergence,  P,   has  been  tabulated  as  well  as 

the  elastic,  Q  _,  and  inelastic,  C^,-,  scattering  cross  sections  calculated 
eB  iB 

from  the  Born  approximation.  The  Bom  results  are  \msatisfactory  for  P<  1, 

The  general  discussion  of  the  Born  expansion  would  lead  tis  to  expect  this 

resTilt.  For  P^l  Q,^^-  improves  steadily  while  Q,   fails  to  approach  the  exact 

0  even  for  large  values  of  the  energy.  These  results  are  confirmations  of 
e 

the  discussion  given  above. 

The  preceding  work  indicates  that  in  general  to  obtain  useful  results 
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higher  order  terms   in  the  Born  expansion  anist  be  calculated.      In  practice 
this   is  very  difficult   to  do,   so   that   other  uethods  mist  be  used.      Among 
these  are  the  variational  procedures  ana  in  particular  the  Hulth^n-Kohn 
and  Schwinger  vsiriatlonal  principles  which  will  now  be  considered. 

2.5     H\il thIn-Kohn  Yaristional  Principle 

Equation   (2.1.3)  with  the   interaction  given  by    (2.1.5)   is  reconsidered. 
Tor  convenience   these  equations  are  again  written  down. 

(2.5.1) 

(k^+k|-E )  f  (k^kg )-  (  B/n  )yf  (k^kg  )dk2+  ( 1/^-n^  )/(^^2 1  ^  1  M^2  ^^^^2^'  ^H^2  ^  =   ° 

and 

(k,kJTlk.«k')  =       A  J ^ + ^^ I  (2.5.2) 

^   '^  ^   '^  1  (k^+iT)(k|-iT)  (k2+iT)(k^-i?)J 

In  terms  of  the  Green's  function,  equation  (2.5.1)  can  be  represented 

by  tlie  following  integral  equation, 

(2.5.3) 
f  (k^kg)  =  f^(k^k2)-(l/^^/G(k^k2;k|k|)dkjJdk^(k^k^|vlk^«'k2«0<ik^''dk2«'f  (k^''k2«' ), 

where  f  (k^k^,)  is  a  solution  of  the  unperturbed  Schroedinger  equation  and 
represents  the  initial  st-ite  of  the  system.   The  Green's  function 
G(l'  k2;k«k«)  is  given  by  equation  (2.2.37)  or  (2.2.38).  With  this  fvmction 
(2.5.1)  can  be  written 

th  K  =  B  +  E  and 


wi...  „^ 
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Table  I 


Exact 


Radius  of 
Gon7erip;ence 


Born 


% 


%1 


MB 


-.9 

1.981103 

.153299 

84.3468 

-.5 

1.855550 

.363870 

14.72996 

-.1 

1.6556U7 

.518062 

7.22660 

-.05 

1.622907 

.536320 

6.74632 

-.01 

1.592725 

.550914 

6.39884 

-.005 

1.588288 

.552750 

6.35743 

-.001 

1.58^3^9 

.554227 

6.32468 

+  .001 

1.581312 

.000^+76 

.554833 

6.30840 

.0016854 

+  .005 

1.57^897 

.002350 

.555733 

6.27604 

.00842096 

+  .01 

1.567^23 

.OOi+652 

.556832 

6.23596 

.0167470 

+  .05 

1.515046 

.022>i2 

.565409 

5.93000 

.0789516 

.1 

1.^58^19 

.042524 

.57573 

5.58064 

.147110 

.5 

1.150270 

.152834 

.65061 

3.67353 

.451456 

.9 

.962272 

.210601 

.71706 

2.62637 

.549148 

5. 

.  33852i^ 

.236493 

1.23844 

.409404 

.356102 

10. 

.15^223 

.165691 

1.73059 

.l?'''c96 

.205592 

100. 

.002788 

.021490 

6.10792 

.001718 

.020981 

lOOC. 

.000029 

.002190 

20.03332 

. 000018 

.002192 

•J 

1     =  E/B   ,   E  is   the   totel  energy. 
Q^  =  exact  elastic  scattering-  cross  section. 
0^   =  exTCt  ineltstic  scattering  cross  section. 
0  ^  =  Bom  elastic  scs.tterin^-  cross  section. 

^     =  Bom  inelEstic   scattering  cross  section. 

2 

p  =  radiiis  of  convergence  of  the  Born  expansion  in  units  of  B  . 
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f(K^)  =  -(l/4Ti2)y(K^2|k|kp*dk|dk^(k|k^|Vlk^>'k2")dk^''(ik2''  fCk^'-k^"),   (2.5.5) 
fdC^K^)  =  -  (l/^TT^)  y(K^K2|kjk^)*dkj<ik^(k|k^|Vlk^«n-^«')dk^"dk2"  f(k^"k«p. 

Associated  with  the  Integral  equation  (2.5.3)  i^  an  adjoint  integral  equa- 
tion which  represents  solutions  of  equation  (2.5.I)  expressed  in  terms  of 
incident  waves  plus  incoming  or  concentrating  waves.  The  Green's  function 
for  these  solutions  will  be  indicated  hy  a  subscript  c.   The  relationship 
between  the  G-reen's  functions  is 

G  (kj^k2;  kjkp  =  G*(k{k^;k^k2)  (2.5.6) 

The  integral  equation  representation  of  these  solutions  is 

(2.5.7) 
f^^(k^k2)=f^^(k^k2)-(l/^VG^(k^k2:k{kpdkjdk^(k{k^|V|k^"k2")dkt^dk2"fi(k^»'k2»') 

where  f  ,(k,k.)  is  a  solution  of  the  unperturbed  Schroedinger  eqxxation. 
ol   1  2 

The  specific  solution  desired  will  be  indicated  later.  This  equation  becomes 

(2.5.8) 

f  i(  k^k^)  =f  ^^(  k^k^) -( l/^^)/f  ^(  k^M  k^i . )  dk^ii dk^' '  ( k^' •  k^ '  I V 1 44) * dkjdkp* ( k^k^ :  k^k2) 

when  use  is  made  of  the  hermitian  property  of  the  interaction  and  the 

relation  (2.5.6) . 

\  39 

The  H\ilthen-Kohn  variational  principle   is  based  upon  the  considera- 
tion of  the  following  functional, 


h^^ 


fv^.v!-   11m   y*v*  (k^k^:  E')/h-  E|v(k^k2;  E)  dk^dk2 .   (2.5.9) 


39.   H.  S.  Moses,  Phys.  Pev.  Q2.  81?  ( 1^53) • 


-  39  - 


in  which 


+   (lATf)y(k^k2|v|k[kpdkjdk^  v(k^k^;E).  (2.5.10) 

The  functions  v(  k  k   ;E)   and  v  (k  k   ;E)   are  arhitrary.      The   limiting 

process   indicated  in  (2.5.9)    is   to  te   executed  after   the   integrations 

have  teen  performed.      This   is  necessary  for   the  T(k  k  ;E)   and  v(k  k   :E') 

12  12 

are  singular  and  to  treat  the  singularities  in  a  consistent  manner  (see 
footnote  39)  the  above  limiting  procedure  must  "be  used. 

Equation  (2.5.9)  ^^iH  ^e  varied  about  those  solutions  of 

£h  -  eJ  fds^kg*  ^)  *  0  •  (2.5.11) 

|h  -  E»}  f^(k^k2;E«)  =  0  .  (2.5.12) 


and 


which  are  represented  by  the  integral  equations  (2.5.3)  ajid.  (2.5.8) 
respectively.   The  variation  of  (2.5-9)  about  these  solutions  yields 

6Ii  =  lim  y6f(k^k2;E)  |E'-EJf^*(k^k2;E')dk^dk2  .  (2.5. 13) 

The  variations  about  the  solutions  of  (2. 5. 11)  and  (2.5.12)  are  now 
restricted  to  have  the  following  form  (2.5.1^) 

6f(kk:E)=  /-(Ki  2|k^k^)dZi6f(Ki)  ^  f      (KiKg|k^kg)dKidXg6f(g^K3) 

i  o  12 

Variations  of  this  type  rec^uire  that  the  inhomogeneous  term  of  the  correct 
solution  is  to  be  unaltered.  These  variations  have  the  further  property 
of  admitting  trial  functions  which  differ  from  the  exact  solution  only  in 
the  amplitudes  of  the  scattered  waves .  Substitute  (2.5.1^)  into  (2. 5. 13) 


-  ^  - 


choose  f^^  so  that  f^;^  =  (Z|  2|k^k2)  with  Kj^  .  b^S«  and  f^^   =  (K^K^lk^]^^) 

with  K«^K«2  =  E'.   Then  we  find  that 
1   2 

6(1^  -  f(Ko))  =  0  ;        k2  =  B^E  (2-5-15) 

and 

id^  -  f(K3_K2)  )  =  0  :      K^  +K^  =  E  (2.5.16) 

These  equations  assert  that  trial  functions,  which  we  denote  by  a  sub- 
script t,  which  differ  from  the  ezact  functions  in  the  first  order  determine 
f(Kjj)  and  f(K,K  )  ,  (K?+ K^  =  E)  correctly  to  the  first  order.  Thus,  for 
trial  functions,  equations  (2-5.15)  and  (2. 5. 16)  take  the  form, 

I    -  f  (K  )  =  -  f(K^);  K^  =  B^E 
It     to       ^0"   0 

and  (2.5.17) 

It  should  be  noticed  that  these  equations  determine  approximations  to 
f(K.)  and  f(K,K-)  only  in  a  restricted  domain  of  their  arguments.  But  the 
values  so  determined  are  J\ist  those  needed  to  obtain  the  asymptotic  part  of 
the  wave  function.   Equations  (2.5.1?)  determine  the  functions  at  the  sing- 
•ularlties  which  occ\ir  in  the  integrands  of  equation  (2.5.^). 

One  disadvantage  of  this  variational  principle  is  apparent.   Since  f(K,) 
and  f(K-,ii_)  are  known  approximately  only  for  a  restricted  domain  of  their 
arguments,  an  iterative  procedure  to  improve  the  trial  functions  by  sub- 
stituting the  variational  results  into  ea^uation  (2.5.^)  would  not  be  satis- 
factory.  The  amplitudes  of  the  scattered  waves  depend  on  the  interaction  of 
the  particles  when  they  are  close  to  each  other  and  this  means  that  f(K^) 
and  f(iLK-)  should  be  approximately  known  thro\;ighout  the  range  of  values  of 
their  respective  arguments.   This  knowledge  is  not  forthcoming  from  the 
variational  principle. 
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The  mathematicel   simplicity  of   the  Hulth"^n-Kohn  veriational  principle 
makes   it  feasible    to   employ  a  very  large   class   of   trial  fijnctions.      However, 
we  shall  be   interested  in  comparing   the   results  with  those   obtained  from 
the  Born  approximation  and  from  the   Schwinger  variational  principle.      The 
trial  functions   are   therefore  restricted   to  be   of   the  Bom  type,    i.e.,   the 
solutions   of   the  unperturbed  equation.      For  such   trial  functions   the 
Kulth^n-Kohn  variational  principle  cpn  readily  be  shown  to  give  the  same 
results   as   the  Born  approximation.      This  follows   immediately  from  equations 
(2.5.17)   when  compared  with    (2. 5. if)  and    (2,5.5).      For  our  purposes   the 
Hulth^-Kohn  procedure   is    the  equivalent   of   the  Born  approximation.     However, 
in  the  following  discussion  of   the   Schwinger  variational  principle  we   shall 
note  certain  similarities   between  these  principles  which  will  be  further 
explored  in  the  succeeding  section, 

2.6     The  Schwinger  Variational  Principle 

The  integral  equations    (2.5.3)  aJid.   (2.5.8)  are  rewritten  in  the 

following  way, 

(2.6.1) 

f  ^ (k^k2)  =  f  (ki^s^"^ (l/4TT^)/'G(k^k2;k|"k^"  )d«^«'dkn«  (k«^»k|«  1 Y 1  kj^i'k^*' ) 

•   dk^««dk2"f(lc^"k:2")» 
and  (2.6.2) 


^oi  ^^^2^  =  ^i  ^^1^2^"^  (l/^TT^)_/f  i  (kj^"k2«'  )cLki"c-k2"  (k^^kg"  JV  |k«««k«« ) 


>» 


•  dK««dk^"G*(k»'k«2";k^k2). 

I^ultiply  these  equations  with   (kik'lvlk^k  )dk^dk     and  integrate  over  the 

indicated  variables.      The   resulting  equations   csn  then  be  written  in  the 

forms 

a(k^k2)  =  K  f  (k^kp  (2.6.3) 

and 

^i^W^  "  ^'  ^i^^i^P  (2.6.i+) 
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wnere 


a^Ck^k^)  =/(k^k2|V|k«k«)dk«dk^f^^(k|Kp.  (2.6.6) 

and  K,  K'  are  integral  operators  with  kernels 

(2.6.7) 

Z(k^k2:k^kp  =y*[6(k^-kJ')6(k2-k'2«)+(l/^^)(k^k2|V|ki|''kn')dkn'dk'|«:}(k!^»ky';k|«k«)] 

•  dk^dkJ(kpc^|V|k^kp, 

aad  (2.6.8) 

K<(k^k2;k|kp«y  (kjk^lVlk^^)*dk^dkn   6(k^-k^)5(k2-k«) 

+  (i/4TT^G*(k"k";k"«k«^dk'"dk"«(k">k"'  'Vlk^k  )*l 
respectively.  The  meaning  of  eouation  (2.6.3)  is  that 

a(k  k  )  =  /K(k  k  ;  k«k')dk«dk«f(k'k') 
12   **  12   12   12   12 

with  a  similar  interpretation  for  equation  (2.6.^).  From  (2.6.7)  ajicl 
(2.6.8)  it  follows  that 

Z«(kj^k2;  k^kp  =  K*(kjk^;k^k2)  (2.6.9) 

which  asserts  that  K'  is  the  hermitian  adjoint  of  K.  From  (2.6.3).  (2.6.^) 
and  (2.6.9)  we  obtain 

(f^.a)  =  (f^.Kf)  =  (K'fi.f)  =  (a^.f).  (2.6.10) 

in  which  the  expression  (x,y)  represents  the  hermitian  inner  prod\ict  of 

the  q\iantities  embraced.  The  dependence  of  the  quantities  in  (2.6.10) 

on  the  variables  k  ,  k  has  been  omitted  for  the  sake  of  brevity  and 
1   2 

simplicity. 

The  Schwinger  variational  principle  (see  footnotes  7,8)  states  that 
the  expression 

X  |v^.  vl   -  (v^.  Kv)/(vj^.a)(a^,v)  (2.6.11) 


-  k3- 

is  stFtionary  for  independent  variations  of  the  arbitrary  f\inctions 

r  and  v.  about  the  solutions  f  and  f .  of  the  integral  equations  (2.6.1) 

and  (2.6.2).  For  v  =  f  end  v  =  f 

X  =  \  ^f^.f}  =  l/Ca^.f)  =  l/(f^.a).        (2.6.12) 
Note  that  for  t  =  f  and  v,  arbitrary 

X  -f^i.f}   =  l/(ai.f).  (2.6.13) 

The  stationary  property  of  (2.6.11)  serves  as  a  variational 
calculation  for  l/(a.,f)  or  l/(f. ,a).  The  term 

(a^.f)  =/f*i(lc^k2^  dk^dk^dc^k^lvlklkpdk^dk^fdclkp,        (2.6. l4) 

corresponds  to  f  (IL)  or  f  (tZ-)  of  (2.5.5)  provided  that  f  AK^2^   can 
be  replaced  by  (K  2|k  k-)  or  (K-Kplk^k-)  respectively.  However,  f   (k  k-) 
must  be  a  solution  of  the  unperturbed  wave  equation  whereas  (K.2|.k^k_) 
and  (K-Kplk^kp)  are  not  solutions  unless  KT  =  i  =  B  +E  for  t/ie  former 
and  B^  +  Kp  =  E  for  the  latter.  If  these  fxinctions  are  used,  with  the 
values  of  their  arguments  so  restricted,  in  place  of  f  .  C-c.k-)  in  (2.6.1^), 
then  -(a.,f)/^Tr   takes  on  the  values  f(K  )  end  f(K^K2),  which  are  the 
values  of  the  ezact  functions  at  the  sing-olerities  of  the  integrands  of 
(2,5.^).  Exactly  the  same  result  was  obtained  in  the  discussion  of  the 
Hulthen-Kohn  varintional  principle.  Hence  the  Schwinger  variational 
principle  yields  approximations  to  the  ainplitudes  of  the  scattered  waves 
only. 

The  exact  expressions  (2.6.1^)  are  equal  to 

(a. ,f )  =  aF  CP^^  +  DP^^  ]  ,   i  =  1,  2  (2.6.15) 


where 

(2.6.16) 

The  left  hpnd  side  of  equation  (2.6.15)  is  approximated  by  the  Schwinger 
variational  principle.  This  equation  then  yields  the  approximate  valuea 
for  the  constants  C  pnd  D.   In  this  way  we  obtein  the  exact  solution 
except  for  the  values  of  the  constants  C  and  D.   The  approximate  values 
for  C  and  D  may  be  inserted  in  the  formulas  for  the  scattering  cross  section 
to  obtain  the  variational  results. 

Before  discussing  the  results  obtained  from  the  Schwinger  variational 
principle  when  Born  trial  functions  are  used,  it  is  worthwhile  to  consider 
the  expression  for  \  j v.  .vh  in  detail. 

Let  P  and  Q,  be  two  integral  operators  with  kernels  l/(k^-iT)  and 
l/(kp-iT)  respectively.  Then  for  any  function  f(k^k2), 

P(f)  =  ydk^dk2f(k^k2)/(k^-iT)  (2.6.17) 


and 


Thsn 


Q(f)  5  ydk^dkgf (k^k2)/(k2-iT)  .  (2.6.18) 


(v^.Zv)  =  A  [p*(v^)Q(v)  +  0*(v.)P(v)  ]  (2.6.19) 

+  (A2/i+TT)[(-iJ^)(P*(v.)q(v)+n*(v^)p(v)  ) 

+  (BJ2^iIj)P*(v.)P(v)  -  (Ji^/T)/(v^)Q(v)  ], 
where  J.  (i  =  1,2, 3, U)  are  defined  in  equation  (2.2.16)* 

(vj.a)  =  A  [p*(v^)q(a)+Q*(v^)P(a)],  (2.6.20) 

(a^,v)  =  A  [^P*(a^)q(v)+q*(a^)P(v)l  .  (2.6.21) 


and 
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If  we  asstuae  that  the  trial  f\aiction  v  is  symmetric  or  antisymmetric, 

then 

P(t)  =  ±     q(v).  (2.6.22) 

The  plus  sign  is  to  be  associated  with  a  symmetric  v  end  the  minus  sign 

with  an  antisymmetric  v.  This  convention  will  he  adhered  to  in  what 

follows.  From  equations  (2.6.19),  (2.6.20),  (2.6.21) , (2.6.22)  and  (2.6.11), 

we  find  that, 

f.  -.     2f(A/iHT)  |l2iJ  +  (BJ^+iJ^)  ±  (-JjT)  J 

X  /  V.  ,v  \  =      z ~~--r ^ — i — ^^^^  (2.6.23) 

L  ^  J      A  l_q(a)  ±  P(a)J  [p  (a^)  ±  Q  (a^)  J 

Here  the  Schwinger  variational  principle  is  independent  of  trial  functions 
which  are  symmetric  or  antisymmetric. 

The  physical  problem  that  is  being  considered  here  cannot  have  either 
a  symmetric  or  antisymmetric  solution  because  as  formulated  particle  1  does 
not  interact  with  the  center  of  force  btfc particle  2  does.  For  this  reason 
the  results  obtained  from  (2.6.23)  wo\ild  be  expected  to  yield  very  poor 
approximations  to  the  scattering  cross  sections.  Certainly  these  trial 
functions  are  not  "close  "  to  the  exact  wave  function  for  which  case  we 
would  expect  to  get  good  results  from  the  vaxistional  calculations.  How- 
ever, the  numerical  calculations  given  in  Table  III  obtained  by  using 
symmetric  trial  functions  give  surprisingly  good  results.  This  indicates 
that  some  criterion  other  than  "closeness  <•  is  necessary  to  determine 
the  appropriateness  of  a  given  trial  function. 

Next  consider  the  principle  when  Born  type  trial  frinctions  are  used. 
Then 


which  becomes,  when  (2.6.7)  is  used. 


{^oi'^o}  =  ^^oi'  ^o>/  ^^oi'-)"' 


=  [(^oi'  ^^0^  *  ^  ^oi'^^^0^3  /  ^^oi-^^0^^^.,,   ^2-6-2^^ 
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where 


and 


VG7f  =  (l/ATT^)y'(k^k2lvlk»'k«0ik«'dk«2»'(J(k«^«'k"|;k«^«k«2»)dk^   • 

(k^"l?  1  k^kpdk|dk^f  ^  (k^kp . 
From  (2.6.24), 

where  equation  (2.6.6)  has  "been  Tised  to  obtain  the  last  equation.  Mow 

f  -  GYf  is  the  Born  approximation  to  the  exact  wave  function  f .  Eauntion 

0        0 

(2.6.25)  shows  (see  footnotes  7,  15)  that  the  Schwinger  variational  pro- 
cedure contains  the  second  Born  approximation  to  the  scattering  amplitude. 
Since  the  expajislon  which  led  to  (2.6.25)  implies  the  existence  of  s   radius 
of  convergence,  it  would  appear  that  the  Schwinger  variational  procedure 
which  is  not  so  rf>stricted  contains  more  useftil  information  then  the  Bom 
approximation.  However  (2.6.25)  is  valid  only  if 

I  (^oi'^^^^oV  ^^oi'  ^^o^  I  <1  (^•^•26) 

This  ratio  is  nothing  else  but  the  r?tio  of  the  asymptotic  values  of  the 
third  to  second  terms  in  the  Born  expnnsion  and  we  have  been  in  section 
(2.4)  that  this  ratio  need  not  satisfy  (2.6.26).   Tnerefore  unless  the 
same  rather  severe  restrictions  of  that  section  rre  satisfied,  the  comparison 
between  the  Schwinger  variational  results  and  the  Born  results  cannot  be 
made.  The  expected  second  Born  approximation  which  we  know  to  be  usef\il 
is  not  obtainable  from  the  Schwinf:er  variational  principle. 

In  Appenaix  II,  the  expressions  for  C  and  D  obteined  from  the  Schwinger 
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varietional  principle  are  given.  For  A  small  end  E  Icrge  the  Schwinger 
results  are  the  srme  as  the  Born  results  which  are  adeq^uate  under  the 
stated  conditions.  But  for  A  of  reasonable  magnitude  we  found,  in  section 
(2. A),  that  this  approximation  had  no  range  of  applicahility.   In  Appendix  II 
we  display  several  Born  type  trial  functions  which,  for  large  vslues  of 
the  energy,  make  the  varigtional  result  reduce  to  the  Born  result.  By 
stating  that  several  Born  type  trial  functions  were  used  we  mean  that 
for  a  specific  energy  E  the  relation  27+  K^  =  E,  which  must  be  satisfied 
by  K.  and  ^  in  (K-Kglk-kg)  =  f  2  OJ^  (2.6.2),  permits  some  leeway  in  the 
choice  of  K.  or  K«.   Thus  there  are  two  parameters  that  may  be  varied,  E 
end  K-  ,  say,  in  the  determination  of  the  constants  C  end  D. 

Table  II  represents  the  results  of  numerical  calculations  for  two 
such  Born  type  trial  functions.  The  third  ana  fourth  columns  represent 
results  from  one  of  the  trial  functions  which  are  very  good  for  small 
values  of  the  energy  but  are  poor  for  very  large  values  of  the  energy. 
The  results  listed  in  the  first  two  columns  are  obviously  unsatisfactory. 
The  interesting  feature  of  these  trial  functions  is  that  neither  one  yields 
results  which  approach  either  the  exact  or  the  Born  results  (Table  I)  for 
large  values  of  the  energy,  which  means  that  (2.6.26)  is  not  satisfied, 
yet  one  of  the  trial  functions  yields  ey.cellent  results  for  small  values 
of  E  for  which  the  Born  approximation  is  unreliable. 

Table  III  contains  the  results  of  the  numerical  calculations  for 
those  trial  functions  (symmetric  or  antisymmetric)  for  which  the  Schwinger 
variational  principle  gave  results  independent  of  eny  of  them.   The 
symmetric  trial  functions  give  surprisingly  good  results  for  small  values 
of  the  energy.  Among  these  functions  are  of  course  the  symn.etrized  Born 
type  trial  functions. 
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The  atove  results  emphasize  the  complexity  of  the  problem  of  choosing 
suitable  trial  functions  for  the  Schwinger  variational  principle.  Though 
for  two  body  problems  (see  footnotes  7,  15)  equation  (2.6.25)  seems  to  be 
satisff.ctory,  we  have  seen  that  this  is  not  necessarily  the  case  for  three 
body  problems. 

In  developing  the  methods  of  the  lest  two  sections,  we  noted  certain 
features  that  were  common  to  the  Hulthln-Kchn  and  Schwint;er  variational 
principles.  We  now  turn  to  a   closer  exsjnination  of  both  methods. 
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TABIE  II 
SCHWINGEE  VAUIATIONAL  PRINCIPLE 


%1 


%: 


^1 


- 

.9 

78.05^3 

1.99102 

> 

.5 

7.82721 

1.87292 

— 

.1 

7.76282 

1.64441 

— 

.05 

7.7^562 

1.60664 

— 

.01 

7.72238 

1.57217 

— 

.005 

7.71808 

1.56716 

- 

.001 

7.71661 

1.56271 

+ 

.001 

7.707^5 

35.0382 

1.55937 

.000420 

+ 

.005 

7.69002 

34.7296 

1.55236 

.002095 

+ 

.010 

7.66928 

34.7536 

1.54415 

.004165 

+ 

.050 

7.51826 

34.9560 

1.48623 

.019870 

+ 

.1 

7.3^926 

35.2473 

1.42358 

.037670 

+ 

.5 

6.40765 

38.1054 

1.07304 

.131780 

+ 

.9 

4.66525 

40.6805 

.857346 

.178022 

+ 

5. 

1.02642 

9.40526 

.212231 

.166743 

+ 

10. 

.384106 

12.7160 

.075576 

.094603 

+ 

100. 

.000089 

.230986 

.000292 

.002346 

+ 

1000. 

.000000 

.322939 

.000000 

.000064 

z  =  ratio  of  total  energy  to  energy  of  bound  state. 
(y^   =  elastic  cross  section  for  first  trial  f\inction. 
0  2  ~  elastic  cross  section  for  second  trial  function. 
^^  =  inelastic  cross  section  for  first  trial  function. 
Q.2  =  inelastic  cross  section  for  second  trial  function. 


(A/4tt  =  1) 
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TA3LS  III 

Synnetric  Solution 


Antisytiaetric  Solution 


s 

0 
«8 

^s 

^a 

Qla 

-.9 

1.98668 

.399036 

-.5 

1.89174 

.143718 

-.1 

1.71993 

.114093 

-.05 

1.67722 

.112869 

-.01 

1.6^*642 

.1124/+1 

-.005 

1.64182 

.112470 

-.001 

1.63771 

.112540 

+  .001 

1.63449 

.000440 

.112611 

.0000788 

.005 

1.62749 

.002196 

.112649 

.0003950 

.01 

1.61934 

.00^370 

.112622 

.0007906 

.05 

1.56166 

.020896 

.111948 

.0039010 

.1 

1.49869 

.039700 

.110622 

,00/6240 

.5 

1.14633 

.141436 

.0979432 

.0314150 

.9 

.927302 

.194124 

.0867578 

.0465544 

5. 

.259336 

.211024 

.0374990 

.0641592 

10. 

.107261 

.135296 

.0215910 

.0490688 

100. 

.00349740 

,0101514 

.00241656 

.00799950 

1000. 

.000254130 

.0008708 

.0002436 

.00034758 

z  =  total  ener^jy  in  units  of  Bound  stpte  energy  (e/b  ). 

0   =  elastic  acatterine  cross  section  determined  by  symmetric  trial 
functions, 

Q.  =  inelastic  scattering  cross  section  determined  ty  symmetric  trial 
functions . 

Q^  =  elastic  scntterin^  cross  section  determined  by  ejitisyiMietric  trial 
functions . 

Q^  =  inelastic  scattering  cross  section  determinea  by  antisymmetric  trial 
functions. 
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2.7  Relation  between  the  Hul then-Kohn  and  Schwinger  Vaxiational 
Principles 

The  Hulthen-Kohn  variation?,!  principle  is  based  on  the  functional 

The  notation  employed  here  and  the  significance  of  equation  (2.7.1)  has 
been  explained  in  section  (2.5).   Let 

v(k^k2;  E)  =  f^Ck^kg^  -  [(^oi'^^^/^^oi'^^^J  ^^^  •  (2.7.2) 

in  which 

Vh  =  (l/4TT^)y  (k^kglvlk^kpdk^dk^  h(k^kp 

and 

G7h  =  (l/4TT^)y  G(k^k2:k^kpdk[dk^(k^k^|vlk^«'k  'pdk^^'  h(k«^'ki|), 

for  any  function  h(k^k_).  The  function  g(k-k-)  is  arbitrary  but  the 

functions  f  ,  f  .  and  f  are  those  defined  in  section  (2.5).  For  g   set 
o    01 

equal  to  f  in  (2,7.2)  v  becomes  eq-oal  to  f  so  that  (2.7.2)  can  be  considered 
as  an  admissible  trial  function  for  the  Hulthen-Kohn  variational  principle. 
The  H\ilthfen-Kohn  variational  principle  grve  the  follov/ing  relations, 

^  =  6  f  (K^)i  Z^  =  B^  +  E,  (2.7.^) 


61 
and 


filg  =  6  f  (K^^Kg);  K^  +  K^  =  E  (2.7.5) 

where  f (K  )  and  f (K-Kg)  are  the  asyaiptotic  parts  of  f(K,)  and  f(K,K2) 
defined  by  (2.5.5).  These  asymptotic  parts  can  be  expressed  as  follows, 
f (Z^)  =  -  (f^^,  Vf );  K^  =  B^  +  E  (2.7.6) 

and 

f (K^Kg)  =  -  (fo2'^^^i  Z^  +  ^  =  E  .  (2.7.7) 

The  corresponding  asymptotic  parts  of  (2.7.2)  are  /„  7  a) 

^^^0^  = "  [(^oi'^^'Z  (^oi'^^a  ^^01'^^)  =  -(^01'^^^'  ^  =  ^^*^ 


~  5^' 


and  similarly 

vCK^K^)  =  -  (fo2'^^^'    4  "^  ^  =  ^  •  (2.7.9) 

These  are  exactly  the  same  ns  (2.7.6)  nud  (2.7.7)  respectively.   It 
follov/s  thnt  for  such  trisl  fxinctions  -7 {^z-^c^^i   E)  tKat  equations  (2.7.^) 
end  (2.7.5)  reduce  to 

bl^   =0,   (i  =  1,2)  .  (2.7.10) 

In  other  words  equation  (2.7.2)  defines  triel  functions  which  hsve 
exactly  the  seme  asymptotic  pprts  as  the  exact  wave  function. 
Now  evaluate  (2.7.1)  hy  substituting  (2.7.2)  for  v.   Then 

i/'v^,vj=  i'(^i.e)=  (vjVf^)- l'(foiVf)/(f^.,V£:)j  |(:v.Vg)+(v.,\fGVg)J.  (2.7.11) 

The  tern;  in  curly  'brE.ckets  may  be  written  (v.  .Eg)  where  K  is  the  integral 
operator  defined  by  (2.6.7).  Then  (2.7.11)  becomes 

I'K.g]  =  (^I'V^o^  -  ^^oi'^^^^^'i'^^o^  ^^^v^.Kg}  (2.7.12) 

in  which 

xfv^,  Kg^  =  (v^,Kg)/(Vj,,Vf^)(Vf^^.g). 

\  /"v-iEgV   is  the  q-uantity  in  ter^s  of  which  the  Schwinger  variational 
principle  is  formvilated. 

Consider  the  variation  of  the  functiohal  if '^'it'^'f  E-^out  the  exact 
solutions  f .  and  f .   This  Is  equivalent  to  varying  v  and  g  independently 
nbout  the  exact  solutions.   We  find  that 


61  fv  ,v|         =  -(f^^.Vf)^  Sxfv^.XgJ  (2.7.13) 

This  eq\iation  asserts   thrt  for  independent  variations  about  the  exact 
solutions  the  stctlonnry  character  of  >^  j'»"j.  e  \      implies   the  same  for 
the  ftinctional   llv    ,vlrna  conversely. 

It  must  be  borne   in  mind  that  g  pnd  v  ere  related  by  equation   (2.7.2). 
Because  of   this   the  same  trial  function  used  in  both  variational  principles 
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will  in  general  give  different  approximations.  This  was  fouad  to  be 
the  case  for  a  Born  trial  field. 

late  relationship  between  these  variptional  principles  which  was 
established  above  needs  but  little  alteration  to  be  made  quite  general. 
The  proof  above  depended  on  the  knowledge  of  the  explicit  form  for 
GCk-kpik^kl).  But  this  is  unnecessary  although  It  serves  the  purpose  of 
giving  motivation  to  the  proof.  Start  with  equations  (2.7.1)  and  (2.7,2), 
Equation  (2.7.11)  follows  as  before  with  no  restrictions  on  V  nor  the 
need  for  an  explicit  exp.ression  for  G.  The  remainder  of  the  argument  is 
the  same  and  leads  to  (2.7,13).  The  interpretation  of  this  result  is  as 
before. 

Appendix  III  contains  the  general  proof  outlined  above.  Both  Kohn 
(see  footnote  12)  and  Kato  (see  footnote  ih)   have  obtained  this  relation- 
ship for  two  body  problems. 

Before  proceeding  to  a  discussion  of  a  problem  involvii^  identical 
particles  it  is  desirable  to  summarize  the  results  obtained  thus  far. 

An  exact  solution  to  a  three  body  problem  consisting  of  distir.guish- 
able  interacting  particles  has  been  found.  The  Bom  expansion  for  this 
problem,  which  consists  of  the  Born  approximation  and  all  higher  order 
terms,  can  be  evaluated  and  it  converges  to  the  exact  solution.  However, 
the  Born  approximation  does  not  have  a  useful  domain  of  applicability  for 
any  value  of  the  total  energy  unless  the  interaction  between  the  particles 
is  very  weak  end  the  total  energy  is  large.  Bat  when  the  second  order  term 
in  the  Born  eacpansion  is  included,  a  useftil  domain  of  applicability  is 
found  to  exist  for  large  values  of  the  energy  and  there  is  then  no 
restriction  on  the  magnitude  of  the  interaction  between  the  particles. 
Ihis   result  impl: es  that  in  general  more  terms  in  the  Born  expansion  must 
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■be  evelTxated  "beyond  the  Born  approximetion  if  useful  results  are  to  "be 
obtained.   In  practice  this  is  difficult  to  do.  Two  variational  procedures 
are  then  considered.   The  Hulth^n-Eohn  veriational  principle  limited  in 
application  to  Bom  trial  ftmctions  leads  to  the  same  results  as  the  Born 
approximation.  The  Schwin^er  variational  principle  similarly  restricted 
can  give  excellent  or  very  poor  results  depending  on  the  particular  Bom 
trial  function  used.   Die  problem  here  is  to  supply  some  giiide  in  choosing 
8  good  trial  function.   In  practice  the  use  of  other  than  Born  type  trial 
f\inctions  in  the  Schwinger  variational  principle  leads  to  very  complicated 
mathematical  computations  which  severely  restrict  the  application  of  this 
procedure.  On  the  other  hand,  the  Hul th'fen-Kohn  variational  principle  is 
more  amenable  to  mathematical  calculation  than  the  Schwinger  variational 
principle.  Limited  to  Born  type  trial  functions,  the  potentialities  of 
the  Hulth^n-Kohn  method  were  not  f-ully  exploited.  However,  the  relttion- 
ship  hetv/een  the  Schwinger  and  Eulth^n-Eohn  variational  principles  which 
was  known  to  exist  for  two  body  problems  was  Indicated  by  the  study  to 
hold  for  three  body  problems  and  a  general  proof  of  this  relationship  for 
many  body  problems  resulted. 


Chapter   III 

IDEWIICAL  PARTICLES  -  SYI>C4ETRIC  FEETURmTIQN 

3.1     statement  of  the  Problem 

The   system  which  was   studied  in  Chapter   II  is   now  altered  hy 
permitting  particle  1  to  interact  with  the  center  of  force  and  changing 

the  energy  of  interaction  "between  t}ie  incident  and  bound  particles. 
These  changes   result   in  replacing   the  former  distinguishable  particles 
with  identical  ones. 

The   Schroedinger   equation  v/hich  describes   this   system  is 

(3.1.1) 

t2  2  2  2  "1 

1     bx^  2     dig  J 

All   the   quantities  appearing   in   (3.1.1)   liave   the   same   significance  ascribed 
to   them  as    in  section   (2.1).      The  momentum  representation  is   obtained  in  a 
manner  completely  analogous   to  the  procedure   of  nection    (2.1).      Retaining 
the  notation  of  that   section,   we  find  that    (3.1.1)   becomes 

(3.1.2) 

(k^  +  kg  -E)f(k^k:2)  -   (B/TT)y  f(kj^lC2)dk2  -   {B/tt)J  f(]c^k^)dk^ 

+    (lATT^)y  (k^kgl^lM^P^l'^^  ^^M4^  "  °* 

To  complete  the  mathematical  description  of  this  system,  the  interaction 

kernel  is  assumed  to  be 

(3.1.3) 


(kj^kglvlkj^kp  =  2ATi 


j    ^1  "1 ^   '^2  -^2 I 

j^(k^+T^)  (k^^+T^)    (kg+T^)  (k^+T^)  j  ' 


where  A,  the  interaction  strength,  and  T  are  real  arbitrary  constants. 

The  interaction,  ps  (3-1.3)  will  tie  designated  in  what  follows,  is  symmetric 
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in  the  interchange  of  particles  1  end  2  and  is  also  hermitian,  i.e., 

{k^k^\7\k*]s:^)   =  (k  k  |v|k'k«)  (3.1.^) 

and  ^ 

(k^k^lvlk^kp  =  (k|k^|7|k^k2)  .  (3.1.5) 

In  coordinate  space  this  interaction  is  represented  by  an  integral  operator 
which  hrs  as  its  kernel  the  expression 

A  [6(i2)7(-z^)  exp(Tx^)  5{xpy{^3^)   exp(Ti^) 

-&{x^)^{x^)      exp(-Tx^)  6(x^)>7(x^)  exp(-TX|) 

+  6ix^)y{-x^)   exp(Tx2)  6(x|)^(-xp  exp(Txp 

-6(x^)37(x2)  expC-Txg)  6(x«)^Up  exp(-Tzp3  • 


in  which 

{    0  X  <  0 

y  M   = 

t  1         X  >  0 

This  interaction  is  siailrr  to  that  used  for  the  preceding  problem 

involving  distinguishable  pnrticles.   The  form  chosen  here  is  based  upon 
the  same  considerntions  es  those  given  for  the  choice  of  the  interaction 
(2.1.5)  of  section  (2.1).  This  latter  interaction  cannot  be  used  here 
because  as  we  hsve  seen  the  Schwinger  variational  principle  is  then  in- 
dependent of  symmetric  or  antisymmetric  triPl  functions  pnd  for  identical 
particles  the  solutions  of  physical  interest  are  just  those  with  this  kind 
of  symmetry.   It  is  for  this  reason  that  the  interaction  must  be  altered. 

The  physical  problem  v/hich  is  ma thematicnlly  represented  by  equations 
(3.1.2)  and  (3.1.3)  differs  from  that  considered  in  Chapter  II  by  the 
occurrence  of  exchange  phenomena.   In  other  reepects  the  two  problems  are 
similar.   The  exact  solution  of  (3.1.2)  is  desired  when  the  initial  stete 
of  the  system  describes  sn  incident  particle  from  the  left  with  aomeutum 
K  and  a  particle  bound  to  the  center  of  forcei  Although  the  physically 
interesting  solutions  are  the  symmetric  fna  nntisymmetric  ones,  so  that 

the  Peuli  exclusion  principle  is  satisfied,  it  is  better  for  our  purposes 
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to  otto  in  the  solution  for  the  initial  state  described  above.  Then 
this  solution  can  be  synmetrized.  Once  the  exact  solution  hs-s  been 
determined,  the  various  scattering  cross  sections  are  to  be  calculated. 

The  approximation  methods  are  also  to  be  applied  to  this  problem. 
This  can  be  done  in  two  distinct  vieya.     The   first  of  these  will  be  con- 
sidered in  this  chapter  rjid  consists  in  taking  the  Ipst  term  on  the 
left  hand  side  of  equstion  (3.1.2)  as  the  pertxxrbe.tion.  For  en   obvious 
reason,  this  procedure  will  be  referred  to  es  the  symmetric  perturbation. 
The  second  method  which  will  be  considered  in  detail  in  the  following 
chapter  regards  the  last  two  terms  on  the  left  hand  side  of  (3.1.2)  as 
the  perturbation  and  in  this  form  it  will  be  referred  to  as  the 
asymmetric  perturbation.   It  is  this  method  which  is  generally  found  in 
the  literature  (see  footnote  l). 

The  approximation  methods  will  be  used  to  determine  the  various 
cross  sections  and  these  v/ill  then  be  compered  with  the  exact  values. 

3.2  The  Exact  Solution 

The  direct  method  of  solution  does  not  readily  yield  the  answer  to 
this  problem.  However,  the  expansion  method  and  the  Green's  function 
technique  can  both  be  used.  The  expansion  method  will  be  used  to  deter- 
mine the  Green's  function  nnd  then  the  latter  will  be  used  to  obtain  the 
exact  solution.   In  this  way  the  details  of  both  methods  are  simulteneously 
e3±iibited  and  previous  results  can  be  used. 

The  Green's  function  G(k^k2;  k'kl)  is  that  solution  of 

2  2  />  />  (3.2.1) 

(fc^+k2-E)G(k^k2;k^k|)-(B/TT)/G(k^k2;k|kpdk^-(B/TT)/G(k^k2;k^kpdk2  = 

=  6(k^-k^)6(k2-kp 
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which  represents  outgoing  waves  in  coordinete  space.  The  eigenf unctions 
in  tenns  of  which  the  solution  of  (3.2.1)  is  to  be  expanded  are  the 
solutions  of  the  uin;.ert\u''bed  equation 

(k^+k2-E)f  (k^k2)-(B/TT)/f  (k^k2)dk^-(B/TT)/f(k^k2)dk2  =  0,        (3.2.2) 

which  represent  an  incident  cind  outgoing  wave  in  coordinate  spree. 
Equation  (3.2.2)  is  separable.   Let 

f(k^k2)  =  t[k^)f^{k^)  (3.2.3) 

then  f-  and  f^  satisfy  the  equations, 

(k^-Z^)f^(k^)-(B/TT)C  =  0. 


(k^-i:^)f2(k2)-(B/TT)D  =  0 


(3.2.4) 


where 


C  =/f^(k^)dk^  and  D=y'f2(k2)dk2  (3.2.5) 

and  KT+Kp  =  E,   KT  is  the  energy  of  particle  1  and  Kl  is  the  energy  of 
particle  2.   The  eigenf unctions  of  the  second  of  equations  (3.2.4)  have 
already  been  determined.  These  were  given  in  equations  (2.2.6),  (2.2.7), 
and  (2.2.23).   The  eigenf unctions  of  the  first  of  these  equations  follows 
by  replacing  k_  and  K^  with  k.  and  K^  respectively  in  the  eigenf unctions 
for  the  second  eq\iation. 
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Hence  the  eigenfttnctions  for  equation  (3.2.2)  ere: 

(3.2.6) 

(iLK-kk.)  =   6(k,-IL)  +  ^  -i — 1- ( 


=  B^+S  , 


(Ulk^k^)  =  N^/(k^+B^)(k2+B^),   E  =  -2B^  . 
The  following  identities  follow  from  (3.2.6) 

(K^X2l3^k2)  =  (K2K^|k2k^),  (3.2.7) 

(K^2|k^k2)  =  (iK^lkgk^). 
(lalk^kg)  =   (I3|k2k^). 
The  orthogonality  relations  satisfied  by  the  eigenfunctions  (3.2.6)  are 

y  (K|K^lk^k2)*dk^dk2(K^K2|k^k2)  =  6 (K^-K|)6(K2-K^),  (3.2.8) 

y(K»2lk3_k2)*dk^dk2(K^2|k^k2)  =  6{K^-  Kj^), 
y  (lK2|k^k2)*dk^dk2(lK2|l5:3_k2)  =  S(K2-K^)  . 
y(12lk^k2)*dk^dk2  (12|k^k2)  =  1 
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All  other  products  vanish.  The  closure  property  is  elso  satisfied,  i.e., 

(3.'^.9) 

f  {lK^\]s:*]LplK^{lK^\k^k^)   +  (I2|:c'k-^)*(l2llc^k:2)  =  6i]c^-^)b{k^-k^) . 

Therefore  the  eigenfunctions  (3.2.6)  form  a  complete  set. 

Now  expena  the  Green's  function  in  terms  of  these  eigenf^jjictioiis. 

(3.2.10) 
S(k:^k2;k|kp  =f  G{K^y:^)dK^dK^iK^Z^\k^k^)  +  f  Q{K^)dS^{K^2\k^k^) 

+  y'G(K2)dK2(lK2|k^k2)  +  G(12)  (I2|k^k2) . 

The  functions  0(1''  E„),  G(K.),  G(K_)  end  G(l2)  are  the  ex[jansion  coefficients 

and  are  distinct  functions  which  heve  been  disting:uiahed  for  purposes  of 

brevity  froa  each  other  by  their  arguments . 

The  procedure  for  the  determination  of  these  coefficients  is  exactly 

the  same  as  that  given  in  section  (2.2.2).  The  result  of  the  calculation 

yields 

V  2  2  (3.2.11) 

G(K^K2)  =  (K^K2|k^kp/(K^+K^-E), 


G(K^)   =  (K^2|k^k^)y  (4-K^);  ^  =   B^+E 
G(X2)   =  (iKglk^kpy  (4-Kj)i 
G(12)   =  -(I2|k^kp7  (K^+B^). 


Froai  (3.2.10)  and  (3.2.11)  the  Green's  f\inction  becomes 

(3.2.12) 


rf.    .     ..vn    /'^^2'M4>^^^^2V2"^l'^2> 
GO^k^jk'ky  =J     ^-- 


■/ 


■/ 
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(K^2llc^Jcp*dK^(K^2lk^k2) 

ilK^\)c*k^)*dX2    (lK2lk^k2) 
(I2lk^kp*(l2|k^k2) 

0 

The  identities  (3.2.?)  coupled  with  equation  (3.2.12)  yield 

G(k^k2;k_«kp  =  G(k2k^;k^k^).  (3.2.13) 

Prom  the  orthogonality  and  clostire  relations  we  see  that   the  complex 
conj-ugates  of  the  eigenfunctions  also  form  a  complete  set.     Thios   the 
Green's  function  of  equation    (3,2.2)  could  be  obtained  as  pn  expansion 
in   terms   of   the  complex  conjugates    of   the   eigenfunctions.      Therefore   it 
follows  from    (3.2.12)   that 

G(k^k2;k|k|)  =  G(k^k|;k^k2).  (3.2.1^) 

In  terms   of  the  Green's  fionction,   the  solution  of  equation   (3.1.2) 

may  be  represented  by  the  integral  equation 

(3.2.15) 

f(k^k2)  =  f^(k^k2)-(l/^^yG(k^k2;k|kpdkJdk|(k|k||vlk^»'k2")dk^"dk2"f(k^"k^0, 

where  f  (k,k„)  is  the  solution  of  the  unperturbed  wave  equation  which 
o  1  ^ 

describes  the  initial  strte  of  the  system.  Substitute  the  interaction 
(3.1.3)  into  (3.2.15)  exid  evaluate.  The  result  is 

(3.2.16) 

f(k^k2)=  f^(k3^k2)-(2ATi/UTT2)  [  IiC^+IgC^-I^D^-Ii^Dg]  . 
where 
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(3.2.17) 


^1  "  ^l^V2^  =fGik^k^ik*k^)]^dklak^/   (k^^+T^), 

h   =  ^2 V2)  =/GV2;k.kp:c.dk^dk.  /  (kf+T^). 

I^  =  I3(k^k2)  =/G(k^k2ik'kpdk^dk^y  (k|W), 

(k^k^)  =/G(k^k2;k|kpdk^dk^/  (k^W). 


1^*=  ^il 


and 

°1  =y*^('^k2)dk^dk^(k^+T^);  C2  =y*f(V^2^  ^^2/  (^2*^^) 

(3.2.18) 
D^  =/f(k^k2)k^dk^dk2/(k^+T^);  ^2  =  f  ii,\^2^^2^^Z  / ^^2^'^^^ ' 

The  constrnts  C  ,  Cp,  D  ,  and  D-  cnn  be  determined  from  equations 
(3.2.16)  and  (3.2.18). 

To  simplify  the  calculations  the  following  integral  operators  are 
defined.  For  arhitrery  functions  f (k-kg)  let, 

P^(f )  =/f  (k^k2)dk^dk2  /(kJ+T^),  (3.2.19) 

P2(f)  E/f(k^k2)dk^dV2/(k2  +T^) 

P3(f)  =y*f(k^k2)k^dk^dk2/(k^+T^), 

and 

P4(f)  =y*f(k^k2)k2dk^dk2/(k2+T^). 
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Then  from   (3.2.18)   aad    (3.2.16)  we   o'btaia 

°1  =  ^l^^o^  ""    ^2^1(13)^1  +  ^1(14)^2]    •  (3.2.20) 

^2  =  ^2(^0)*  ifi(V^i-^^i(V^2J  • 

D,   =  Po(f.) 
and 


'1-3-0'-  i[V'i^*  W°2] 


^2=  V^^-i[V'2^*  V'l^'2j 


2 
in  which    ^  =  ZATi/Urr     .     This  siinplificatioii  results  when  the  symmetry 

properties   of   the  Sreen's  function,   equations    (3-2.13)  sJid.    (3.2.1^),  axe 

used  ejid  the  fact  that  ^^{l-^)  =  ^1(12^  "^  ■^3^^3^  ^  ^3^^4^  "^  °* 
The  values   of  the  constants  are 

°1  =  °11  -^  °12   •      ^2  =  \l  -  ^12  (3.2.21) 

where 

(3.2.22) 


2°11 


.p  (3.2.23) 

A^    =!+£''  Pi(i3  -  I4)  P3  (I1-I2). 


2 

and 

^1  =  ^11  -^  ^12      •       ^2  =  ^1  -  ^12   •  (3.2.24) 

where  _^  (3-2.25) 

2D 


=       h^'o^  *  h^'o^  -J2P3(VI2)    fl(^o)-^^2(^o0}  ^1 


2°12  =      1^3(^0^  -  h^'o^  -i^P3CIi-l2)    lfi^fo^-^2^^oM  I  ^2- 

Hence  eqi^ation    (3-2.16)  may  "be  written 

f  (ic^k^)  =  f^dc^k^)  -i  [Cii(Vl2)  ^  ^z^h-h^  (3.2.27) 

-  V^3*V-^2W]  • 
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Prom    (3.2.17)  end    (3-2.13)    it  follows    that 

Igdc^k^)  =   11(^2''^^     ^""^     ^k^h^z'^   "   13(1^2^1)  (3.2.28) 

Since  f  (k^kg)  may  be  written 

(3.2.29) 

fo^^2)  =  I  [^o(^^2^^^o(^2^0'-  2  [^o^'^2K(VlO  • 

eauation  (3.2.27)  represents  the  solution  f  (ic.k2)  as  a  stun  of  a  symmetric 
and  antisymmetric  fiinction.  The  symmetric  function  is 

fgCk^ks^  =  I   [^o(^^2^*^o(^2^1^]  -  ^  [°11^V^2^  -^Il(l3-'l4^]»   (3-2-30) 
and  the  antisymmetric  function  is 

^a^^^2)  =  I  to  V2)  -  'o('^2^l)]  -^['^12^^r^2^  -^2^^3-^4>]  »  ^3.2.31) 

The  initial  state  f  (ic-k^)  has  not  been  specified.   If  f  (k-k^)  is 

synmetric  then  C^2  =  ^12  "  ^  "because  P^^Cf^)  =  ^2^^©^  ^^  ^3^^o^  "^  ^k^^o^' 
Then  for  a  symmetric  initial  state  the  solution  is  just  f  (ic^kp)  given 
by  (3.2.30).   Similarly  if  the  initial  state  f  (k^k-)  is  antisymmetric 
the  solution  f (k^kp)  reauces  to  f  (k^k2). 

The  solutions  of  physical  interest  are  the  symmetric  a nu  antisymmetric 
ones.   In  particular,  our  interest  centers  around  the  scattered  waves 
associated  with  these  solutions.  The  scattered  waves  rre  the  asymptotic 
parts  of  the  solutions  and  these  are  the  following  expressions: 

(3.2.32) 

-1/2 
f(k^k  )  ^(2)    rfc(u-K  )  -N—  ±  6(k„-K  )  ^S_ 
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Ltt(2)^/^     E  -Bi         ^     |k^+T^       T(IH-B)   K -Bi  ^K^+T^ 


N 


^  (k|+B2)(k^-E^) 


Ckg 


r        B  0     _  NnO 

"       tt(2)^  K^-Ei        B       |k?+T^       T(T+B)      K„-E1   (^  K^+T^ 


-j2 


[     \(kJ+T2)(l 


kg  Ik  J 


2.„2, 


)(lk2l-Bi)  (k^+r)(lk^|-Bi) 


-J  T 


1^! 


N 


I)    (k^+B^)(k2-K^) 


(IkJ-BDdkgl-Bl)    ^^kJ+T^ 


,2  „2 
kg+T 


,  2  ,2  „ 
k^+kg-E 


In  this   equption,    the  plus   sign  in  -  refers   to  the   symmetric  solution 
and  the  minus   sign   to  the  putisymmetric   solution.      C  and  D  represent  the 
values  for  the  constants  C      ,   C.^*  ^it.   ^^^  ^i?  wnich  imist  he  associated 
with  these   solutions   and  which  are  given  by  equations    (3.2.30)  and   (3«2.3l)« 
The  initial  state  for  which   (3.2.32)   is   the  solution  is 


f^(k^k2)  =    (2) 


-1/2 


^6(k^- 


(3.2.33) 


V       B      ^0 
K  )  + 

O'  TT 


N 


K  -Bi 
o 


A<hh' 


Z 


t       6(kp-K  )  +  ^ 


N 


2     0' 


B     

""     K-Bi     kf-K^lk^+B^ 
o  c.     oj    1 


The  first  term  on  the  right  hand  side  of  (3»2.32)  represents  that 
portion  of  the  incident  wave  which  is  unaltered  hy  the  interaction  with 
the  "atom.  "  The  second  term  represents  the  elastic  scattering  of 
particles  1  which  is  composed  of  the  direct  scattering  of  the  incident 
particles  1  plus  the  exchajige  scattering  of  the  incident  particles  2  with 
the  bound  prrticlas  1.   The  third  term  represents  the  elastic  scattering 
of  particles  2  and  describes  the  same  physical  phenomena  just  given  for 
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particles  1.   The  fourth  term  on  the  right  hand  side  of  (3.2.32)  represents 
the  inelastic  scattering  of  the  incident  particles  1  and  2. 

Following  the  same  procedure  given  in  Chppter  II  for  distingulshsble 
particles  we  calculate  the  various  cross  sections—elastic  and  inelastic. 
The  elastic  scattering  cross  section  for  particles  1  or  2  is  given  "by 


%   =  (tt'/K^  ) 


|C(Z^)1^  -.  |C(-K^)|^ 


(3.2. 3^) 
(3.2.35) 


where 

B       K       >r„  «  A   CK       D 


0 


°'  "  tt(2)^  Z  -Bi     ^   1"  K^+T^    T(B+T)  K 

^     '  O  V_    0  0 


for  the  syametric  wave  function,   C  (^K  )  for  the  ant  isi'm;ae  trie  solution 
Is  given  hy  the  nejne  expression  except  that  +1,  the  last  term  in 
parenthesis  in  the  curly  hrnckets ,  must  he  replaced  "by  -1.  The  cons  tents 
C  and  D  niust  be  those  thnt  are  to  be  associated  with  the  sy-iimetric  or 
antisymmetric  solutions.  The  inelastic  scettering  cross  section  is 

Q^  =  .{2n^i^/K^)    r  Icl^J^  +  |d12(T-B)^(J2/T^)1  ,  E>0  (3.2.36) 

Here  o  2 


TK?      /   ,      ..2^vn1/2 


"  —     0      ^   1   _   (T  +^)'       )   ,  E^  (3.2.37) 
0  \  b 


'i '  ^W^  (i^.z)V^  [7  '  '    •  ^ 


which  is  the  same  for  both  the  symmetric  and  antisymmetric  solutions  end 

^^     (E.4)l/\^d.,       ,  .2     ^'-''^^ 

J2  =   /    = /—I t  -^ )  .  K>0 


The  +  sign  hfs  the  same  meaning  as  before. 
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Equrtions  (3.2.3?)  and  (3-2.38)  show  that  if  E  =  0,  then  0=0 
whatever  the  values  C  and  D  may  have.  This  of  covorse   is  physically- 
necessary  and  is  similar  to  the  situation  which  arose  in  the  discussion 
of  the  prohlen  for  distinguishable  particles. 

The  continuity  equation  leads  to  the  relation 
Q^  +  0^2  =  k;^  Im  {^^(2)^/^  ^^^0^}  •  (3.2.39) 

This  equ£.tion  is  used  to  check  the  nvunerieal  calculations.   Its  signifi- 
cance has  been  fully  discussed  in  section  (2.3)« 

3.3  13ae  Born  Approximation 

In  analogy  to  the  development  given  in  section  (2.4),  it  is  expedient 
to  consider  the  solution  of  (3.1.2)  as  a  power  series — the  3orn  expansion— 
in  the  interaction  strength  A, 

Let 

00 

f  (k^kg)  =  7"   a""  ^n^^^2^-  (3.3.1) 

n=0 

Sabstitute  (3.3.1)  into  (3.1.2)  and  equate  the  coefficients  of  like 

powers  of  A.  We  obtain  the  following  set  of  equations: 

For  n=0 


(k^+k2-E)f^(k3^k2)-(B/TT)/'f^(k^k2)dk^  -  (B/TT)y  f^(k^k2)dk2  =  0. 


(3.3.2) 


and  for  n  >  1 

(kJ+k|-E)f^(k^k2)  -  (B/TT)yf^(k^k2)dk^  -(B/TT)yf^(k^k2)dk2     (3.3.3) 

=  -  (l/4TT^)yV'dk^dk^  Vl^^^^ 

where 

AV»  =   (k^k2|v|k^k^). 
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Equation  {3.3.2)   represents  the  unperturbed  system  as  e  comparison 
with  (3.2,2)  reveals.  Ihe  physical  situation  descrihed  by  this  equation 
is  th?t  in  which  particles  1  and  2  interact  with  the  center  of  force  but 
not  with  each  other.   'Ihus  the  solution  of  this  equation  represents  the 
initial  state  of  the  system.  Equations  (3.3«3)  represent  the  contributions 
to  the  wave  function  that  arise  from  the  interaction  between  the  particles. 
Physically  the  solutions  of  these  equations  must  represent  outgoing 
particles  in  coordinate  space. 

mie  Born  approximation  consists  in  the  evaluation  of  the  terms 
n  =  0,  1.   This  will  be  considered  below. 

Equations  (3«3«3)  are  solved  by  the  Green's  function  technique.  The 

result  for  n  ^  1   is 

(3.3.^) 
f^dc^k^)  =  -  [  Vi(f,_i)  -   V2(V2)-V3(Vl)-V4(^"n.lO  • 

/I  2 

where  m  =  2Ti/4TT   end  the  notation  of  the  preceding  section  has  been 

used.  From  this  equation  and  the  relations  P  (I^)  =  P  (l„)  =  P_(l„)  = 
P  (I^)  =  0,  we  find  that 

P,(fJ  =  m  [P,(l3)P3(f^_,)-.P,(I,)P^(f^_,)]  . 

P^CfJ  =  n>  [P,(I,)P3(f^_,)-.P,(l3)P^(f^_,)]  . 

P3(fJ  =  -  [^3^h^h^'n-l^   ^  ^3^h^hK.l^]   • 

P^(fJ  =  -m  [P3(l2)P,(f^_,)  ^  P3(li)F2(f,.i)]  . 


(3.3.5) 


and 


These  equations  (3.3.^)  and.  (3.3.5)  niay  be  written  in  matrix  notation 
and  become 

fj^dc^kg)  =  "'^^n-l  "  -'"I'^'^'X  •  (3.3.6) 

\  =  ^Vl  =  '^\   '  (3.3.7) 
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in  which 

I  =  (I,.  I,.  -  I,,  -I.), 


1  n 


E   = 
n 


and 


0        P^(I^)    ^l^^i+^ 


Substitute  (3.3.6)  into  (3.3.I).  Then 

fCk^kg)  =  ^0(1^^2^  "  ^'^^^  -  -^^^'"^  ^0*  (3.3.8) 

The  evaluation  of  this  equation  gives  the  exsct  solution — -equation  (3.2.16). 
Therefore  the  Born  expansion  converges  to  the  exact  solution.  However, 
the  series  representation  of  the  exact  solution  is  valid  only  within  its 
radius  of  convergence.  The  latter  is  determined  hy  those  values  of  the 
elements  of  the  matrix  (l-Ao()  such  that  its  inverse  exists.   In  other  words 
the  smallest  in  ahsolute  value  of  the  roots  of  det.(l-Aa<)  =  0  will  yield 
the  redius  of  convergence. 

Explicitly  the  radius  of  convergence  of  the  Born  expansion  (3.3.1) 
is  given  hy  the  smallest  in  absolute  value  of  the  roots  for  A  of 

where  J(_  =  ZiAT/hrt     .   In  the  notation  of  the  preceding  section 

A  ~  ^  A  1  A?  ^"*  •  ^°  restrictions  have  been  placed  on  the  initial  state 
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f  (k-k-)  in  the  derivation  of  (3.3.8).   The  result  i3  then  applicable 
if  f  (k^lC2)  is  chosen  symmetric  or  antisymmetric.   The  series  representing 
either  of  these  solutions  h.ive  different  radii  of  convergence.  Thus  for 
the  symmetric  solution  P^(I„-  I^)  =  P_(l^-l2)  =  0  ajid  the  radius  of  con- 
vergence is  determined  hy 

^l    =  1  +£^P^(l3-^I^)  P3(l^+l2)  =  0. 

For  the  antisymmetric  solution  it  is  determined  hy  ^    =0.   It  is  apparent 
that  in  general  the  radius  of  convergence  of  a  series  representation  of  a 
general  solution  of  (3.3'l)  is  the  smaller  of  the  raaii  of  converg'-uce  of 
the  symmetric  and  antisymruetric  solutions.   It  follows  that  if  such  a 
solution  is  symmetrized  the  resulting  function  is  li.aited  hy  tjie  smaller 
of  the  two  radii  of  convergence.  On  physical  grounds  one  would  expect  the 
antisymmetric  solution  to  have  a  larger  radius  of  convergence  then  the 
symmetric  solution  because  for  the  former  the  wave  f\inction  is  small  when 
the  interaction  between  the  particles  is  large.   The  net  result  is  to  have 
an  effective  interaction  which  is  weaker  tlian  the  actual  interaction.   If 
we  form  sxi   antisymmetric  solution  from  a  general  solution  (in  series 
representation),  the  resulting  function  may  be  valid  over  a  larger  range 
of  values  than  the  general  solution.   In  general,  when  spin  is  taken  into 
account  both  symmetric  and  antisymmetric  fxonctions  are  needed.   If  reliable 
eSMperimental  data  is  on  hand,  then  a  better  estimation  of  the  effect  of 
the  symmetric  wave  function  can  be  obtained  by  calculating  the  contribution 
from  the  antisymmetric  wave  function  and  subtracting  its  value  from  the 
experimental  data. 
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The  Born  spprozimation  is  given  by 

(3.3.10) 

For  a  symmetric  or  antisymmetric  initial  stft.te  f  (k  k„),  i.e.,  equation 
(3. 2.33) »  the  above  equs.tion  gives  the  following  approximations  to  the 


exact  constants  C^,,  D^.  ,  C^^*  ®^<^  ^12*  ^°^   *^®  symmetric  solution 


IIB 


=  P,  (f  )  =  (2)-^/^  SE 
1^  os'   ^  '      BT 


K  -Bi 
o 


T-B   _^   1 


T4-B 


K^+T^ 


(3.3.11) 


K 


^IIB   *^3^  03''    '''^''      B     -^ 


K^^T^ 


For  the  antisymmetric  solution. 


K 


o 


T-B 


"+1' 


2M-B 


(3.3.12) 


12B 


3  oa 


(2) 


-1/2  Ntt 


K 


^  Z^+T^ 


Here  f   and  f   are  the  symmetric  and  antisymmetric  initial  states  of 
08     oa 

the  system.   These  valiies  for  the  constants  C  and  D  are  equal  to  the 

exact  values  provided  that  A  is  very  small  or  that  the  products 

AP^(I_  +  I^)  and  AP_(I^  ±  I^)  are  very  small.   Biis  situation  is  similar 

to  that  found  for  the  problem  of  distinguishable  particles.   It  implies 

that  the  Born  approximation  is  valid  if  the  interaction  is  weak  and  tlie 

total  energy  is  large. 

The  dependence  of  tlie  radius  of  convergence  on  the  total  energy,  E, 

of  the  system  is 

p  =  T(E)^/2/  TT  (3.3.13) 


for  large  values   of  E.     For  convenience  we  have  taken  p  exx^ressed  in 
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terns  of  |a/4tt|  .  This  relation  holds  for  tjoth  the  symmetrical  and  anti- 
symmetrical  solutions.   It  is  also  identical  with  (2.^.23)  which  was 
obtained  for  distinguishable  particles.   Therefore  the  conclusions  drawn 
there  are  valid  here  also. 

The  general  features  of  this  iteration  procediare  are  the  same  as 
those  discussed  in  section  (2.^)  and  the  results  of  that  discussion  can 
be  carried  over. 

The  Born  approximation  for  nxbitrary  values  of  A  can  readily  be 
considered.  For  large  values  of  the  total  energy  E,  the  values  for  the 
exact  constants  C  and  D  are 

c  =  i  {Z)~^^^^      Ntt/BT(T+B)  (3.3.1^) 

^   r+T^   \     ^"^     T('J>4-B)    / 

Under  the  saae  conditions,  the  scattering  cross  sectioas  become 

q^  =  (8t7T/E)(a/4tt)2  [Cl^  , 

and  2 

(^  =  B  /E  . 

Die   Born  approximation  yields 

C^  =  C,  (3.3.15) 

_  /•o\-l/2     Ntt  0  . 

O 

Since  the  constant  D  plays  no  role  ia  determining  the  scattering  cross 
section,  we  see  that  the  Bom  result  will  be  a  good  approximation  to  the 
exact  result  for  this  range  of  the  energy.  However,  the  scattering 
amplitudes  depend  on  C  and  D  and  in  this  esse  w€  3ee  that  unless 
A/^tt  «  T(14-E)/tt  tiie  Bora  result  will  not  give  good  agreement  with  the 
exact  scattering  amplitudes  even  for  large  values  of  E. 
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The  second  Born  approximation  gives  the  following  values  for  the 
constants  C  and  D  f or  large  E, 

°2B  =  °  •   ^2B  =  ^  •  (3.3.16) 

Here  we  see  that  both  the  scattering  cross  sections  and  scattering 
amplitudes  obtained  from  the  second  Born  approximation  will  give  good 
agreement  with  the  exact  values  for  large  values  of  E.   This  result  is 
again  similar  to  that  obtained  in  the  preceding  chapter  and  the  analysis 
there  holds  in  the  present  case. 

Table  IV  consists  of  the  results  of  a  numerical  calculation  of  the 
scattering  cross  sections  for  the  exact  symmetric  solution  and  for  the 
associated  Born  approximation.  The  parameters  have  been  chosen  so  that 
A/4tt  =  1  and  T  =  2B.  For  this  value  of  the  interaction  strength  A,  the 
table  indicates,  as  is  to  be  expected  from  general  considerations,  that 
the  Born  approximation  yields  useful  results  only  when  the  radius  of 
convergence  of  the  Born  expansion  for  the  symmetric  solution, a      ,  is 
greater  than  1.  For  large  values  of  the  total  energy  E,  the  Born 
approximation  approaches  the  exact  values. 

Table  V  contains  the  same  information  as  Table  IV,  except  that  tiie 
calculations  are  for  the  antisymmetric  solution.   In  this  case  the  value 
assumed  for  A  lies  well  within  the  radius  of  convergence^  and  the 
corresponding  values  for  the  quantitites  calculated  are  very  good  approxi- 
mations to  the  exact  result  for  all  values  of  tiie  energy. 

These  table  s  verify  what  has  been  said  about  the  radii  of  convergence 
of  the  Born  expansion  for  the  symmetric  and  antisymmetric  solutions.   It 
is  also  clear  that  as  the  value  for  A  increases  larger  values  of  E  must 
be  taken  in  order  to  obtain  useful  results  from  the  Bom  approximation. 
The  Born  expansion,  however,  converges  and  here  we  meet  again  the  same 
problem  discussed  in  Chapter  II.  The  analysis  given  there  holds  here  also. 
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To  utilize  the  Born  expansion  pdditional  terms  must  be  calculated  end  in 
prsctice  this  is  very  difficult  to  do.  Variational  procedures  offer  en 
approach  which  suhst-mtially  reduces  the  mathematical  calculations  involved 
and  we  now  turn  to  a  study  of  such  methods. 
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Table  IV 
SIMMETEIC  PEETUKEATIOK  -  STMKSTKIC  SOLUTION 


Exact 


Bom  Approximation 


z 


>S 


is 


Cs 

*e3 

1.067Vf 

1.47820 

.903667 

2.27123 

.925606 

2.32585 

.930891 

2.31255 

.935332 

2.30006 

•93589^ 

2.29841 

.9363^7 

2.29703 

.936573 

2.29633 

.937026 

2.29500 

.93759^ 

2.29326 

.942112 

2.27868 

.9^7773 

2.25885 

.990^3 

2.06239 

1.03888 

1.84770 

1.47143 

.651187 

1.91763 

.289203 

6.01935 

.013138 

20.0416 

.001033 

''IB 


-.9 
-.5 
-.1 

-.05 
-.01 
-.005 
-.001 
+  .001 

+.005 

.01 

.05 

.1 

.5 

.9 
5. 
10. 
100. 
1000. 


.853805 
1.05210 

1.10402 
1.10477 
1.10476 

1.10473 
1.10468 
1.10467 
1.10464 

1.10459 
1.10364 
1.10154 

1.06 137 
1.00185 

.518503 
.284078 

.015933 
.001053 


10-7.147577 

10-6.366396 

10-^.107206 

10-3.286340 

10-2.107717 

.017392 

.037187 

.1065 81 

.107917 

.036397 

.004302 


10-7.175339 

10-6.437171 

IO-5.258I33 

10-3.649928 

10-2.247901 

.429422 

.993868 

.307757 

.257^11 

.041558 

.004359 


z  =  total  energy  E  in  units  of  the  iDOund  state  energy  B  . 

0   =  elastic  scattering  cross  section  for  symmetrical  wave  function. 

0   =  elastic  scattering  cross  section  for  Born  approximation  to 
symmetrical  wave  function. 


Qjo  =  inelastic  scattering  cross  section  for  symmetrical  wave  function. 
LB 


Q.  =  inelastic  scattering  cross  section  for  Born  approximation  to 


symmetrical  wave  function. 


fs  =  li;^ 


radius   of  convergence  of  Born  expansion  for  symmetric  wave 
fionction. 
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Table  7 
SIWMETEIC  PERTDEBATIOH  -  ANTISYMMETRIC  SOLUTION 


Exact  Solution 


Born  Approximation 


z 

Q 

ea 

-.9 

.9^8379 

-.5 

.786795 

-.1 

.678528 

-.05 

.667623 

-.01 

.659256 

-.005 

.658233 

-.001 

.657397 

+  .001 

.657011 

.005 

.656205 

.01 

.655196 

.05 

.6^7282 

.1 

.637722 

.5 

.576029 

.9 

.520376 

5. 

.259111 

10. 

.147364 

100. 

.012500 

1000. 

.001016 

""ia 


10-7.123372 

10"°. 308733 
io"5. 123441 

10"^.  322985 

10-3. 1185 26 

10-2.298766 

10-2.879515 

.0877559 

.II8I74 

.038254 

.004322 


1.75739 
1.78538 
1.82632 
1.83194 
1.83644 
1.83700 

1.83757 
1.83766 
1.83810 
1.83866 
1.84295 
1.84820 
1.87063 
1.92261 
2.28234 
2.67852 
6.46166 
20.2075 


::sE. 


.922751 

.730850 

.630407 

.621099 

.614017 

.613154 

.612465 

.612121 

.611440 

.610589 

.603945 

.596012 

.543805 

.504394 

.290243 

.174831 

.0127781 

.00103251 


!iB_ 


10-^.109669 
10"°. 27425 6 
10-5.109565 
10"7.297264 

10-^.995573 
lO-'-.  24895  7 
10-2.730707 
.074559 
.103738 
.0369564 
.0043067 


z  =  total  energy  E  in  -units  of  the  bovmd  state  energy  B  . 

Q,   =  elastic  scattering  cross  section  for  antisymmetric  wave  f-unction. 

Q  -    elastic  scattering  cross  section  for  Bom  approximation  to 
antisymmetric  wave  fvmction. 

Q,   =  inelastic  scattering  cross  section  for  antisymmetric  wave  function. 
ia 


0   =  inelastic  scattering  cross  section  for  Born  approximation  to 
antisymmetric  wave  function. 

I  A  1 
=  |r~|=  radius  of  convergence  of  Born  expansion  for  antisymmetric 

wave  function. 


P 
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3.^  The  SchwJDger  7ariatioiis.l  Principle 

Consider  the  integral  equation  representation  of  (3. 1.2)  which 

satisfies  the  initial  conditions  of  the  prohlem.   This  is 

(3.^.1) 

f  (k^k^)  =  f^(k^k2)-(l/'M^)yG(k^k2;kj|kpdk^dk^(k^k^lv|k^««k2««)(ik^"di:2"f  (k^«'kj«« ). 

Here  f  (k^kg)  describes  the  initial  state  of  the  system  and  is  a  solution 
of  the  uaperturhed  equation,  i.e.,  equation  (3.1.2)  with  the  interaction 
set  equal  to  zero. 

In  addition  to  (3«^-l)f  we  consider  the  integral  equation  adjoint  to 
it.  This  integral  equation  represents  solutions  of  the  wave  equation  in 
terms  of  incoming  or  concentrating  waves.   It  is  erpresaed  in  terms  of  the 
Green's  ftonction  G^Ck^kpik-'ki)  which  is  related  to  that  occurring  in 
(3.4.1)  as  follows: 

Gg(k^k2;fc»k|)  =  G*(k^k^ik^k2)  (3.4.2) 

This  relation  implies  the  same  properties  for  G-  which  equations  (3«2.13) 
and  (3.2.14)  express  for  G.  The  adjoint  integral  equation  is 

f^(fc^k2)  =  f^^(k^k2)-(l/4Tr^)y6(,(k^k2;k^kpdk|dk^(k^k^|V|k^"k2")<ik^"d^  "^(k^"k2") 

which  may  he  written 

(3.4.3) 

f  ^ (kll^=f oi  (lCi''C2)-(l/4/) 7f  i (k^'«k^)dfc«' 4k2"  (k^"k«  17 [ k^kp*dk|dkp* iXi^^^^^^, 

when  (3.4,2)  is  used  and  the  hermitian  chara-cter  of  the  interaction  is 
taken  into  account.  ^ n^^'^o)  ^^^   solutions  of  the  unperturbed  equation 
which  will  he  specified  below. 

The  integral  equations  (3.4.1)  and  (3.4.3)  are  treated  in  the  same 
way  as  in  Chapter  II,  section  (2.6).  IQie  equations  (2.6.3)  to  (2.6.13) 
inclusive  may  be  taken  over  directly.  The  only  changes  that  mist  be  made 
in  the  application  of  these  equations  to  this  problem  is  the  replacement 


-  78  - 


of  the  interaction  and  Green's  function  with  the  corresponding  quantities 
associated  with  the  present  problem. 

The  Schwinfer  varietional  expresBion  is 


\j^,   -^j   =  (vj^,Kv)/(v^^)(a^,v). 


(3.^.^) 


Now 


(a^.f)  =  2ATi  |C^P3(f^^)  +  C 


[°1^3^<i^  *  ^A<i)-Vl^C)-V2(<i)]'      <3.^.5) 


is  ohtained  from  (2.6.1^)  by  inserting  the  interaction  ana  carrying 
thro-ugh  the  integrations.  The  notation  is  thst  of  section  (3»2).  This 
equation  represents  the  exact  value  for  (a.  ,f ).   The  elements  of  (3.^«'^) 
can  also  be  evaluated  and  the  results  of  the  calcolation  are, 

(3.^.6) 
(v^.Kv)  =  2ATi   P3(v*)P^(v)-P^(v*)P3(v)  +P^(v*)P2(v)-P2(v*)P4(v) 

-(ATi/TT)2|i^(v*)|p^(l3)P^(v)  +  Pi(li,)P4(v)l 

+  P^U*)  f  P3(l^)P^(v)  +  T^{1^)T^{t)\ 
+  Pi^(v*)  |p^(lj^)P3(v)  +  P^(l3)P^(v)  1 

+  P2Uj)  P3(l2^^1^^^  *  Fj{l^)-P^(:r)\  J  , 

(3.^.7) 

2ATi  [Pi(v)P3(<iHP2(v)P4(<i)-r3(v)P,(f;,)-P^(v)P2(f;,)]  . 


(ftjlV) 


and 


(3.^.8) 


(v^.a)  =  2ATi  _P3(v,)Pi(fJ-Pi(v^)P3(f^)*P,,(v^)P2(f^)-P2(v^)P4(f^) 

Equation  (3.'+.^)  gives  an  approximation  to(a.  ,f)  which  in  turn  yields 
approiimntions  to  the  exact  constants  C,  ,  C-,  Di»  sJid  ©2. 
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The  eqtiations  (3.^.6),  (3.^.7)  and  (3.4.8)  are  general  and  no 
restrictions  have  "been  placed  on  the  trial  fiinctions  r   and  v..  The 
principle  is  to  be  used  to  obtain  the  scattering  cross  sections  for 
the  syBimetric  and  antisymmetric  solutions  so  that  the  trial  functions 
T  will  be  chosen  to  have  the  same  symmetry  properties  as  the  solutions. 
There  is  no  a  priori  reason  for  choosing  v  in  this  way.  It  might  be 
argued  that  since  the  exact  solutions  have  these  symmetry  properties  then 
to  obtain  good  approximations  the  trial  functions  should  be  "close  "  to 
the  exact  solutions.   Biis  implies  that  the  trial  functions  should  have 
the  same  symmetry  properties  as  the  exact  solutions.  However,  in  the 
problem  for  distinguishable  particles,  we  found  that  "closeness  ••  was 
an  inadequate  criterion  and  in  fact  found  trial  functions  which  could  not 
be  considered  "close  "  to  the  exact  solutioh  but  which  nonetheless  gave 
excellent  results.  Since  Born  trial  fionctions  of  the  symmetric  or  anti- 
symmetric type  are  to  be  considered  and  since  the  exact  solutions  have 
these  symmetries,  we  shall  restrict  v  in  this  way  but  v.  is  to  be  arbitrary. 

With  this  restriction  on  v,  we  find  that 

(3.4.9) 

(v.,Kv)  =  EAIi  [p^(v)  rP3(v*)  ±  P^(v*)|-P3(v)  |"Pi(v*)  i  ^^i^l)\ 
-(ATi/n)^  ^^(v)  |P3(v*)  ±  P4(''i)]Pi(l3  ±  Ii+) 

+  P^(v)|^P^(v*)  ±  P2(v;)j  P3(I^±  I^)]  . 
(a,.v)  =  2ATi  [Pi(v)[p3(f;;i)  t  F^ifl,)^     -T^iv)!^?^^^)  t  ^z^Ki^}]    ' 


and 


(v^.a)  =  2ATi 


fl(^o^{^3W*^  -  ^^^i^}  -^3^^o){^^^*^  *  ^2^<j\    • 


-  80  - 

The  plus  sign  in  i  is  to  be  associated  with  the  symmetric  trial  function 
V  and  the  minus  sign  with  the  antisymmetric  trial  function.  This  con- 
vention will  "be  adhered  to  in  what  follows,  f  (ic-kg)  is  taken  symmetric 
or  antisymmetric  since  it  represents  the  initial  state  of  the  system. 
Notice  that  the  relations  (3.^.9)  in  conjunction  with  (3.2.19)  resiilt  in 
the  symmetrizatlon  of  v..  Hence  we  may  assumed  that  v.  has  the  same 
symmetry  properties  as  v. 

If  use  is  made  of  (2.6.13).  equations  (3.^.9)  and  (3.^.^)  yield  the 
exact  values  for  the  constants  C. ,  C-.  D, t  and  D,.  For  let  v=  f ,  then 

X-T   v^.fj  =  y(a^,f) 
for  arbitrary  v..  Bie  left  side  of  this  equation  can  be  found  from  (3.^.9) 
and  the  right  side  from  (3.^.5).  This  gives 

P^Cf )  Tji^l)  -   P3(f )Pi(v*)  -d  rP3(f)P3(v*)P^(l3  ±  I^> 

+  P^(f)P^(v*)P3(l^  -  l2^]  =  ^l^^o^^3^^^  -^3(^0)^1(^1^  • 

This  relation  must  "be  valid  for  arbitrary  symmetric  v.  and  hence 

P3(f^)  =  P3(f)  +  iPj^(f)P3(Ij^  ±  Ig)  (3.^.10) 

Pl(fo)  =  Pi(f)  -  lP3(f )Pi(l3  ±  I4) 

For  if     V. (kikg)   ^^  chosen  such  that 
^i(-^'^2)  =  ^i(^^2)' 

Dien  Po(v. )  =  0     and  P^ (v . )     need  not  vanish.      In  this  way  the  first  of 
equations    (3.^.10)  is   obtained  and  then  the  second  one  follows.     From 
(3.2.18)  and    (3.2.19)  it  follows  that     C^  =  Pj^(f )  and  D^  =  P3(f ). 
Solviiig    (3.^.10)  for  C     and  D^ ,  we  obtain  the  exact  values  for   tiies© 
c  ons  tants , 
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The  restilts  for  C  and  D  obtained  by  using  Born  type  trial  fvmctions 


are, 


(a)  for  symmetric  solution 


C.      =  D. 


1^0^ 


V3^^^     _2Bi_     _1 


K  -Bi  K  -Bi 


1  ^A(l3«,)  ^\ 


K    -F3(f„) 


P,  (f  ) 


3^o 


-1 


(3.^.11) 


/I  K  +Bi 
=     1-i     ° 


Bi 


Po(fJ 


P,  (f ^)     K  -Bi 


^^^3  *  h^  pSt)  *  ^3^'i-^'2)  Fin  T^ 

_  X     0  J    o         0 


where         ^  =  2ATi/im^  ; 

(b)  for  antisymmetric  solution 


^"o) 


-T    -1 


P,(f    ) 


p,^f„;  1    r  ^^ 


(3.^.12) 


A> 


=    1 


,Z  +Bi 
^   2Bi 


F3(fQ)  Pj^Cf^)        K^-Bi 

P^C  13-14)  p  (f  )    +  ^3(V^2^P-(f  )      kTbT 

10  J     0  0 


If  (3.^.11)  and  (3.^.12)  are  expanded  to  the  first  order  in  A  we  obtain, 
C^  =  \if^)   +  -^Pi(l3  -   I4)  P3(fo)  .  (3.^.13) 

where  f  must  be  replaced  In  the  above  equations  by  the  symmetric  or  anti- 
0 

symmetric  initial  st?te  in  conformity  with  the  convention  on  ^  »  Comparing 
(3.^.13)  with  (3.2.22),  (3.2.23).  (3.2.25)  and  (3.2.26),  we  see  that  to 
the  first  power  in  A  these  equations  are  the  same.  This  is  the  second  Born 
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approximation  and  is  an  instance  of  the  general  reeult  olitained  in  section 
(2.6)  and  discussed  below  equation  (2.6*26). 

The  validity  of  the  expansion  which  led  to  (3»^»13)  can  be  investigated 
for  Iprge  values  of  E.  The  values  for  the  constants  become 

\  -1 


_      CoN-l/2    Ntt       ^o   /      A     TT 


and 

D_  =   (2)"^/^  Ntt  K^B(E^+T^). 


-1 


(3.^.14) 


(3.^.15) 


Compsring  these  results  with  the  exact  values  for  C  and  D  in  the  seme 
range  of  energy  values,  equation  (3«3»l^)t  we  find  that  only  if 
A/^tt  <5C  T(I4-B)/tt  can  the  results  be  said  to  be  in  good  agreement.   In 
other  words  the  Schwinger  variational  principle  for  the  Born  type  trial 
functions  chosen  is  valid  provided  the  interaction  strength  A  is  small 
and  the  energy  E  large.  Tot   the  calculation  of  the  scattering  cross 
sections  in  the  case  of  the  symmetric  solution  this  result  is  less 
satisfactory  than  the  Born  result  (3«3»15).  The  antisymmetric  solution 
(3. ^.15)  is  the  same  as  that  obtained  from  the  Born  approximation,  (3. 3»15)« 

Table  VI  represents  the  results  of  the  calculations  for  the  scatter- 
ing cross  sections  for  the  values  of  the  parameters  a/4ti  =  1  and  T  =  2B. 
The  results  obtained  for  the  antisymmetrical  wave  function,  0  „  oJid  Q.  , 
are  interesting.  Table  Y   shows  that  for  this  case  the  Bom  approximation 
gives  results  in  good  agreement  with  the  exact  ones  over  the  entire  energy 
range.  Except  for  very  large  values  of  the  energy,  the  Schwinger  variational 


-  83- 

principle  does  not  f;iTe  useful  results.  Here  we  have  an  example  of  the 
point  stressed  in  section  (2.6).  Altho^ogh  there  exists  values  of  the 
energy  in  which  hoth  the  Born  results  and  those  obtained  from  the  variap- 
tionel  principle  ere  in  good  agreement  yet  the  variational  principle  fails 
to  extend  the  usefulness  of  tliese  results  beyond  this  common  domain  of 
applicability.  As  a  matter  of  fact  here  the  Born  approximation  is  valid 
over  the  entire  range  of  the  energy. 

The  results  of  the  calculations  for  the  scattering  cross  sections, 
0       and  Q^  ,  for  the  symmetric  solution  ere  fax  from  satisfectory  for  both 
the  Born  epproxiniation  and  the  Schwinger  variational  principle.  For  large 
values  of  the  energy,  the  inelastic  scattering  cross  setttion  is  sensitive 
to  the  values  of  the  constant  C  whereas  the  elastic  scattering  is  inde- 
pendent of  the  values  of  the  constants  C  end  D.  This  is  clear  from  the 
expres 'lions  for  the  scattering  coefficients  given  below  equation  (3.3. 1^). 
Therefore  the  fact,  that  Q.  does  not  approach  the  exact  value  for  large 
values  of  the  energy  but  that  the  Born  result  does ,  implies  that  the 
expansion  (2.6.25)  which  relates  the  Schwinger  results  to  the  Born  results 
is  invalid  ?nd  that  the  condition  (2.6.26)  is  not  satisfied.   This  can  be 
seen  in  greater  detail  if  equations  (3«3.1^)»  (3.3«15)  aad  (3.^.1^)  are 
compared.   It  is  then  evident  that  when  A/^tt  <3c;  T(m-B)/Tr  that  tlie  exact. 
Born  and  Schwinger  results  all  agree  for  large  values  of  E.  However,  if 
this  relation  is  not  satisfied  then  the  Schwinger  results  for  the  scattering 
cross  sections  and  acattertng  amplitudes  will  not  be  good  approximations 
to  the  exact  values  and  in  fact  will  give  poorer  results  than  the  Bom 
results.  Using  Born  trial  functions  in  the  Schwinger  variational  principle 
gives  no  assurance  of  obtaining  usefizl  results  or  even  of  obtaining  an 
improvement  over  the  Born  results . 
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Table  VI 

SCHWINGER  VAEIATICWAL  PRINCIPLE 


Symmetric  Solution 


Antisymmetric  Solution 


z 


-.9 
—-5 
-.1 

-.05 

-.01 

-.005 

-.001 
•f.OOl 

+.005 

+  .01 

+  .05 

.1 

.5 

.9 
5. 

10. 

100. 

1000. 


\i 


%. 


6.91562 
15.^52 
27.5780 
29.5696 
31.3452 

31.4602 

31.6319 
31. 7186 
31.8932 
32.1102 
33.8527 
36.0211 
57.4091 
52.5209 
5.68760 
1.50297 
.0279385 
io"2. 113839 


10"5.380808 

lo-^;.  954174 
10*3.529686 
.0139936 
.0569853 
1.65236 

3.86133 
3.56132 
1.85602 
.167764 
.0192152 


.875615 
.559420 
.394678 

.379955 
.368888 
.367546 
.366484 

,365949 
.364890 
.363594 

.353365 
.341343 
.268081 
.221184 
.327190 
.194063 

.0139329 
.00121054 


10*^.362062 
10-7.907172 

io"f .363278 
10"7 .105366 
10-^.328266 
10-2.101452 
W^.  347454 
10- J. 525699 
io"i. 749875 

10"  •'•.340824 
10-2.426691 


z  =  total  energy  in  units  of  the  bound  state  energy  =  S/B^. 

Q   =  elastic  scattering  cross  section  for  symmetric  solution. 
eS 

Q,   =  inelastic  scattering  cross  section  for  symmetric  solution. 
is 

Q,   =  elastic  scattering  cross  section  for  ant isyimne trie  solution. 


Q,   =  inelastic  scattering  cross  section  for  antisymmetric  solution. 
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The  work  in  this  chapter  has  elements  in  common  with  that  of  the 
preceding  chapter.  The   Born  expansions  in  both  have  been  found  to  converge 
to  the  exact  solution.  Por  large  vsluee  of  the  energy  their  radii  of  con- 
vergence are  also  the  seme.  Furthermore  the  condition  A/Att  «  T(!W-B)/tt 
which  when  satisfied  assures  good  agreement  among  the  exact.  Born  and 
Schwinger  results  for  the  scattering  amplitudes  and  cross  sections  is  the 
seme  in  both  problems.   Ihese  similarities  impose  the  same  limitations  on 
the  approximations  as  discussed  at  some  length  in  the  preceding  chapter. 
Bie  symmetric  problem,  however,  gives  rise  to  other  results.   These  are 
best  understood  in  terms  of  the  symmetric  and  entisymmetric  solutions  of 
the  problem.  ISiuB  the  Born  expansion  for  both  these  solutions  converges 
to  the  exact  solution  but  have  different  radii  of  convergence  which  approach 
each  other  for  large  values  of  the  energy.  The  antisymmetric  solution,  as 
one  would  expect  on  physical  grounds,  converges  more  rapidly  than  the 
symmetric  solution.  The  Bom  approximation  here  in  contrast  with  that  for 
the  problem  of  distinguishable  particles  has  a  domain  of  applicability. 
This  domain  is  larger  for  the  antisymmetric  solution  than  for  the  symmetric 
solution.  This  is  to  be  expected  from  the  fact  that  the  Born  expansion  for 
the  former  converges  much  more  rapidly  than  for  the  latter.  The  behavior 
of  the  Schwinger  variational  principle  when  applied  to  Born  type  trial 
functions  is  such  that  though  the  theoretical  relation  (2,6,25)  between, 
the  second  Born  approximation  and  Schwinger  res-ult  holds  for  large  values 
of  the  total  energy  in  1iie  case  of  the  antisymmetric  solution,  still  the 
variational  principle  does  not  extend  the  result  beyond  this  commoii  domain, 
Per  this  same  solution  the  Bom  approximation  gives  excellent  results 
through*ut  the  range  of  energies,  Obe  important  characteristic  of  the 
Bom  expansion  for  the  symmetie  problem  with  the  symmetric  perttirbation 
is  that  the  effects  due  to  exchange  appear  in  all  terms  of  the  expansion 
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end  in  particular  in  the  Bom  approximation. 

The  symmetric  perturbation  considered  above  is  not  readily  applicable 
to  actual  physical  problems  because  it  leads  to  mathematical  quantities 
which  are  extremely  difficult  to  evaluate.  For  this  reason  the  problem 
is  usually  formulated  in  terms  of  an  asymmetric  perturbation.  We  now 
turn  to  a  consideration  of  this  problem  described  in  this  way. 


Chepter  lY 
Identical  Particles  -  Asymmetric  Perturbation 

In  general  many  body  problems,  for  ei£>iaple,  the  interaction  of  an 
electron  with  s  hydrogen  atom,  are  not  mathematically  tractable  when 
considered  from  the  point  of  view  of  the  symmetric  perturbation.  However 
when  these  problems  are  formulated  in  terms  of  en  asymmetric  perturbation 
the  mathematical  difficulties  are  considerably  reduced.  Since  this  letter 
procedure  is  the  one  most  often  found  in  the  literature,  it  is  important 
to  apply  it  to  the  problem  considered  in  the  previo-us  chapter. 

There  are  some  difficulties  thft  center  around  the  calculation  of 
exchange  scattering  when  this  problem  is  treated  by  approximation  methods. 
For  example,  there  is  some  question   whether  en  iterative  procedure,  the 
Born  expajasion,  can  account  for  exchange  scattering.  A  method  due  to 

Mott  ejid  Massey  (see  footnote  ih)   which  wes  developed  for  this  purpose 

^1  kZ 
has  recently  been  subject  to  study  and  criticism  *   (see  also  footnote 

10).  Both  of  these  procedures  will  be  studied  in  some  detail.   In  addition 

the  method  given  by  Borowitz  and  Friedman  (see  footnote  10)  will  be 

considered.  The  latter  extends  the  Schwinger  variational  procedure  to 

three  body  problems.   In  what  follows, these  methods  will  be  studied  for 

their  ability  to  handle  exchange  scattering. 


kO,     W.  Kohn,  Private  Communication. 

kl.      S.  Altshuler,  Phys.  Eev.  21.  1167  (1953). 

hZ.     H,  E.  Moses,  Phys.  P.ev.  £1,  185  (1953). 
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4,1  An  Iterative  Procedure 

The  Schroedinger  equation  may  "be  written  in  the  following  way, 

-  -^     /(k,k2|7lk'k')dk^dk|f(k'k«). 

The  left  hand  side  of  this  equation  equated  to  zero  is  the  -unperturbed 
equation.   It  represents  the  physical  situation  in  which  particle  2 
interacts  with  the  center  of  force  hut  particle  1  neither  interacts  with 
particle  2  nor  with  the  center  of  force.  The  solution  of  the  unperturbed 
equation  is 

f^(k^k2)  =  6(k^  -  E^)  N/(k|  +  B^),     -       (4.1.2) 

and  describes  the  initial  state  of  the  system  in  which  particle  1  is 
incident  from  the  left  with  momentum  K  and  particle  2  is  bo\md. 

0 

Equation  (4,1.1)  is  to  be  solved  by  inserting  (4.1.2)  in  place  of 
f(k^k2)  on  the  right  side  of  (4.1.1).  The   solution,  f  ^■''^  (k^k2),  of  the 

resulting  equation  is  then  used  to  replace  f (k^k^)  on  the  right  hand 

(2) 
side  of  (4.1.1).  This  determines  the  function  f^  '^(k^k^)  and  the  above 

process  can  be  repeated  so  that  f   (k^kg),  n  >  1  is  a  solution  of 

The  solution  of  this  equation  is  /^  ^    r» 

f^''^(k^k2)  =  ^0(^1^2^  *  |/^^^1^2'M4^^4'^4^^''"'^^^^l"*'^2  ^*^^l' 

-(2ATi/4TT^)  IP^(G)P^(f^*»-^b-P^(G)P^(f^^^')+P^(G)P2(f^^-^^) 


-  P2(G)P^(f^^-^^)  ] 
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in  which  the  Green's  function  for  the  unperturbed  equation  is  used  and 
its  explicit  form  is  given  hy  (2.2.37)  or  (2.2.38).   The  notetion  of 
section  (3.2)  has  "been  used.  Note  that  P.(f  "^  ),  i  =  1,2,3,4  are 
constants  whereas  the  terms  involving  the  Green's  function  are  fiinctions 
of  k^  BJid  kp  for  the  integral  operators  P.  will  be  understood  to  operete 
only  on  the  primed  variables  of  the  Green's  function. 
The  evalur.tion  of  P.  (G)  yields 

1     Ve^        1 

^1^^^  =   22    -^== .  (^.1.5) 

K7+T      /   2  -        2  2 
^L       TE-k^    Bl    k^+k^-E 


S  1 + 

T 


2  „2 


( VE-k^  -Bi )  ( VE-kJ  +Ti  )      kg+T 


T^(G)   = 


P-(G)  =   — p — ^ —   7-   ■>■ -5 — 5 .  and 

^  (k^+T  )   -yE-k"  -Bi    k^+k^-S 


^^^^^  "  ,2,™2     ,2.,  2 


2  2     ' 

k^+k^-E 


Prom  equations  (i*.1.4)  and  (4.1,5)  it  is  clear  that  exchange  terms  can 
only  arise  from  the  second  member  en  the  right  hand  side  of  (4.1.4), 
No  other  member  contains  a  singularity  on  the  real  axis  of  ^2   which  at 
the  same  time  is  associated  with  a  bound  state  in  particle  1.  The 
initial  state  has  particle  2  bound  and  particle  1  incident  so  that  we 
seek  a  state  in  which  particle  1  is  bound  and  particle  2  is  outgoing. 
Such  a  state  indicates  that  particles  1  and  2  have  been  exchanged.   If 
exciiange  has  talcen  place  it  can  only  arise  from  the  second  member,  the 
other  members  in  (4.1,5)  do  not  contain  the  singularities  appropriate  for 
an  exchange  state.  Thus  consider 
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f  /^(Va'M^^'^i^a  f^'^"'^ki"4^i"  = 


(4.1.6) 


B 

TT 


2  2 


B 


"^  -/i^-  Bi 


,(n-l) 


(n-1), 


(k^)  +  D2--^'(k2) 


in  which 


and 


B^'^'fl^) 


D<-^'(k,) 


-=s 


D 


(n-1) 


(Icpdk^ 


2   2 


(4.1.7) 


.  /  f^-^)(k,.2^dk^ 


In  (4.1,6)  the  only  term  capable  of  giving  rise  to  the  proper  singularity 

.(n-1), 


on  the  real  axis  of  kg  ^^  ^®  term  containing  D^    (^2^* 
D^-)(k2)=/dlc,f(-)(k,k2) 


(4.1.8) 


-L+       Bi         D(n-l)(k  )  _   (2ATi/4TT2)p(^-l)(k  ) 


(!)7: 


dk  -/E-k^ 


2.,  2 


kT+kg-E       yE-k7  -  Bi 


2 


D{"-'^k^). 


where 


^^kg)  =/dk^[p3(G)P^(f^'"-^^)-P^((j)P^(f^''-^^)+P^(G)P2(f^''"^^ 

-P3(G)P^(f^°-^^^ 


(4.1.9) 


If  now  we  let  n  ->  oo,    then,   in  an  obvious  notation, 


D^Clc^)  = 


V5 


E-k~ 


Ys-kg     -Bi 


N 


k^.B^ 


/      x2    ^      dk^        Ve^ 


(4.1.10) 


k^+k|-E   Vs-k^    -Bi 


I'^^L' 


-(2ATi/4TT  )P(k,) 


] 
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The  singularity  on  the  real  axis  of  k^  is  evident  in  this  equation.   It 
is  given  "by  the  factor  outside  the  square  braxikets.  However  this  singularity 
only  appears  in  the  limit  as  n.  -J>  oo.  For  finite  n.  (4.1.8)  does  not 
exhihit  this  singularity.  Hence  this  iterative  procedure  when  represented 
"by  a  series  does  not  contain  exchange  terms.   Only  when  the  analytic 
function  represented  by  the  series  is  determined  are  the  exchange  terms 
properly  accounted  for. 

Equation  (4,1.8)  may  be  looked  upon  as  e  recurrence  relation  for 
Dp(kp),  Employing  it  in  this  way  we  findlhat  the  leading  term  is 


N 


1  + 


YE-k2   yi^4) 


The  series   in  the  brackets   is  the  es^jansion  of 


Ve-2 


VE-kf  - 


Bl 


1  -(Bi/yi-kl  ) 


-^   . 


But 


and  this  expansion  is  invalid  at  the  singijlprity,  i.e.,  k^ 
this  singularity  is  just  the  one  needed  to  obtain  the  scattered  wave. 
This  is  the  reason  that  the  iterative  procedure  fails  to  give  exchejoge 
scattering.   The  Born  expansion  breaks  down  at  the  singularities  of  the 
wave  function.  The  series  does  not  converge  to  a  solution  which  contains 
exchange  terms. 

Let  the  interaction  strength  A  be  set  equal  to  zero.   Then  the 
iteration  procedure  yields 

n 


K-o)^^[-^-"H^^y 


V  0 


N 


2  2  * 

k^+B 
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which  in  the  limit  as  n  -v*  oo  is 


6(k-X  )  +^     ° 


'l  0     "  K-Bi  Ic5-E^ 

0        i   0 


N 


2  2 

k2+B 


This  is  the  exact  solution  of  the  -onpert-urbed  symmetrical  proljlem.  Not 
only  does  the  iteration  procedure  in  the  limit  account  for  exchange 
scattering  but  also  gives  the  correct  initial  state  of  the  system  which 
is  obtained  when  the  problem  is  considered  to  be  the  unperturbed  symmetrical 
one.  This  latter  result  is  important  as  we  shall  see  when  the  Mott  and 
Massey  procedure  is  discussed  below. 

It  must  be  emphasized  that  only  when  the  analytic  function  represented 
by  the  series  is  determined  do  we  get  the  exact  result.  Otherwise  we 
obtain  an  approximation  to  the  exact  initial  stste  and  excliange  scattering 
is  not  included. 

Tlie  Born  expansion  (3. 3.1)  obtained  from  the  symmetric  perturbation 
is  the  result  of  an  iterative  procedure  similar  to  the  one  discussed 
above.   In  section  (3.3)  it  was  foiind  that  the  exchange  terms  were  present 
in  ell  orders  of  approximation.  The  reason  for  this  can  be  found  by 
comparing  the  Green's  functions  for  these  two  iterative  procediires.   In 
the  case  of  the  sym:aetrical  perturbation,  the  associated  Green's  ftmction, 
(3.2.12),  contains  ab  initio  all  the  singularities  that  contribute  to  the 
scattered  wave.  However  when  the  Green's  function  (2.2.37)  associated 
with  the  asymmetrical  perturbation  is  examined  we  find  that  not  all  the 
singularities  of  interest  are  present.  Furthermore  the  use  of  the  initial 
state  in  the  iterative  procedure  fails  to  supply  the  missing  singularities 
30  that  if  the  iterative  procedure  is  to  converge  to  the  correct  solution 
it  must  supply  the  missing  singularities.   The  individual  ter:ns  of  the 


-93- 


series  are  free  from  these  singularities.  Therefore  the  siun,  i.e., 
the  analytic  function  represented  by  the  series,  must  contain  them. 
Mathematically  this  means  that  the  series  is  not  uniformly  convergent. 

The  iteration  procedure  evaluated  to  the  second  order  of  approximation 
yields  the  following  asymptotic  part  for  the  approximating  wave  function. 


f  ^^^k^k2)  ^  6(ki-^o)  N/(k2+B2) 


gJAT  Ntt 


-A{l-H^-\-)] 


T(BfT) 


F 


2  2 


kf-K^ 
1  0 


(4.1.11) 


2iAT 
4tt 


.  2 


Cl^illCjl 


ifc^|/    (kj+i^jdkji  -Bi) 


2   2 

k^  +T^ 


BD, 


Ikgl-Bi 


I  kg!  -i^ 


(k?+T2)(lkJ-Bi) 


2  2 

k^+k^-E 


When  this  result  is  compared  with  the  asymptotic  pert  of  the  exact  non- 
symmetrized  solution  we  find  that  the  exchange  term  is  absent  as  we  expect 
from  the  above  general  discussion  but  that  the  other  terms  in  (4.1.11) 
are  the  expansion  of  the  exact  terms  up  to  the  first  order  in  B. 

A  procedure  which  takes  into  accotint  all  the  singularities  oY  interest 
will  now  be  considered. 


4,2  Mott  and  Massey  Procedure 

The  Schroedinger  equation  (3.1.2)  is  written 

(4.2.1) 
(k^+k2-E)f(k^k2)-(B/TT)y*f(k^k2)dk2  =J'ytik^^^i]s:^kpdk*d}c^fik*)s:p 

where 

w(k^k2;k^kp  =  (B/TT)6(k2-kp  -  (i/4tt^)  (k^k2l  VJk^kp  .  (4.2.2) 

The  unperturbed  equation  is  the  left  hand  side  of  (4.2.1)  set  eqiml  to 
zero.   This  eqviation  in  a  slightly  different  notation  has  "been  considered 
in  (2.1)  and  (2.2).   It  is  separable.  Let  f(k^k2)  =  f^(k^)f2(k2).  Then 
^2^^2^  must  satisfy 

(k|  -k|)  f2(k2)  -  iB/v)ff^{k^)ik^   =  0.  (^.2.3) 

The  eigenf unctions  of   this  eqijation  have  been  found  and  appear  in   (2.2,23) 
and   (2.2.6).     Explicitly,    they  are 

(2|  kg)  =     N/(k2  "  ^^^'     ^  =    (2/tt)^/^  B^/^  ,  E  =  -B^  (4.2.4) 

B        1^2'  1 

izJ  fc,)  =  6(k^-K2)  +  J -p~ ;  c=(>  0.  E  >  0. 

^  ^  2  "2         TT      j^^|_3^     kg-^-^"^ 

The  orthof-onality  relations  for   these  eigenfunctions  are 

y  (2|k2)*  dk2(2lk2)  =  1.  (4.2.5) 

/(2lk2)*dk2(K2|k2)  =  0. 
/(K'|k2)*dk2(K2lk2)  =   6(E2-Kp. 

Biey  also  satisfy  the  closure  relation,  i.«., 

/  {K^\]c^)*dK^{Y.^\k^)   +  (2|kp*(2|k2)  =  5(^2-^]^  •  (4.2.6) 
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Therefore  these  eigenfunctions  form  a  complete  aet. 

The  solution  f  (ic-kp)  of  (4.2.1)  considered  as  a  function  of  k_ 
may  be  expanded  in  terms  of  these  eigenfunctions. 

where  f. (k^2)  and  fgCk^Kg)  are  the  expansion  coefficients  which  are  to 

he  determined  as  fvuictions  of  k^ .   Insert  (4.2.7)  into  (4.2.1).  We  obtain 

(k^  -K^)  f^(k^2)(2|k2)  +/(k^-K^)  f2(k3^K2)dZ2(K2|k2)  (4.2.3) 

=  y*  w  (k^k2 ;  k^kp  dk^dk^f  (k^k| ) , 

where  K^  =  B  +E  and  KT  =  E-S^  .  Multiply  this  equation  by  (2|k2)*dk2 
and  (Kilk-)  dkg  in  turn  and  integrate  over  the  indicated  variable.  The 
following  two  equE.tions  emerge, 

(kJ-K^)f^(k^2)  =/  (2|k2)*dk2w(k^k2;k|kpdk^dk^f(k^kp,         (4.2.9) 

and 

(k^-.K^)f2(k^K2^  =y  (K2|k2)*dk2w(k^k2;k|k^)dk^dk^f(k^kp.        (4.2.10) 

The   initial  state  of  the  system  describes  particle  1  incident  from  the 

left  with  momentum  K  and  particle  2  bound.  Hence  the  solutions  of 

(4.2.9)  and  (4.2.10)  are 

(4.2.11) 
f^(lc^2)  =  6(k^-K^)  +  -^   y  (2|lc2)*<ik2w(k^k2ik^kpdk^dk^f(k|kp, 

and 

f2(k^K2)  =  -2~F  J    (^2\^2^*^2^^h^2''H^2'^i^2^^HH^'  (^.2.12) 
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In  both  of  these  equations  the  singularities  on  the  real  k^  azis  are 

explicitly  represented.  These  equations  are  immediately  applicable  to 

the  calculation  of  the  elastically  and  inelastically  scattered  waves. 

An  explicit  representation  of  the  singularities  on  the  real  k^  axis 

is  desirable.  This  csn  be  obtained  by  rewriting  equation  (4.2.1)  as 

follows, 

(4.2.13) 

(kJ+k2-E)f(k^fc2)  -  (B/TT)yf  (k^k2)dk^  =yv(k^k2;k|k|)dk^dk|f  (k^kp 

where 
▼  (k^kgiki^kp  =  (B/TT)6(k^-k^)  -  (l/4TT^)(k^k2|V|k|kp.  (4.2.14) 

The  left  hand  side  of  (4.2.13)  equated  to  zero  is  now  considered.  It  is 
treated  in  a  manner  completely  analogous  to  that  given  equstion  (4.2, 3) • 
The  complete  set  of  eigenftuxctions  which  are  obtained  are, 

(l|k^)  =  N/(kf  +  B^)  ;   N  =  (2/it)^/2  B^^^  (4.2.15) 

B    1^1  ^ 


These  eigenfunctions  satisfy  similar  orthogonality  and  closure  relations 
as  do  the  eigenfunctions  (4.2.4)  and  form  a  complete  set  for  fvuactions 
which  depend  on  k, . 

Consider  the  solution  f (k^k^)  of  (4.2.1)  as  s   function  of  k,  and 
expand  it  in  terms  of  the  above  eigenfunctions.  Then 

f(k^k2)  =  g^(lk2)(l|k^)  +  J^  g^{K^k^)aX^(K^\k^) ,  (4.2.16) 


in  which  g  (ik^)  and  g2(K^k2)  are  the  coefficients  of  the  expansion  and 
are  functions  of  k^  which  are  to  be  determined.  Following  the  procedure 
used  in  the  evaluation  of  f. (k-2)  and  f2(k^K2),  we  find  that 
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^1^^2^  "  -~2  y  (l|k:^)*dk^v(k^k2;k^k|)dkL«dk^f(lc'kp.        (^.2.1?) 


kZ-K' 
2  0 


and 


ggCK^kg)  =  -—2"  /(^l'^)*^l^(^^2''^l^i^*^'4^F^4^2^*       (^.2.18) 

Here  ^  =  E-ET  and  K  =  B  +  E.  These  equations  do  not  contain  any 
inhomogeneouB  terms  since  the  physical  sitiiation  requires  that  particles  2 
mast  "be  outgoing.   In  these  equations  the  explicit  dependence  f (k-k^)  <^ 
the  singularities  on  the  real  axis  of  k2  is  exhi'bited. 

Equations  (A. 2.11),  (4.2.12),  (^1.2.1?)  and  (4.2.18)  are  exact.   To 
solve  them  approximately,  f (k^kg)  in  the  integrands  on  the  right  hand  side 
of  these  equations  is  replaced  "by  the  initial  state  of  the  system  which 
is  given  "by 

f^(k^k2)  =  6(k^-K^)  N/(k|  +  B^)  (4.2.19) 

The  results   of  the  calculrtions  are 


f^(k^2)  =   6(k^-K^)  +  -2^ 

V"o 


B  p/NTT 

TT       "*Ib 


k-,-X 

1       0 


^^   '    CkJ+T^)(K^+T^) 


(4.2.20) 


and 


g^dkg)  -     2^^ 


2     o 


K^+B 
o 


si_    of^f     1 

:^g2     "•^l^B  J       T(T+3) 


r 


K 


n 


^4..^ 


z2.-2 


>T 


(4.2.21) 


The   terms  £2(^22)  and  g2(E-k2)  are  not  considered  since  these  terms  give 
the  contri"bution8  to  the  inelastically  scattered  waves. 

The  excliange  scattering  terms  are  given  by  (4.2.21).  Both  equations 
give  the  asymptotic  part  of  f(k^k2);  (4.2.20)  for  the  singularities  on  the 
real  axis  of  k,  and  (4.2.21)  for  those  on  the  real  axis  of  k2.  Combining 
these  results  we  fina  that  the  elastic  scattering  amplitudes  for  the 


symmetrized  solutions  are 


izr^i' 


o 
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(U.2.22) 


„2     ,„   s2     lc,-K  ^ 

B  ^   N     /)/Ntt)   T. 


K^+B^   "^V^  /  (K^+T^)   T(T+B) 


tCih-e) 

^  0 


±  1 


Except  for  the  term  N  /(& +B  )  this  express  ion  is  an  approximation  to 
the  elastic  scattering  amplitude  of  the  exact  solution  in  which  K  /(K  -Bi) 
is  replaced  hy  unity  (valid  for  K  and  hence  E  large),  and  the  constants 
C  and  D  are  replaced  by  the  values  obtained  for  them  by  the  Born  approxi- 
mation in  the  symmetric  perturbation  formulation  of  this  problem. 

The  term  ir/(K^+B  )  which  appears  in  (4.2.22)  is  disconcerting.   It 
occurs  in  (it.2.2l)  snd  iii5>lies  that  exchange  scattering  takes  place  even 
if  the  interaction  between  the  particles  vanishes,  for  when  A  =  0, 
g^ (ifcp)  does  not  vanish.  On  physical  grounds  there  cannot  be  exchange 
scattering  when  the  particles  fail  to  interact. 

Equation  (4.2.21)  is  admittedly  the  res\ilt  of  using  an  approximation 
to  the  exect  wave  function.   If  we  assume  A  to  be  zero  and  iterate  using 
fdsLkg)  determined  by  (4.2.17),  (4.2.18)  and  f^dc^k^),  we  find  that  the 
n   iteration  gives 

(n)  f,^   .      /Bi\''    /    _i (4.2.23) 

0      0        c.      0 


and 


«2 


(a) 


'Va'-f^y-'^jV^- 


These  expressions  converge  to  zero  for  BiyK  |<1  and  Bi/|K2ll<l. 

Eor  real  values  of  K  and  K-  other  than  those  satisfying  these  inequalities 

these  dxpressions  diverge.   Ihus  the  physically  acceptable  solutions  vanish 

as  n  — >  00.  The  contributions  to  exchange  scattering  given  by 

i  TSr/iVT  +B  )  in  the  first  approximr^tion  is  spurious. 

The  exact  mathematical  calculations  justify  the  physical  requirements. 

PractiCF-lly,  however,  the  calculations  are  restricted  to  first,  second  or 
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possitly  third  order  epproximstions.   In  this  situation  the  contributions 
to  the  scattering  froci  the  spurious  terms  assume  en  unwarr&nted  importance. 
The  difficulty  arises  from  the  use  of  the  initial  state  as  e  first  approxi- 
mation to  the  exact  solution.   Hhe  physical  basis  for  taking  the  initial 
state  as  a  first  order  approximation  is  that  for  energies  siif ficiently 
large  it  is  expected  that  the  interaction  of  the  incident  particle  with 
the  atom  will  have  a  negligible  effect  on  the  initial  state  of  the  system. 
When  exchange  phenomena  are  possible  this  assijmption  is  incorrect.  Por 
exchange  phenomena  to  occur  two  conditions  must  be  met.  First,  the  incident 
particle  must  interact  with  the  center  of  force  so  thet  it  can  be  bound 
and  secondly,  it  must  interact  with  the  boTind  particle  so  that  an  energy 
transfer  is  possible.  When  approximation  methods  ere  to  be  used,  the 
relative  order  of  magnitude  of  these  interactions  is  also  important.  Naive 
considerations  lead  us  to  expect  that  the  interaction  between  the  incident 
and  bound  particles  is  smaller  on  the  average  than  the  interaction  between 
the  Incident  particle  and  the  center  of  force.   Consider  the  problem  in 
which  electrons  ere  incident  on  hydro£en  atoms.  The  distortion  of  the 
bound  electronic  charge  by  the  incident  electrons  wooild  oh  the  average 
result  in  the  distance  between  the  electrons  being  larger  than  that  between 
the  nucleus  and  incident  electrons.  Furthermore  the  fact  that  the  incident 
electron  can  be  bound  implies  that  the  attractive  forces  which  exist  would 
result  in  the  enhancement  of  the  distortion  of  the  electronic  cloud  about 
the  nucleus  and  thus  increase  the  influence  of  the  nucleus  on  the  incident 
electrons.  Hence  a  first  order  approximation  which  takes  into  account  the 
relative  magnitudes  of  the  interactions  involved  is  one  which  contains  the 
effect  of  the  interaction  of  the  center  of  force  on  the  incident  particles. 
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The  atove  physical  arguinent  can  "be  reinforced  "by  considering 
equations  (^.2.11)  and  (4.2.17).   In  (4.2.11)  replace  f(k^k2)  with  its 
expansion  (4.2.7).  Then 

{h.z.zh) 


f^(k^2)  =  6(k,-K^)  *  I  ^  fh^h^^^^l 


1     0 
'  t^    'i^^  (2lk2)*dk2(k^k2lvlk.k.)dk.dk.f  (k'k.). 

The  first  two  terms  on  the  right,  since  /  f^(k-2)dk-   is  e  constant, 

have  the  form  of  (K  [k^ )  which  is  a  solution  of  the  unperturhed  equation 

that  takes  into  accoxint  the  interaction  of  the  incident  particle  with  the 

center  of  force.  Now  take  (K  |k^)(2|k2)  as  the  first  order  approximation 

to  f (k^kg)  and  insert  it  into  (4.2,17).  We  obtain 

(^.2.25) 

g^dkg)  =  -  ^  -^  /(llk^)*dk^(k^k2lvlk.k^)dk«dk'(Kjk^)(2lkp. 
4Tr   k_-K^ 

4S    O 

Here  v;e  see  that  the  exchange  scattering  vanishes  when  the  interaction 
"between  the  Incident  and  'bound  particles  vanish.  The  result  of  using 
an  approximating  function  of  the  type  considered  has  been  to  eliminate 
the  spurious  terms  previously  obtained  in  (4.2.21)  and  to  make  g^(lk2)  of 
(4.2,25)  depend  on  the  interaction  between  the  incident  and  botind  particles 
only.  The  interaction  between  the  incident  particle  and  the  center  of 
force  has  been  incorporated  in  the  trial  function  used.   In  this  way  the 
order  of  magnitude  of  the  two  interactions  has  been  properly  treated  and 
the  weaker  interaction  is  now  the  perturbation. 

A  detailed  calculation  using  the  trial  function  (K  }k-)(2|k2)  gives 
the  following  asymptotic  part  for  the  elastic  scattering  coefficient  for 
the  symmetric  and  antisymmetric  wave  fxinctions. 


(2)-^/^ 
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(^.2.26) 
2.   ^0    -Cf^A^f   ^1    /T-B   1         1   \   ^o 


,"  K^-Bi    ^"^  '  (I^+T"   \  ^   K^+T*-    T(m-B)  /  Z  -Bi 


o 


o 


T(T+3)  V  K^^+T     'V. 


If  this  result  is  compared  with  (4,2.22),  which  was  obtained  "by  using 
the  initial  state  as  a  first  order  approximetion  to  the  exact  solution, 
we  find  thet  the  spurious  term  is  absent  and  that  the  factor  K  /Z  -Bi) 
appeexs  instead  of  its  approximc tion  by  unity.  Furthermore  the  constants 
C  and  D  are  much  better  approximated  by  (4.2.26)  than  by  (4.2.22). 
Equation  (4.2.26)  is  except  for  the  last  term  in  the  curly  brackets 
identical  to  the  result  obtained  by  the  Born  approximation  in  section  (3.?), 
This  is  significant  because  the  Bom  approximation  in  that  section  is  the 
leading  term  in  the  Bom  expansion  which  is  known  to  converge.  Thus  the 
Mott  and  Massey  procedure  with  trial  functions  which  take  into  account 
the  effect  of  the  interaction  between  the  incident  particles  and  the 
center  of  force  will  yield  results  comparable  to  those  obtained  by  treating 
the  problem  by  the  much  more  mathematically  difficult  symmetric  perturba- 
tion procedure. 

There  is  yet  another  argument  from  the  computer's  point  of  view  for 
the  use  of  the  altered  initial  state  trial  functions  discussed  above. 
Namely,  if  the  interaction  term  were  to  be  neglected  both  representations 
(4.2,11)  and  (4.2.16)  of  the  exact  solutions  must  give  the  initial,  state 
of  the  system.  This  is  easily  verified  to  be  so.  But  if  the  procedvixe 
leads  to  this  result  which  can  be  obtained  exactly,  surely  an  approximation 
method  should  avail  itself  of  this  fact  from  the  start. 

The  above  discussion  is  not  restricted  to  the  problem  under 
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consideration  "but  carries  over  directly  to  actw^l  physical  problems. 
The  physical  arguments  are  of  general  epplicatility  and  we  have  applied 
them  to  this  problem  so  that  they  could  be  illustrated  and  examined  in 
detpil. 

Table  YII  lists  the  numerical  results  obtained  for  the  elastic  and 
inelastic  scattering  cross  sections  from  the  Mott  and  Massey  procedure. 
The  first  two  columns  were  obtained  for  the  symmetric  solution,  and  the 
last  two  for  the  antisymmetric  solution.  The  column  headed  0       is  the 
elastic  scattering  cross  section  obtained  with  the  initial  stete  of  the 
system  taken  as  the  approximation  to  the  exact  solution  and  0  „  heads  the 
column  for  which  the  altered  initial  state,  i.e.,  the  state  which  takes 
into  account  the  effect  of  the  center  of  force,  was  used  as  an  approximation 
to  the  exact  ware  function. 

The  first  two  columns  give  exceedingly  poor  approximations  as  com- 
parison with  Table  IV  will  show  except  for  large  values  of  E.  This  is 
expected  since  for  the  symmetric  solution  we  have  already  seen  that  the 
Born  approximation  with  the  symmetric  perturbation  is  unsatisfactory  for 
all  but  laxge  values  of  E  and  the  results  here  are  approximations  to  this 
approximation. 

The  last  two  columns  when  compared  with  Table  7  indicates  that  0  „ 
is  a  much  better  approximation  than  0  g.  Here  the  comparison  is  much  more 
significant  tht-n  the  previous  one.  For  the  Born  approximation  with  the 
symmetric  perturbation  for  the  antisymmetric  wave  function  yields  excellent 
values  for  the  scattering  cross  sections  as  Table  V  indicates.  As  a 
matter  of  fact  the  Born  approximation  is  part  of  the  Born  expansion  and 
this  is  well  within  its  radius  of  convergence  as  is  rlso  shewn  in  Table  V. 
The  0       and  0  _  of  Table  VII  ere  approximations  to  the  0   of  Table  V  as 
was  shown  in  equations  (^.2.22)  and  (^.2.26)  and  the  fact  that  Q^j^  is 
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Table  VII 


Symmetric  Solution 


Antisymmetric  Solution 


1^ 


m. 


!fi£- 


!&1L 


-.9 

-.5 
-.1 

-.05 

-.01 

-.005 

-.001 
+  .001 

+.005 

+  .01 

+.05 
+.1 
+.5 
+.9 

5. 

10. 

100. 

1000. 


80.6882 
20.3238 
11.91^5 
11.3275 
10.9017 
10.3451 
10.8040 
10.7836 
10.7493 
10.6927 
10.3053 
9.85636 
7.23167 
5.60363 
I.I73&6 
.429427 
10-1.138501 
IO-2.IO38I6 


155.576 
76.0175 
43.5778 
41.1821 
39.3207 
38.8173 
38.6335 
38.5421 

38.3599 

38.1334 

36.3899 

34.3628 

22.6594 

15.8669 

1.99032 

.560365 

10-1.143099 

10~2. 104165 


7.14250 
.414772 
10-1,767417 
10-J-.651345 
10"*!. 577614 
lo-i. 569441 
10*1.561920 
10"  ■'•.559906 
IO-I.5537I8 
10-1.546157 
10*1.492562 
10-1.440770 
10-1.362542 
10*1.509812 
.141181 
.119278 
10*1 . 123460 
10*^.102935 


8.84130 

3.05059 
1. 33780 
1.22707 
1.14778 
1.13846 
1.13105 

1.12737 
1.12005 
1.11102 
1.04245 
.965360 
.579865 
.409900 
.184651 
.125789 
10-1.114606 

10*2.101717 
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superior  to  0  _  confirms  vhat  has  "been  said  about  the  altered  initial 
state  as  an  approximation  to  the  exact  wave  function. 

As  in  previous  considerations,  it  is  desirable  to  study  the  Schwinger 
variational  principle. 

it. 3  Ihe  Schwinger  Variationel  Principle— Ecrowitz-Friedman  Procediire 

The  solution  of  the  Schroedinger  eq-uation  (3.1.2)  is  written  es  a 
sum  of  the  incident  wave  end  en   outgoing  wave.  Thus 

f  (k^k2)  =  ^0(^1^2^  *  ^(^^2^  (^-3.1) 

where  f  (k-k^)  is  the  solution  of  the  unperturbed  equation  and  is  given 

by  (ii.2.19).  The  function  h(k  k^)  represents  the  outgoing  particles  1 

end  2.  The  Borowitz-Friedman  (see  footnote  10)  procedure  constructs  en 

integral  equation  for  h(k^k„)  from  which  the  exchange  scattering  ajaplitudes 

may  be  computed  by  means  of  the  Schwinger  variational  principle.   This 

procedure  can  be  formulated  so  that  the  direct  as  well  as  the  exchange 

scattering  amplitudes  ean  be  treated  symmetrically. 

Substitute  (^.3.1)  into  the  Schroedinger  equation  (3.1.2).  Then 

h(kLk  )  must  satisfy 

(^.3.2) 

(k^+k^-E  )h  (k^kg )-  (B/tf )  y* h  (k^k^  )dk^-  (b/tt)  f  h  (k^^k^  )dk2 

-(B/n)/f^(k^k2)dk^+(l/^TT^)/(k^lc2l^l^^2''^^2  [^o^'M^^+^^M^^  |  =  ° 
This  equation  may  be  written  in  the  following  two  weys 

(^.3.3) 

(kJ+k2-Z)h(k^k2)-(B/n)  /h(k^k2)dk^  =/w(k^k2;k|k^)dk|dk^f  ^(k|k|) 

+y  V  (k^kg :  k^kpdk|dk^  (k|k| ) , 
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or  (^.3.^) 

(k^+k^-E)h(k^k2)-(B/TT)yh(k^k2)(ik2  =y  w(k^fc2;k^kpdk^dk|f  ^(k|k|) 

+  y  w  (k^k2;k^k^)dk^dk^  (k^kp . 

in  which 

w(k^k2;k^kp  =      (B/TT)6(k2-kp  -    (l/^/)  (k^k^l?  |k|kp  (4.3.5) 

and 

v(k^k2;k^kp  =  (B/TT)6(k^-k^)  -  (l/i+TT^)(kj^k2|v|k|kp.  (4.3.6) 

Note   that 

w(k^k2;ki^kp  =  v(kjk^ik^p.  (4.3.7) 

Denote   the  Green's  f\mctions   of   the   left  hand  sides   of    (4.3.3)  and    (4.3.4) 
hy  G        (k^k2;kJki)  end  G(k^k2;k'kA)  respectively.     The  relation  between 
these  Green's   f\inctions   is 

G^^^(k^k2;k|kp  =  G(k2k^ik|k^).  (4.3.8) 

This  follows  from  the  fact  that  the  left  hand  sides  of  eq-oations  (4.3.3) 
end  (4.3.4)  differ  only  in  an  interchenge  of  k^  and  kg.  The  explicit 
expression  for  6(k,k2;k'k»)  is  given  by  (2.2.37)  or  (2.2.38). 

V/e  may  now  represent  h(k^k2)  by  either  of  the  following  two  integral 

equations : 

(4.3.9) 

h(k^k2)  =  0C(k^k2)  +yG^^^k^k2;k|k^;dk^dk|v(k|k^;k^«'k2")dk^"dk2"h(k^«k2") 

or  (4.3.10) 

h(kj^k2)  =  ?i^^2^   +yG(k^k2;k^kpdk^dkj^w(k^k^;k'^«k«|)dld^«dk«^(fc!^«kip, 
where  (4.3.11) 

o((k^k2)  =yG^^\k^k2;fc[k^)dk|dk^w(k|k^ik^"k2")dk^"dk2"f^(k^»'k2"). 
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"^  (4.3.12) 

Pdc^kg)  =y  G(k^k2;k^kpdlc^dic^v/(k^lc|;k^»':c2»»  )dk^»«dk2"f  ^(k^'^k^"). 

The  integral  eqiiation  v;hich  represents   the   solution  of    (3.1»2)   is 

(4.3.13) 
f(k^k2)  =  fgd^kg)  +yG(k^k2;k|kpdk^dfc|w(k|k^;k^"k2")dt|"<ik2«'f(k^««k2")    . 

When  f  (^,^2)  is  substituted  for  f (k-kg)  in  the  above  integral  then  the 
resulting  integral  is  ^(iz^t^).     Hence  p(k^k2)  is  the  correction  to  the 
initial  state  that  is  obtained  froni  the  Born  approxiaiation.  The  expression 
0((k^k2)  cannot  be  similfLTly  interpreted. 

The  integral  equations  (4,3.9)  and  (4.3.10)  have  the  property  of 
explicitly  expressing  the  asymptotic  part  of  h(k^k2)  in  the  momenta  k2 
and  k-  eespectively.  This  is  due  to  the  fact  that  the  respective  Green'o 
functions  ere  explicitly  given  at  the  singularities  on  the  real  axis  of 
ko  and  k,  respectively. 

Equation  (4.3.9)  which  expresses  the  asymptotic  part  of  h(k-k2)  on 
the  real  axis  of  kg  explicitly  does  not  vanish  with  the  interaction  between 
the  incident  and  bound  particles.  This  is  as  it  should  be  for  the  eq^uation 
represents  a  solution  of  the  problem  and  when  the  interaction  vanishes  there 
are  still  outgoing  waves,  namely,  those  waves  that  result  from  the  inter- 
action between  the  incident  waves  sind  the  center  of  force.   It  is  easily 
verified  that  with  the  interaction  between  the  incident  and  bound  particles 

set  equal  to  zero  that 

(4.3.14) 
h(k,k  )  =  ^    0    N     1 

^^^    "   Z-Bi  k^.B^  k?-K^   • 
0       2       1   0 

With  this  expression  inserted  in  (4.3.1)  we  indeed  obtain  the  solution 
of  the  unperturbed  problem.  This  result  implies  that  h(k^k2)  nnist  not 
be  set  equal  to  zero  on  the  right  side  of  (4.3.9),  for  then  in  this 
approximation  with  the  interaction  between  incident  and  bound  particles 

set  equal  to  zero,  v/s  obtain 
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R     ^^^2  ^         _J5 1 

which  has  for  its  asymptotic  part 


hCk^fcg^^ 


jf_   N     1   ^  _B   '^1     N 


K^+B^  k^+B^  k|-^     TT  |k^!-Bi  kg+B^  ^i+'^a"'- 

This  last  expression  indicates  that  there  is  exchange  scattering  even 
tho\]gh  physically  we  know  this  to  "be  impossihle.  Bie  contribution  to 
exchange  scattering  is  just  the  spurioufl  term  found  in  the  preceding 
section.  Secondly,  trial  functions  for  hCk-kp)  must  be  chosen  such 
that  it  consists  of  the  term  (4.3.11)  plus  a  function  that  depends  on 
the  interaction  between  the  incident  and  boimd  particles.  Such  a  choice 
not  only  describes  what  is  occurring  physically  but  also  insures  that 
exchange  scattering  will  vanish  when  the  intersiction  does.  The  need  to 
choose  h(k^k2)  in  this  way  is  a  confirmation  of  the  analysis  of  this 
problem  given  in  the  preceding  section. 

The  asymptotic  part  of  p(k^k2)  which  is  independent  of  the  inter- 
action between  the  incident  and  bound  particles  is 

B     H        1 

"  k^  +  b2   k2   ^   • 
^         X    o 

Excdpt  for  a  constant  factor,  this  expression  is  the  same  as  (h.^.l^). 

Dius  $(k  k-)  has  just  the  form  necessary  to  serve  as  a  trial  function 

for  hCk^kp).  Calculations  have  been  carried  out  using  (k.J.lk)   as  an 

approximation  to  the  exact  solution.  The  results  are  identical  to  (4.2.26) 

which  was  obtained  from  the  Mott  and  Massey  procedure  with  the  altered 

initial  state  as  the  trial  function.  This  is  to  be  expected.  (See  footnote 

10).  But  now  the  integral  equations  are  in  such  a  form  that  we  can 
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apply  the  Schwinger  variational  principle  to  them. 

The  integral  equations  (^.3»9)  and  {>^,J.10)   may  he  written 

aCs^kg)  =  K  hCt«k^),  (^.3.15) 

and 

aiCk^kg)  =  L  h(k^k|)  (4.3.16) 

respectively,  where 

aCk^kg)  =yv(k^k2;fc{kpdk^dk:]o^(k^k|)  (4.3.1?) 

a«(k^k2)  =y  w(k^k2ik^kpdk^dk^  P(44)»  (4.3.18) 

and  K  and  L  are  integral  operators  with  kernels 

(4.3.19) 

ECk^k^rk^kp  =y*  r6(k^-k^««')6(k2-k*2»')-v(k^k2jk|«fc^')ikWdk|«G^^^k^«k«^««k«t  ) 

.  dk«««dk*2"v  (k[««k^' '  ;k|kp , 
and  (4.3.20) 

L(k^k2;k^kp  =/  r6(k^-k?|»')6(k2-k2'")-w(k^k2;k^««k2")dk^"dk«2«  G^k^Hkg"  ikj»«k*2"  )  J 

•  dk}"dk««w(k"'k«2";k^kp. 

Each  of  the  equptions  (4.3.15)  end  (4.3.16)  can  he  treated  in  a  ms-nner 
completely  annlogoMS  to  the  development  of  similar  eq-oations  in  section 
(2.6).   This  variational  principle  has  heen  exhaustively  discussed  in 
Chapters  II  and  III. We  shall  not  repeat  the  details  here  except  to  say  that  the 
computations  have  heen  carried  out  only  for  the  determination  of  the  direct 
end  exchange  scattering  amplitudes.  From  these  we  ohtained  the  scattering 
cross  sections  hy  pursuing  the  procedure  used  in  calculating  the  exact 
cross  sections  in  section  (2.3). 
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Tatle  VIII 


SCHWINGEE  VAEIATIOEAL  PRINCIPLE 
BOROWITZ-PEIEDMAN 


Symmetric 


Antisymmetric 


-.9 

89.8027 

-.5 

19.6501 

-.1 

12.0166 

-.05 

11.^933 

-.01 

11.1119 

-.005 

11.0686 

-.001 

11.032it 

+  .001 

11.0164 

♦  .005 

10.9796 

+  .01 

10.9^20 

+.05 

10.6158 

+.1 

10.2279 

+.5 

7.91153 

♦.9 

6.40569 

5. 

1.36008 

10. 

.467606 

100, 

10-1.137221 

luOO. 

10-2.104397 

6.83310 
,.195340 
10- J-.  197900 
lO"-'-.  204200 
10-^.206700 
10",. 208400 
10"-. 208900 
10-1.209600 
10"1.214300 
10-1.214300 

lo-i. 302700 
io~i. 263500 

10-1.717500 

.124870 

.205875 

.151366 

icrl  126091 

10-". 102327 
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Table  VIII  contains  the  resiilts  of  the  numerical  calculations.  Here 
again  we  find  that  for  the  Born  trial  function  the  results  are  far  from 
satisfactory.  This  failure  is  due  to  the  fact  that  in  this  case  (2.6.26) 
holds  for  the  symmetrical  and  antisymmetrical  solutions  only  for  large  E. 
The  variationfil  principle  does  not  extend  the  usefulness  of  the  results 
to  smaller  values  of  E.   This  difficulty  has  been  met  with  before. 

To  recapitulate:   The  application  of  the  iterative  procedure— the 
Born  expansion— leads  to  results  which  differ  from  those  obtained  in 
Chapters  II  and  III.   The  series  representation  of  the  solution  does  not 
converge  to  a  solution  which  accounts  for  exchange  scattering.  Therefore 
this  iterative  procedure  cannot  yield  any  information  about  exchange 
scattering.  The  Mott  and  Massey  procedure  does,  however,  give  exch^joge 
scattering.  But  if  the  initial  state  is  used  as  a  first  order  approximation 
to  the  solution,  spurious  terms  are  obtained.   If  the  initial  state  is 
altered  by  including  the  effect  of  the  interaction  of  the  center  of  force 
on  the  incident  particles  (the  altered  initial  state  is  to  consist  of 
incident  waves  plus  outgoing  spherical  waves  and  is  a  solution  of  the  un- 
perturbed wave  eaurtion  in  which  the  interaction  between  incident  and  bound 
particles  is  set  equal  to  zero),  then  the  spurious  terms  do  not  appear  and 
in  addition  the  approximation  to  exchange  scattering  is  improved.  These 
results  are  directly  applicable  to  actual  physical  problems.   'Hie  Borowitz- 
Friedman  approach  to  this  problem  is  to  constnict  integral  eo^uations  for 
the  scattered  wave  only.   In  this  formulation  an  iterative  procedijire  is 
possible  for  the  calculation  of  the  exchanged  scattered  wave.   This  is  not 
true  of  the  Mott  and  Massey  (see  footnote  10)  procedure.  Moreover  the 
analysis  of  this  proced-ore  confirms  the  advisability  of  using  the  altered 
initial  state  as  a  trial  function.   In  addition  this  method  can  make  use 
of  the  Schwinger  variational  principle  for  the  calculation  of  the  exchange 

scattered  amplitudes  which  is  not  a  property  of  the  Mott  and  Massey 


-  Ill  - 

formulation  of  the  problem.  The  application  of  the  Schwinger  variational 
principle  leads  to  results  which  have  already  heen  noted  in  the  preceding 
chapters. 


Chapter  7 
CONCLUDING  EEMAEKS 

In  th.e  preceding  chapters,  we  have  considered  the  Born  approzimation 
with  the  associated  Born  expansion,  the  Schwinger  and  Hulthen-Kohn  varia- 
tional principles,  and  where  occasions  required  it,  some  variations  of 
these  approximation  methods  such  as  the  Mott  and  Massey  procedure  and  the 
Bo row itz- Friedman  extension  of  the  Schwinger  variational  principle.  These 
methods  have  been  applied  to  specific  problems  with  defined  interactions 
and  the  conclusions  which  were  reached  must  now  be  examined  in  the  light 
of  their  general  applicability. 

Specific  problems  serve  at  least  three  important  purposes.  First, 
they  may  be  counter  examples  to  conjectures,  to  seemingly  plausible 
assumptions  or  to  conclusions  which  have  been  extrapolated  from  known 
results.  Secondly,  they  can  suggest  ways  and  means  of  obtaining  general 
results.  Finally,  they  may  point  out  certain  aspects  of  general  significance 
which  are  worth  further  investigation.  The  problems  studied  here  have  been 
fruitful  in  these  respects. 

Consider  the  Bom  expansion  which  incorporates  the  Bom  approximation 
and  which  was  studied  in  Chapters  II  and  III.  It  was  shown  that  the  Born 
expansion  converged  to  the  exact  solutions.  In  the  first  of  the  problems, 
dealt  with  in  Chapter  II,  the  Born  approximation  does  not  have  any  useful 
domain  of  applicability  whatever  the  value  of  the  total  energy  unless  the 
interaction  between  the  incident  and  bound  particles  is  vex'y  weak.  Here  the 
meaning  of  weak  is  related  to  the  shortness  of  the  range  of  the  interaction. 
That  the  Bom  approximation  may  not  have  a  usefiol  domain  of  applicability 
is  not  a  consequence  solely  of  the  interaction  employed  in  this  problem. 
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Jo3t  and  Pais  (see  footnote  5)  obtain  a  similar  result  in  an  analysis  of 
the  Born  approximation  for  a  Yiilcawa  potential.   Furthermore  there  is  an 
additional  result  that  their  problem  has  in  common  with  the  one  studied 
here  which  is  that  the  second  Bom  approximation  furnishes  a  substantial 
correction  to  the  Born  approximation  even  for  extremely  large  values  of 
the  energy.  For  the  problem  dealt  with  here  this  substantial  correction 
defines  a  usefxil  domain  of  applicability  for  the  results  so  obtained. 
Moreover,  this  example  indicates  that  in  general  the  second  Born  approxima- 
tion must  be  calculated  if  useful  results  are  to  be  obtained.  Physically 
the  need  for  such  additional  calculation  can  be  understood  qualitatively 
by  considering  the  Born  expansion  of  the  exact  solution  which  may  be 
represented  in  general  by 

^  =  ^0  ^  f;  ^n  =  ^0  "  ^'  ^^''-^ 

n=l 

where  f  is  the  exact  solution,  f  the  initial  state  of  the  system  and 

0 

g  =  /  f   represents  the  scattered  wave.  Here  the  Bom  approximation  is 

given  by  f  +  f  ^  and  the  second  Born  approximation  by  f  +  f -,  +f  o  • 

The  initial  state  of  the  system,  f  ,  does  not  denend  on  the  interaction 

0 

which  is  taken  as  the  pert\irbation.   However,  the  scattered  wave,  g,  does. 
Calculating  the  Born  approximation,  f  +  f,  yields  the  first  term  in  the 
expansion  for  g  which  of  course  contains  no  information  about  the  usef\il— 
ness  of  this  expansion  for  g  though  it  does  furnish  a  measure  of  the  effect 
of  the  perturbation  on  the  initial  state.  But  for  the  Born  approximation 

to  be  of  use  the  effect  of  the  pertTirbation  on  the  successive  terms  of  the 
expansion  for  g  must  be  determined  and  this  involves  at  least  a  calculation 
of  the  second  Bom  approximation  ajid  a  comparison  of  this  with  the  Bom 
approximation.   Dalitz  (see  footnote  5)  in  a  study  of  the  Yutewa  potential 
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in  the  limit  as  it  approaches  the  Coulomb  potential  also  found  it  necessary 
to  go  "beyond  the  Bom  approximation.   Actually  he  calculated  terms  of  the 
third  order  in  the  Born  expansion  although  the  second  order  terms  already 
furnish  information  ahout  the  convergence  of  the  series. 

The  calcuJation  of  higher  order  terms  in  the  Born  expansion  is  difficult. 
Thus  we  have  considered  the  Schwinger  variational  principle  restricted  in 
application  to  Born  trial  fields  "because  otherwise  the  mathematical  compu- 
tations become  extremely  complicated.   However,  for  such  trial  f\mctions 
we  find  that  the  relation  of  the  second  Born  to  Born  approximation  discussed 
above  again  plays  a  significant  role.  la  Chapter  II  it  was  shown  that  if 

<1  (5.2) 


(f  , ,  VG7f  )/  (f  , ,Vf  ) 
^  oi'     0/    oi*   o' 


then  the  variational  results  include  the  second  Bom  approximation.   The 
numerator  of  this  ratio  is  the  asymptotic  part  of  f^,  and  the  denominator 
is  the  asymptotic  part  of  f  both  of  which  appear  in  the  expansion  (^.1), 
For  the  problem  considered  in  Chapter  II,  the  second  Born  approximation 
is  useful  for  all  values  of  the  interaction  strength  A.  However,  in  this 
case  (5»2)  is  not  satisfied  \inless  the  values  of  A  are  restricted.  Thus 
we  see  that  the  variational  results  do  not  include  the  second  Bom  approxi- 
mation which  is  known  to  be  good  unless  the  interaction  is  weak  and  weak 
here  means  that  the  Bom  approximation  is  good.   Obviously  under  these 
conditions  there  is  no  point  in  using  a  variational  principle.  In  Chapter 
III  the  sit-uation  is  different.  Both  the  variational  and  Bom  results  for 
the  antisymmetric  solution  are  good  approximations  to  the  exact  values  for 
large  values  of  the  energy.  But  as  the  energy  becomes  smaller  the  varia- 
tional results  get  poorer  whereas  the  Born  results  are  good  throughout  the 
energy  range.  Here  it  is  clear  that  the  variational  principle  fails  to 
extend  the  Bom  results  beyond  their  common  doiuain  of  agreement.  Although 
the  variational  and  Bom  results  may  not  be  good  approximations  to  the 
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exact  values  for  large  values  of  the  energy  yet,  as  an  example  in  Chapter 
II  reveals,  the  variational  results  may  be  in  good  agreement  with  the 
exact  values  for  small  energies  where  the  Born  results  are  inadeauate. 
From  these  resiolts  enumerated  above,  it  is  clear  that  the  Schwinger  varia- 
tional principle  restricted  in  application  to  Born  trial  functions  does  not 
necessarily  yield  results  superior  to  the  Born  or  second  Born  approximation. 
The  basic  problem  here  is  the  determination  of  some  guide  for  the  selection 
of  appropriate  trial  functions.   We  vill  return  to  this  point  after  a  con- 
sideration of  the  Hulth^n-Kohn  variational  procedure. 

This  variational  procedure  is  more  amenable  to  mathematical  calcula- 
tion than  the  Schwinger  variational  principle  because  it  is  not  formulated 
in  terms  of  a  G-reen's  function.   For  comparison  puposes  its  application  in 
the  preceding  chapter  was  restricted  to  Born  trial  fxinctions  but  then  the 
results  so  obtained  were  shown  to  be  equivalent  to  the  Born  results. 
Though  the  mathematical  simplicity  of  the  Eulth^n-liohn  variational  principle 
was  not  fully  exploited  yet  a  general  relation  betv/een  this  principle  and 
Schwinger' s  was  established.  The  significance  of  this  relationship  is  due 
to  the  fact  that  Kato  (see  footnote  9).  for  two  body  problems,  has  been  able 
to  establish  upper  and  lower  bounds  for  the  Schv^inger  and  Hulth^n-Kohn 
variational  principles.   This  was  accomplished  by  suitably  restricting  the 
class  of  trial  functions  to  be  used.  The  general  relationship  between  these 
principles  and  the  relative  mathematical  simplicity  of  the  Hulthfen-Kohn 
variatioiia.1  principle  indicates  the  possibility  of  extending  Kato's  result 
to  many  body  problems  in  terms  of  trial  functions  to  wliich  the  Eulthfen-Kohn 
variational  principle  may  be  directly  applied. 

Chapter  IV  dealt  with  the  problem  of  exchange  scattering.   The  formu- 
lation of  the  symmetric  problem  using  a  symmetric  pertiorbation  led  to  the 
result  that  each  term  in  the  Born  expansion  gave  evidence  of  exchange 
scattering  and  in  particular  the  Born  approximation  takes  this  effect  into 
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account.  This  same  problem  using  an  asymmetric  perturtation  leads  to 
difficulties  in  the  calculation  of  exchange  scattering.   The  Bom  expansion, 
in  this  case,  gave  rise  to  a  series  representation  of  the  exact  solution 
which  did  not  contain  exchajige  scattering.  This  is  due  to  the  fact  that 
the  series  is  not  valid  at  the  singularity  of  the  wave  function  which 
yields  the  exchange  scattered  wave.   This  behavior  illustrates  and  supports 
a  conjecture  due  to  Kohn  (see  footnote  4o)  in  which  it  is  stated  that  the 
Born  expansion  does  not  converge  to  a  solution  which  accounts  for  exchange 
scattering.  The  latter  physical  process,  of  course,  mast  he  possible  for 
the  problem.  Though  this  conjecture  has  been  verified  only  for  the  specific 
problem  considered  here,  yet  the  result  throws  considerable  doubt  on  the 
applicability  of  this  iterative  procedxire  to  a  determination  of  exchange 
scattering  in  actual  physical  problems.  The  Mott  and  Massey  procedure 
which  was  developed  to  handle  this  problem,  in  turn,  gave  rise  to  difficul- 
ties which  have  only  recently  fsee  footnotes  40,4l)  been  cleared  up.   In 
Chs-pter  IV  this  procedure  is  examined  when  the  Born  approximating  procedure 
is  used,  i.e.,  the  initial  state  is  used  as  a  first  order  approximation  to 
the  exact  function.   The  result  of  this  study  was  that  this  method  does 
account  for  exchange  scattering  but  that  sptirious  terms  arise  in  the  calci*- 
lation  of  the  exchange  scattered  wave.  These  spurious  terms  assume  unwarranted 
importance  especially  if  the  calculated  approximation  is  used  in  an 
iterative  procedxire.  However,  if  the  initial  state  is  altered  so  that  the 
effect  of  the  interaction  between  the  incident  particles  and  the  center  of 
force  is  taken  into  accoxint  (describing  this  in  coordinate  space,  the 
altered  initial  siate  is  the  incident  wave  plus  outgoing  spherical  waves) , 
then  no  stiurious  terms  arise  and  the  approximation  so  obtained  is  comparable 
to  the  Born  approximation  for  the  problem  formulated  with  a  symmetric  pertur- 
bation. These  results  are  of  general  validity  and  are  not  restricted  to  the 

particular  problems  considered  here. 
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APPENDIX   I 
The  explicit  expressions   for  the  qxientities  appearing  in  equations 
(2.2.16),    (2.2.17),  and   (2.2.18)  are 


r    ,    f- ^^1 , 

■"•       V       (k- +iT)(P-Bi)(m-iT) 


(lc^+iT)(P-Bi)(lH-iT) 


P  =    (E  -  ic^)^/^ 


i 


^    inB  -  22i  (tt  +  i  i'AuS(B)    )  -  i(T^+E)^/^(TT  -iAs(T)    ) 


K^+T^ 


-     0 


E 


] 


+  -4—   J  ttT  +        }^   .  y.      (tt  -  iJ?^S(T)    ) 


2T^+E 


(T^+E)^/^ 


=  / 


^ 


P(P-Bi)(PfiT)' 


(BfT)(ir+E)  K^(BfT)  (B+T)(T+E)   '    ISfT       T  +E 


^      P(k?+T 


^F-2-     =  2       1/2      (--iAs(T)    ). 


=/ 


dk. 


(k^+iT)(P-Ei) 


j^^^2 


>T 


7TB  -  2|i  (tt+   i£n.S(3))  -  i(T^+E)^/^(TT  -i^S(T)) 


S(x)=    Hx2.E)l/^.xl/("(x2.E)^/2-x 


In  the  above  expressions  E  >  0.     For  E  <  0  the  logPTitiimic  terms  are 
defined  by 

.feU)  =    ^((x^+E)^/2  +  x}    -A((x^+E)^/^  -  xj 

where  for  each  term  the  principal  vElijee  of  the  logarithms  are  to  be  tpken. 
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For  very  large  positive  values  of  E,  we  find  that 


J^   w^  tt(B  -  2T)  (E)"^/^  -  ttCE)"-"-! 

Jp  '^  (V3)  (B  -2T)E"^  -  (TT/2)E~^/^i 


Jo  rV 


(tt/t)E~^/^  -  2(E)"^i 


J^^  ,^     (ttB  -2T)E~-'-  -  TtE"-'-/^! 
For  large  values   of  E,   the  exact  constants  C  and  D  become 

C  '^  NTTi/B(Bfl) 

D^   (ntt/bCK^  -iT)_y  1^1  +    (A/^Tt)    (^7/T(^^-T))  | 
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APPENDIX  II 
Using  Born  trirl  functions  f  ,  =  (K  2lk:,k„);(ir  =  B  +E),  end 

01  O  1    <i  o 

f       =    (K_Kp;lc,k:2)  I    (KT+KZ  =  E),      in  the   Schwinj^er  veriatione.1  principle 
of  section   (2.6)   yields   the  following  values  for  the  constnnts  C  and  D. 

C    =     bTTi/B(B+T)       1    -     (A/^TT)j^i    +     U/^){BJ^+iJ^)(>    +     (A/4TT)(J^/T)h]    A'"'' 

D  =    [NTT/B(K^-iT)][l  -    (A/Mj^i  +    (A/4tt)  (BJ^+iJ-j)  €  +    (A/itrr)  (J^/T)t  |   ^^ 

^=    fl  -.    (A/4TT)j^i  +    (A/4Tr)(BJ2+iJ3)(»^T)(2Tr^-(A/iiT7)(j^/T)(K^+T^)(2Tr^(E+Tr^J 
Tl  -    (A/i+TT)j^i  -    (A/iiTT)(BJ2+iJ3)(B4-T)(K^-iTr^(p/i<X)+(A/iHT)(j^/T)(v/o<)j 


where 


OL    = 


iKgl  i(BfT) 


Bi 


^• 


+iT        iK^l+Bi        K^-iT     Ug+iT  (lK2l+Bi)(|K2l-iT) 


P     =     - 


.K2I 


K|_+iT        iKgl+Bi 


V     = 


V^iT 


Bl 


(|KJ+Bi)(lKJ 


-iT)    / 


£     = 


h     = 


B-l-T 

K  -iT 
0 


B+T 


2T  I 


(BfT)  _o 

0<  "     2T 
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e  = 


K  -IT    I  2T 
0 


£  =  - 


K  -IT 

0 


/  K  -IT      .  a  \ 
2T    / 


(^* 


4.4  = 


E  and  K^  =  B^  +  E. 


Three  limiting  cases  may  be  immediately  evaluated. 
Case  1  —    KI'^  E,  £,  small  and  E  large. 
Then     C  <^   NTri/B(BfT)  and 

D  ^  Ntt/B(K  -iT). 
These  are  the  Born  results. 
Case  2  —   ^'^  ^»  ^2  ^°^^-'-  ®^^  -^  large  . 


Then 


C  ^^ 


Nrri 


and 


B  r^ 


B(BfT) 


B(K  -iT) 
0   ' 


^       k-n       T 


(B±B)-T 


T^+(BfT)  {_(BtB)-T}^ 


-1 


■'■  "  Jm   T   „2. 


T^+(BfT)  jjBfcB)  -T^ 


J 


-1 


The  plus  sign  in  (B±B)  is  to  he  associated  with  positive  vslues  of  Zg 
and  the  minus  sign  with  negative  values  of  K^. 

These  results  do  not  reduce  to  the  Born  results  unless  T  =  2B  and 
Kj,  is  taken  positive. 

Case  3  —  KT  =  KZ  /v/  e/2  ;  E  large.  The  results  for  C  and  D  are 


those  of  the  first  crse  considered. 
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APEBHDIX  III 

An  Altern&tive  Proof  of  the  Relation  Between  the  Hulth'^n-Kohn  and 
Schwinger  Variations  1  Principles 

•   I^^i»^^  =  1^   ^i(E').  [h-e]t(E)J  (1) 

H  in  (l)  represents  the  Hamiltonian  operator  associated  with  the  problem 
under  consideration.   The  Hajniltonieji  H  may  be  resolved  in  the  following 

way: 

H  =  H  -  3  +  T 
0 

Here  the  operator  H  -  E  may  be  considered  to  represent  the  system  when 
the  component  parts  do  not  interact,  i.e.,  in  the  region  in  which  V  is 
negligible.  The  Schroedinger  equation  is  then 

(H  -  E  +  V)  f  (E)  =  0  (2) 

or 

(H^  -  E)  f  (E)  =  -  Tf  Cb).  (3) 

Equation  (3)  may  be  solved  symbolically  and  the  solution  is 

f  (E)  =  f^(B)  -  ^H^  -  E  -  ioi^  "^  Vf  (2),  (4) 

where  f  (E)  satisfies 
0^  ' 

(H  -  S)  f  (E)  =  0 

0         0 

and  satisfies  the  initial  conditions  of  the  problem.  Equation  (k) 
represents  the  solution  f (E)  as  an  incident  wave  plus  outgoing  waves. 
Another  solution  of  (3)  which  represents  an  incident  wave  plus  incoming 
waves  is  given  by 

f^CE)  =  f^^(E)  -  (H^  -  E  +  i<^r^  V  f^(B),  (5) 

where  f  .(E)  Is  any  solution  of 

(H^  -  E)  f^^(E)  =  0.  (6) 
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Multiply  eqiiations    (h)   and   (5)  by  V,    then  these  equetions  mey  "be  expressed 

as 

Tf    (E)   =  Vf  (B)  +  V  {h^  -  E  -  kI""''  7f  (E)  (?) 

=    (V  +  T  [h     -  E  -  io^H  V)f  (S)  =  K  f  (E) 

and 

T  (8) 

Vf^^(E)  =    (V  +  V  {h^  -  E  +  iaC]f"     V)  f^(E)  =  K«f^(E). 

IPhe  definitions  of  the  operators  K  end  Z*  iEtaedietely  lead  to  the 
relation  that 

U,   Ky)  =   (Z'x.y)  (9) 

i.e.,  that  Z'  is  the  hermitian  adjoint  of  K.   Therefore  from  (?)  and  (8) 

(f^(E).  Vf^E)  )  =  (f^(E),  Zf(E)  )  =  (K'fj^(E),  f(E))      (lO) 

=  (VfQi(E),  f(E)  ). 
Take  the  following  trial  fimction  for  v(B) 

r(3)  =  f^CB)  -  [(foi'^^^/^^oi'^S^jl^^o  -  ^  -  i<^^"^«]         (11) 
where  g  is  an  arbitrary  function.  When  g  is  equal  to  f  then  v(E)  =  f (E) 
from  (k).     Substitute  (ll)  into  (l). 

(12) 

^{^i»^}  =  ^]^  (i^^')'  {v^+tI  l/o^^)-  Ki'^^M^oi'^^^lJ  v^^'^^'^^li) 

=  lim   ^^(E»),  Yf^(B)-  ((f^^.Vf)^(f^^,Vg)|j^Vgf7(H^-B-i(0~V3  | 


=   ^v^(3),  Vf^(E)  -  {(foi»^^)/(foi'^«)]  ^) 
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where  the  definition  of  K  of  equation  (?)  has  been  uBed.  The  limit 
may  be  taken  in  (12)  because  the  potential  function  is  a  well  behaved 
function  and  the  limiting  process  may  be  interchanged  with  the  operation 
defined  by  the  inner  product.  Then,  omitting  the  dependence  of  the 
quantities  on  E, 

(13) 
I  ^^v^.  v]  =  I'  [v^,sj     =  (v^.Vf^)  -(f^^.Vf)(v^.Vf^)  x{t^.6)^ 

=  (v,.Vf^)  [l  -  (f,,.7f)X  \v^,e}] 


where 


{v^.sj  =   (T^.K€)/(v^.Vf^)(7f^^,g)  (14) 


Equation  (l4)  is  the  expression  associated  with  the  Schwinger  variational 
principle. 

Consider  the  variation  of  (13)  about  the  solution  of  (4)  and  (5). 

(15) 


'v=f         V     ^g=i 

^i=^i  ^i=^i 


-(f^^.Vf)x(v^4(6v^,Vf^)] 


=  -(f.i.Vf)^  6X  [v^.g] 


For  from  (lO)  and  (l^J), 

(f^.Vf^)  =  (Vf^j.f)  and  ^ff^.f}  =  l/CVf^^.f)  respectively. 
Equations  (13)  and  (l5)  establish  the  relations  sovight. 
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